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DIELECTRIC BREAKDOWN INDUCED BY 
PICOSECOND LASER PULSES 

W. Lee Smith, J, H. Bechtel and N* Bloem.be r gen 

ABSTRACT 

This report contains the consolidated results of investigations, 
carried out in the period January 1974 - June 1976 on the damage 
threshold of transparent optical materials. Single picosecond pulses at 
1. 06 fxm, 0. 53 (im and 0. 35 pm were obtained from a mode-locked 
Nd-YAG oscillator -amplifier -frequency multiplier system. The 
pulses are Gaussian in space and time and permit the determination 
of breakdown thresholds with a reproducibility of fifteen percent. It is 
shown that the breakdown thresholds reported here are characteristic 
of the bullc material, which include nine alkali halides, five different 
laser host materials, KDP, quartz, sapphire and calcium fluoride. 

The extension of the damage data to the ultraviolet is significant, 
because some indication is obtained that two- and three -photon absorption 
processes begin to play a role in determining the threshold. Throughout 
the visibile region of the spectrum the threshold is still an in<:reasing 
function of frequency, indicating that avalanche ionization is the 
dominant factor in determining the breakdown threshold. This is 
confirmed by a detailed study of the damage morphology with a high 
resolution microscope just above the threshold. A detailed discussion 
of the influence of self -focusing is given, and evidence for beam distortion 

iii 



below the power threshold for complete self-focusing is presented, 
confirming the theory of Harburger, 

This report constitutes essentially the material of the Ph. D. 
thesis of Walter Lee Smith, submitted in partial fulfillment of the 
requirements for the degree of Doctor of Philosophy in the subject of 
Applied Physics, Harvard University, September 1976. 

Chapter I describes briefly previous work on laser induced break- 
down and the relationship of the present work to that described in 
technical report 643, Laser Induc«?d Damage in Solids , May, 1973, 
by D. W. Fradin. Chapter II gives a description of the picosecond 
laser system. Chapters HI-V give the experimental results of the 
breakdown at three different wave lengths, in the near infrared, the 
green and the near ultraviolet. This material is presented mostly in 
the form of reprints, or of preprints of papers submitted for 
publication. 

The four appendices include reprints of papers on two -photon 
absorption in semiconductors and multiphoton photo-electric emission 
of metals, which have also been made possible by the existence of the 
reproducible picosecond pulse generator. 

All this material is assembled in the present report because of 
its relevance and potential usefulness for high power pulsed laser 
systems, including damage criteria for windows and other components 
in the ultraviolet. 
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CHAPTER I 


INTRODUCTION 


The investigation of the properties of matter as a 
function of the intensity of electromagnetic fields at 
optical frequencies has been stimulated by hhe advent of 
lasers, notably of pulsed lasers with high peak powers. 
Many nonlinear optical properties have been studied, but 
the ultimate nonlinearity of any condensed phase is a 
breakdown process, which transforms the material into a 
plasma. The damage threshold of optical components is an 
important limiting factor in the design of high-power 
laser (HEL) systems. The incipent processes of plasma 
formation are responsible for the establishment of damage 
thresholds in transparent optical materials. 
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The process of irreversible structural alteration of 

an initially transparent material substance, upon excitation 

from a sufficiently intense light wave, is known as laser- 

induced dielectric breakdown. 

Although the first pulsed laser operation was reported in 
1 

1960 by Maiman, no observations of breakdown processes were 

made. The early pulsed ruby lasers of 1960 operated in a "burst 

mode," with an optical output consisting of a collection of many 

relaxation oscillations of the system. It was not until the 

2 

technique of Q-switching, which allows temporal concentration 

of all the stored optical energy into a single output pulse of 

— 8 

about 10 sec duration, was implemented that breakdown was in- 
duced by a laser. At the Third International Quantum Electronics 

3 

Conference in 1963, Maker et al . , who used a Q-switched ruby 
laser of about 10 MW peak power, made the first communication 
of laser- induced breakdown. They reported the breakdown of air 
by focused pulses with an estimated breakdown field strength of 

Q 

10 V/cm. In addition, they reported (without elation) the first 
laser-induced breakdown in solids, relating that the "experimental 
factor causing the greatest difficulty was breakdown or burning 
of the sample" used for their third-harmonic generation experi- 
ments, 1965 brought a great improvement in laser peak power 
production, as Mocker and Collins^ demonstrated the first simul- 
taneous mode-locking and Q-switching of a ruby laser by use of 
a saturable absorber. Mode-locking techniques resulted in an 


i 
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immediate hundred-fold advancement in the peak power obtainable 
from a laser oscillator. Power levels exceeding the gigawatt 
mark were then obtainable with modest amplification. Again the 
blight of breakdown was noted in an original paper, as Mocker 
and Collins reported that "considerably higher mirror damage was 
found" when the laser was mode-locked. 

Throughout the remaining years of the 1960 *s, laser-induced 
breakdown (or laser-induced damage) thwarted laser experimenters, 
in particular those involved in the development of high-power 
laser systems. The reason for such longevity of the state of 
confusion concerning laser-induced damage was the complexity 
and variety of the involved processes. During those years, 
progress was slow and investigations were largely empirical in 
nature. But by 1970, the elements of laser- induced damage prob- 
lem were being unraveled, due in large part to improvements in 
laser system mode-control and reproducibility. Laser-induced 

5 

damage due to absorbing impurities had by that time largely 
been relegated to a quality control production problem. It was 
realized that self-focusing of powerful pulses leads to inten- 
sification of the light wave in a material, but it was not under- 
stood whether self-focusing caused the actual material disruption, 
as well. Indeed, many different processes were still being dis- 
cussed as possible breakdown mechanisms, among them photo-exci- 

g 

tation of populous impurity levels, stimulated Mandel'shtam- 

7 8 

Brillouin Scattering, bulk mu.Ttiphoton ionization, and electror 


REPRODUCIBILITY OF IHL 
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9 

avalanche ionization. 

In 1971/ several advances were made. The introduction of 

10 . 11 

short-pulse TEA CO 2 lasers enabled Yablonovitch to present 
conclusive evidence that avalanche ionization is the intrinsic 
physical mechanism which conveys disruptive quantities of in- 
frared light energy to the transparent crystal lattice in laser- 
induced breakdown. His experiment was conducted with a wave- 
length (10.6 ym) long enough that self-focusing and multiphoton 
ionization were absent. Impurity effects were minimized by 
tight external focusing and further eliminated by microscopic 
observation. With the identity of the intrinsic breakdown 
mechanism in hand, further progress in understanding the inter- 
play between self-focusing and breakdown was rapid. Guiliano 
and Marburger reported direct streak photography observation 

of pairs of moving self- focal spots in sapphire, verifying the- 

13 

oretical interpretations of time-dependent self-focusing theory. 

. 14 

In 1972 Yablonovrtch and Bloembergen explained the ubiquitous 
diameter of the filamentary self-focusing damage track, which is 
formed by the moving self-foci, with an argument based on ava- 
lanche electron production. Fradin in 1973 demonstrated 

that intrinsic laser-induced breakdown at 1,06 ym and at 
6943 gl6 

is also an avalanche ionization process, and they^"^ 
measured the laser pulse duration dependence of the breakdown 
threshold in NaCl. 

The experiments listed above (beginning with Yablonovitch) 


t 



placed laser-induced breakdown on a firm experimental footing 
for the spectral region {10.6 pm to 6943 S) and pulse duration 
region (~80 nsec to 5 nsec) investigated. However, other areas 
of research remained uninvestigated and many questions were still 
unanswered. Laser- induced breakdown in the picosecond pulse 
regime was not yet studied, except for Tradings NaCl datum, and 
it was not known whether avalanche ionization would in fact be 
effective for such short pulse durations. It was argued that 
submicro scopic inclusions would dominate the breakdown process 
for picosecond pulses. And, most intriguing of all, it was not 
known what would be the light frequency dependence of the break- 
down processes. Although avalanche ionization was established 
as the intrinsic low-optical-frequency breakdown mechanism, 

multiphoton ionization was known to be an additional intrinsic 

19 

process in matter that was expected to xnfluence the breakdown 
process as the light photon energy increased beyond 1.7 eV (the 
ruby laser value) toward the material optical band gap. These 
questions addressed the basic physical behavior of transparent 
matter under intense illumination, yet they also derived con- 
siderable practical importance from various high-power laser 
projects underway then and now. This report recounts 
three years of experimental effort to answer them. 

It is appropriate next to describe the avalanche ionization 
process that was introduced in the above chronology of laser- 
induced breakdown. Avalanche ionization breakdown may be 
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separated for discussion into three stages; 1) initiating 

electron production; 2) avalanche ionization plasma production; 

3) transfer of a damaging quantity of energy from the light 

wave to the material lattice by the dense, absorbing plasma. 

The first stage, at low photon energies {<1 eV for a typical 

insulator) , relies on the presence of shallow electron traps 

(color centers, inclusions, etc.) to yield starting electrons 

early in the laser pulse. As the photon energy is increased, 

such shallow traps may be ionized by a single photon, or direct 

interband multiphoton transitions may be invoked to supply the 

initiating free electrons. The second stage of the avalanche 

breakdown process concerns the evolution of an initiating 

electron in the strong optical light wave. Energy is supplied 

to the oscillating electron by the optical field. At the same 

time, the electron undergoes inelastic collisions with the 

-1 15 -1 

lattice at a rate i of roughly 10 sec . If the light 
field is sufficiently strong, the electron will succeed in 
attaining an energy that will allow it to produce a secon- 
dary electron by impact ionization of a lattice atom. The 
two electrons then may repeat the energization process to 
produce four electrons, those four may produce eight, and the 
avalanche has begun. The growth of the avalanche is described 
by an avalanche ionization coefficient n: 

N(t) = Hq exp 

op 

PMS !S POOE 


g{E(t') } at' 


( 1 ) 
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Here N(t) is the number of free electrons per unit volume, 

is the starting density of electrons, and E{t) is the 

laser pulse electric field distribution. The function r\ is 

14 

highly nonlinear in the electric field E(t). Fig. I.l 

illustrates this point. This graph displays E(t) , the quantity 

log[N(t)/N^] computed using experimentally derived data for 

n in Nad, and also the computed lattice temperature T(t) for 

a typical pulse from our experiments in NaCl. One observes 

that the electron distr , -lution is formed predominantly at the 

peak of the laser pulse, indicating strong nonlinearity and 

the threshold character of avalanche ionization. The quantity 

10 -3 

Nq typically has the value 10 cm (discussed in Chapter III, 

Part C, in detail) and it is often observed that breakdown 

18 ““3 

occurs if N(t) exceeds ~10 cm 

The production of the plasma at the very peak of the pulse 
serves to separate, for purposes of discussion, stage 2 from 
stage 3. In stage 3, the plasma absorbs energy from the final 
portion of the light pulse. The rate of energy absorption by 
the plasma at a time t, W(t) , may be approximated with the 
Joule heating expression of 


W(t) = 


e^ T„ |E(t) r N(t) 


m^ 


1 + W 


( 2 ) 


Here e and m* are the electron charge and effective mass, 
is the electron momentum relaxation time, and w is the light 
angular frequency. As described in Part C of Chapter III, 



t, in psec 


Fig. I.l Time dependence of the plasma 
growth [N(t)/NQ] and lattice temperature 
[T(t)] during breakdown induced by a pico“ 
second laser pulse with a Gaussian electric 
field distribution [E(t)]. 
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with a simplifying neglect of diffusive or radiative thermal 
loss during the picosecond heating, the temperature of the 
lattice may be calculated by performing an integration of Eq. 

(2) over time and using the heat capacity of the material. The 
T{t) curve of Fig. I.l is the result of such a numerical cal- 
culation. Although the calculation is not accurate above the 
sample melting temperature (800 ®C) , it nevertheless illustrates 
the manner in which lattice heating occurs subsequent to plasma 
production at the laser pulse peak. If the lattice heating for 
a given pulse succeeds in melting the lattice over a material 
volume sufficiently large to be observed under a high-power 
microscope (as related in Chapter III) , then production of 
breakdown would be recorded for that pulse. 

The qualitative picture of the avalanche ionisation process 
given above will serve to acquaint the reader with the basic 
characteristics. The notions that were introduced will be 
discussed in detail in following chapters. 

The most recent review article on laser- induced breakdown 

19 

is that by Bloembergen in 1974. That review and references ' 

. 20 

therein, the proceedings of the Conferences on Laser Damage 

21 

of 1969-1975, and a paper by Holway and Pradin will introduce 
the reader to the published literature on the subject. References 
to several other relevant publications which have appeared since 
1974 are included in the bibliographies at the end of each of 
the chapters in this dissertation. For information on DC 
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22 

breakdown,, the reader is referred to the book by O' Dwyer, 
and to a recent experimental paper and its references. 

Chapter II of this report describes the mode-locked 
laser oscillator-amplifier system which was constructed in 
our laboratory beginning in late 1973. Construction details 
which contributed to the stability of the system are reported, 

A full discussion is given there of the employed techniques 
for characterization of individual laser pulses. These tech- 
niques enabled us to measure optical frequency electric field 
distributions in absolute units with an accuracy of ± 15 
Several new observations were made during the calibration of 
our laser system. They were subsequently published and repro- 
ductions of those publications are included here as appendices . 
Appendix I describes the effects of laser pulse depletion, by 
linear and nonlinear absorption, in the common two-photon 
fluorescence method of measuring the duration of picosecond 
laser pulses. If not properly treated, depletion effects 
lead to erroneous pulse duration measurements. Appendix II 
describes a nonlinear optical method - the tA technique - for 
the extraction of an individual picosecond laser pulse duration 
with only the use of two detectors, an oscilloscope, and a 
harmonic generation crystal. This technique offers important 
advantages over other picosecond pulse duration schemes in- 
cluding streak camera photography, two-photon fluorescence, etc. 
Using the tA technique, we measured the degree of fluctuation 



in the pulse duration of our laser which is mode-locked 

by a saturable absorption dye. That experiment is recounted 

in Appendix III. Such fluctuation is intrinsic to the sat« 
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urable absorber mode- locking process^ but it had not been 
considered as a source of error in experiments prior to ours- 
As demonstration of the utility of the laser system and the 
characterization techniques described in Chapter II, papers 
concerning two topics other than breakdown are included as 
Appendix IV. Part A of Appendix IV describes the observation 
of a nonlinear photei=>lectric emission effect from a tungsten 
metal surface which involved two-photon absorption by the 
photoelectrons. Careful measurements of the two-photon 
absorption coefficient in several semiconductors are related 
in Part B of Appendix IV. The experiments of Appendix IV 
derived important benefits from the use of low-energy but 
high-peak-power laser pulses. 

The description of the central experiments of this 
research begins with Chapter III. That chapter describes 
our research at the 1.064 urn wave length, and it is composed 
of three parts. The first (Part A) is a paper recounting a 
measurement of the breakdown threshold, with our nominally 
30 picosecond pulses, in fourteen transparent solids. These 
solids include several practically important optical materials 
(KH 2 P 0 ^, a laser glass, a laser crystal, Si 02 and CaF 2 ) in 
addition to nine alkali-halides. The data for the alkali- 
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halides is compared with previous data taken, with nanosecond 
laser pulses by Yablonovitch^^ and Fradin ^ to 

show that avalanche ionization of electrons is indeed the break- 
down agent at 1.064 lim, even for picosecond pulses. Simultaneous 
measurements of the nonlinear refractive index are obtained and 
compared with other published values, to demonstrate that our 
experiment accounted properly for self -focusing. Part B of 
Chapter III is a follow-'up experiment in four additional state- 
of-the-art laser component materials. Part C of Chapter III 
describes one final aspect of the 1.064 ym laser-induced break- 
down study, its morphology. It was discovered that the pico- 
second breakdown morphology consists of a collection of micron- 
sized, spatially distinct sites, localized along the focal axis. 
Each of the microsites is the vestige of a tiny plasma v^hich 
grew from a small number, perhaps only one, of electrons which 
were present in shallow traps when the laser pulse traversed 
the focal region. Because of the short laser pulse duration, 
the plasma sites did not grow to the point of overlap. The 
observed density of sites allowed us to measure for the first 
time the actual density of free or slightly bound electrons 
which may serve to launch an avalanche inside a transparent 
dielectric. Photographs of the damage sites are presented in 
this part of Chapter III. 

Chapter IV recounts our effort to extend the work of 
Chapter III to a shorter wavelength than had been studied 
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previously. By frequency doubling the 1.064 ym pulses in a 
noncritically phase-matched crystal, we obtained 5321 ^ 
pulses with a spatio-temporal waveform as noiseless and re- 
producible as at the fundamental wavelength. Self-focusing 
constraints in this experiment were severe. Nevertheless, 
breakdown thresholds were obtained in six materials. Compar- 
ison of the thresholds at the two wavelengths showed a most 
exciting result: for the small band gap material KH 2 PO^, 

the threshold at 5321 S was less than the infrared threshold, 
indicating for the first time a quantitative observation of 
multiphoton (three-photon, in this case) absorption assistance 
to the avalanche. For the other materials, the 5321 S thres- 
hold equaled or exceeded the infrared threshold, in a systematic 
fashion, providing conclusive data also for the first time for 
a frequency-dependent increase in the threshold of avalanche 
ionization. This increase was predicted for some time, as 
was the effect of multiphoton absorption on the threshold. 
Morphology of this breakdown showed consistently that the 
operative process at 5321 S is avalanche ionization for the 
large band gap materials. Even for KH 2 PO^, the morphology 
confirmed that an avalanche dominates the breakdown process 
in stages two and three, as discussed earlier. The exper- 
imental results concerning the frequency dependence of the 
two intrinsic breakdown mechanisms constitute the highlights 
of Chapter IV. 
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Chapter V describes our investigation of breakdown at 
3547 S in the ultraviolet wavelength range. For the necessary 
frequency summing process to produce this wavelength, no material 
is known which allows noncritical phase-matching to be achieved. 
Severe difficulties were therefore caused by pulse spatial dis- 
tortion due to walk-off, and by self-focusing, in the nonlinear 
generation crystal. These distortions, time- and power -dependent , 
coupled with the smallness of the critical power for self-focusing 
in the samples, did not permit precise breakdown threshold inter- 
pretation to be made from the data at this wavelength, as discussed 
in Chapter V. Nevertheless, the data obtained at 3547 S verified 
qualitatively the frequency dependent behavior first discussed 
in Chapter IV. 
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CHAPTER II. 


LASER SYSTEM 

A. Introduction 

' This chapter contains a description of the diffraction- 

limited/ mode-locked Nd:YAG laser system which was constructed 
in our laboratory beginning in late 1973. The heart of this 
1.064 ym laser system is a passively Q-switched and mode-locked 
Nd;YAG oscillator, which produces a train of roughly 20 optical 
pulses, each of 30 picosecond duration. A single pulse is 
separated from the pulse train by a fast Pockels cell electro- 
optic shutter. The use of a single picosecond pulse optical 
source, rather than an entire pulse train, is greatly preferred 
for straightforward data analysis and attainment of conclusive 
results. The final elements of the laser system are two single- 
pass, NdiYAG amplifiers. These amplifiers provide an energy 

8 

boost of up to roughly 100, enabling a peak power of 3x10 
watts to be achieved in a single pulse. Peak powers of this 
Q magnitude make possible the efficient conversion of this 1.064 

ym light into the visible and ultraviolet harmonic wavelengths, 
c. We have observed power conversion efficiency of ~505 to 5321 2, 

'^20% to 3547 2, and ~15% to 2661 2, limited by the damage prop- 
erties of the nonlinear crystals. Because both the spatial and 

o 

temporal profile of these optical pulses is Gaussian (with some 
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reservation at 3547 S ) , calculations for extraction of results 
from data, whether analytic or numerical, are often simplified. 

The Gaussian spatial profile is convenient for focusing the 
light pulses to achieve higher intensity since the spatial profile 
remains Gaussian while only its width changes. The focused 
pulses from this laser system make available intensity levels 
higher than that sustainable by any solid, liquid, or gaseous 
form of matter: all are converted to the plasma state at such 
excitation levels. However, the utility of this laser system 
stems not from the available raw power, but from the reproduc- 
ibility and careful characterization of the waveform. These 
two qualities will be discussed further in the last section 
of this chapter and in later chapters as well. 

The following sections of this chapter contain details 
of the laser system components, operation, and performance. 

B. Mode-locked Laser Oscillator 

1. Resonator Design 

The Brewster configuration resonator is diagrammed 
schematically in Fig. II. 1. Between the two reflectors are 
positioned the laser rod, mode-locking cell, and spatial mode 
selectiiig aperture. The front and rear reflectors are di- 
electric-coated glass substrates with reflectivities of 46 % 
and 100%, respectively. The reflecting surface of the rear 
mirror is concave with a 3m radius of curvature. This 
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Fig- II. 1 Schema-tic illus-tra-tion of the 
oscillator configuration of the YAlG:Nd 
mode-locked laser system. 


curvature relaxes appreciably the alignment tolerances found 
in a flat-flat resonator configuration and reduces the need 
for frequent realignment. 

2. Oscillator Head 

The oscillator head design is a simple single 63° 
ellipse configuration. A cross-sectional drawing of the 
head is shown in Fig. II. 2. The single flash lamp is posi- 
tioned \t one focal axis of the ellipse and the laser rod 
at the other. The elliptical cross-section has axes of 2,750 
and 2.450 inches, and the cavity is 2.0 inches long. The head 
pieces were machined of brass and the inner cavity surfaces 
were then hand polished with successive grit sizes down to 
jeweler's rouge, until a mirror smoothness was obtained. The 
cavity pieces were then electroplated with electroless Ni 
and then with Au. Gold was chosen over Ag because of its 
greater durability and slightly higher reflectivity for light 
in the absorption bands of the Nd active ions. In addition, 
the slightly lower reflectivity of the Au compared to Ag in 
the ultraviolet reduces solarization of the laser rod. 

The laser rod itself is a Yttrium Aluminum Garnet rod 
doped with Nd, purchased from Allied Chemical Co. The laser 
wavelength is 1.064 ym. The rod ends are cut at Brewster's 
angle with respect to the rod cylindrical axis, for three 
reasons. 1. The Brewster configuration alleviates the fre- 
quency-selective properties of pairs of parallel, normal 
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Fig. II. 2 Illustration of 
the single'-ellipse configu- 
ration of the laser heads. 
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surfaces inside the resonator. Frequency selectivity would 
be deleterious for the production of bandwidth-limited output 
pulses. 2 . It removes the need for low-reflectivity dielectric 
coatings / with their attendant damage susceptibility, on the 
rod faces. 3. It results in 100% linearly polarized output 
radiation. 

The laser rod is surrounded by a tightly fitted Pyrex 
glass tube, or cladding, which serves to increase the pumping 
efficiency. The Pyrex glass also absorbs ultraviolet light 
and reduces rod solarization, thereby lengthening rod lifetime. 

An explanation is in order for our choice of Nd:YAG as a 

laser material rather than less expensive Nd-doped glass. As 

mentioned in a previous paragraph, our emphasis in this work 

has been on pulse reproducibility and characterization, in order 

to maximize control in our experiments. Ndrglass mode-locked 

lasers are now well known to be less suited for that class of 

+3 

applxcation, for the following reason. Because the Nd ion 
fluorescence linewidth is broadened substantially more by host 
interaction with glass than with crystalline YAG, the Nd: glass 
oscillator is capable of producing pulses correspondingly 
shorter in duration than a Nd;YAG oscillator. The Nd:glass 
pulses arc typically, in initial states, shorter even than 
the response time of the dye used to mode-lock the laser. This 
response time is reportedly about 7 psec. ^ Therefore the Nd; 


glass pulses are not smoothed and transformed, by the many 



successive passages through the bleachable dyer into a 
standard final waveform to the degree of Nd;YAG pulses. The 
result is that the temporal waveform of Nd: glass mode^locked 
laser pulses generally is not smooth but has noise-like 
structure of picosecond scale which fluctuates from shot to 
shot. Further problems of Nd: glass involve temporal broaden- 
ing due to self-phase modulation which progressively increases 
the width of successive pulses in the same pulse train. These 
factors render Nd: glass markedly unsuited for our purposes. A 
further inconvenience of Nd: glass is that its lesser thermal 
conductivity limits the firing rate to about 1 pulse per second 
whereas our Nd:YAG system can be operated at 30 pulses per 
second. This fast rate is beneficial for alignment procedures, 
as well as for data recording. 

Our Nd:YAG rod is cooled by distilled H 2 O which flows 
through a 1 mm thick annular space between the cladding and 
an outer Pyrex glass tube. The flow velocity is approximately 
10 cm/sec and the H^O temperature is maintained at 20 °C ± 

0.1 by a Lauda -Brinkman K2/R thermo-regulating circulator. 

No other elements are cooled with this circuit. This minimal 
thermal load allows the thermo-regulator to operate very sen- 
sitively to control the rod temperature for increased stability 
The Pyrex glass tubes were bonded to the gold-plated brass end 
pieces with a transparent epoxy (Eccobond 24) . 

The oscillator flash lamp, and those for our amplifier 
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heads, were manufactured by ILC Technology (#L"1419) and are 
xenon-filled. Xenon lamps are superior to krypton- filled lamps, 

an alternate for pumping Nd, for use in high-current, pulsed 

23 6 

operation. They also offer longer lifetime, approaching 10 

shots. The lamp dimensions are 4 mm inner bore diameter and 

2.125 inch cathode-anode separation. 

The oscillator flashlamp and both the rods and lamps of 
the dual amplifiers are water cooled by a circuit entirely 
separate from that for the oscillator rod, as mentioned 
above. A high-capacity industrial water cooler (Blue M #pcc- 
3455C) maintains a 20 “C temperature. This cooling system 
has the capability to service several additional amplifier 
stages or adjacent lasers. 

3 . Mode-locking Dye Cell 

A schematic illustration of the construction of the dye 
cell is shown in Fig. II. 1. The cell is of so-called "contacted" 
design, meaning that the flowing sheet of dye is in physical 
contact with the rear, 100% reflective mirror. This arrange- 
ment is optimum for the suppression of satellite structure in 
the mode-locked output pulse train. The sheet of dye solution 
is pumped upward against gravity and flows in a thin (0.5 mm 
thickness) sheet between the two confining surfaces of the rear 
reflector and the other uncoated substrate. The dye solution 
itself is Kodak #9860 saturable absorber diluted with laser- 
grade 1, 2 - dichloroethane. The dye has a relaxation time of 
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7 picoseconds and a saturation (bleaching) intensity of 60 
MW/citi^. The proper concentration (22 ml of dye with 100 ml 
dichloroethane) results in a low intensity transmission factor 
(through the cell) of about 0.7 at 1.064 ym. The other elements 
of the dye circuit, which are not included in Pig. 11.1, are 
the dye pump, dye reservoir, filter, and flow lines. Black, 
opaque polyvinyl flow tubing was used to shield the dye from 
ambient UV room light which would hasten its degradation. The 
Millipore teflon filter, installed in the circuit just before 
the dye cell, removes particulate matter from the pump gears, 
etc. Without the filter, the 0.5 mm channel in the dye cell 
is the most restrictive part of the circuit, and the light- 
absorbing, particulate matter collects and leads quickly to the 
damage of the dielectric coating. The stainless-steel reservoir 
is never completely filled, allowing the air volume to act as a 
cushion to damp fluctuations in the flow velocity. 

4. Flashlamp Pulsed Power Supply 

The two power supplies including the pulse-forming net- 
works used to fire the flashlamps in our laser system were 
manufactured by ILC Technology. This supplier also makes our 
flashlamps, and in fact, the power supplies are designed from 
the beginning for the exact firing characteristics of the L-1419 
xenon flashlamps. Two power supplies are used here, one for the 
oscillator lamp and one for the series pair of amplifier lamps. 
The maximum repetition rate is 30 pulses per second, limited 
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by the onset of "latch-on” or faulty continuous lamp ener- 
gization. This rather high repetition rate is made possible 
by special recharge delay circuitry designed to inhibit "latch- 
on". The circuitry begins recharge only after a sufficient time 
has elapsed for the ionized gas in the lamp to recombine to the 
neutral state. The oscillator (amplifiers) power supply charac- 
teristics are a maximum energy per pulse of 20 (50) joules, a 
pulse width of 150 (200) ysec, a maximum continuous power of 
1500 (2000) watts, and a charging voltage of 0 to -2hV. 

C. Pockels Cell Single Pulse Selector 

As mentioned in the introduction to this chapter, in 
order to use mode-locked lasers for many applications, it is 
desirable to employ only one of the pulses from the charac- 
teristic train of pulses. If the entire train is used in an 
experiment, the output signal is usually an average or accu- 
mulation over the many pulses. Unambiguous and precise con- 
clusions are at least more difficult and often impossible to 
draw from such data. Therefore we employ an electrooptic 
shutter after our oscillator to extract from the pulse train 
a single picosecond pulse. 

The general technique for single pulse selection makes 
use of the different transmission factor of an optical polar- 
izing element for the two orthogonal polarization states of 
the laser light. An electrooptic device is positioned between 
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a pair of crossed polarizers. Unless a proper voltage is 
applied to the electrooptic device, no light passes through 
the second polarizer. When given the proper voltage pulse, 
the electrooptic element rotates the light polarization by 
1/2, allowing it to proceed through the second polarizer 
unattenuated . 

The electrooptic device may employ either the Pockels 
electrooptic effect or the Kerr electrooptic effect. Because 
of the toxicity of the nitrobenzene fluid normally used in 
a Kerr cell and the excessively high required voltage, we 
have employed a Pockels cell electrooptic shutter in our 
apparatus . 

The Pockels effect is a linear electrooptic effect that 
occurs in all piezoelectric crystals. The ideal transmission 
T of a Pockels cell is given by 


T = T^ sin^(im^rg^V/X) 


( 1 ) 


Here n^ is the ordinary refractive index of the electrooptic 
material (in our case, KD 2 PO^) , r^^ is the material linear 
electrooptic coefficient, X is the light vacuum wavelength, 
and V is the applied voltage. From this expression, one 
observes that the condition for full transmissio? (that for 
X/2 polarization rotation) is obtained with an applied voltage 
of 


^X/2 ” ^/^”o^63 


( 2 ) 
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known as the "half-wave” voltage. This quantity for KD2P0^ 
(KD*P) at 1.06 iim is approximately 6000 volts. 

The Pockels cell shutter used in our experiments is 
illustrated schematically in Fig. II. 3. A high voltage power 
supply (HV) charges a transmission cable (Cl) through a 

g 

charging resistor (Rl) of 1.1 x 10 ohms. This voltage V 
( 2 i 12,0 00 volts) is applied to one side of the electrode pair 
in the spark gap (SG) , the other side of which is at ground 
potential. High-pressure nitrogen gas was used to increase 
the dielectric breakdown voltage limit of the gap. When the 
focused initial pulses of the pulse train enter the shutter, 
their polarization is rotated by \/l by the half-wave plate 
(WP) , allowing them to pass through to polarizer P2, There 
the pulses are deflected and focused into the electrode gap 
region where they initiate breakdown. A voltage pulse of V/2 
then travels through cable C2 for a duration 2fi,^/v, where 
is the length of cable Cl and v is the propagation velocity 
of the cable. (Cl, C2, C3 are all 50 ohm, RG-8/U cable.) 

While this voltage pulse is traveling across the Pockels cell 
PC, itself a 50 ohm impedance, traveling-wave device, the 
polarization of any light passing through it is rotated by 
A/2. The length of the cable Cl is made so that only one 
pulse can be switched out from the pulse train, as expressed 
by the condition 2£,^v~^ ^ 2L/c. Here 2L is the optical round- 
trip path length of the laser oscillator, and c is the vacuiim 



Pig. II. 3 Schematic diagram of the Pockels cell single pulse 
selector used in the laser system. Notation; WP, half-wave 
plate; P, and P^/ Gian prism polarizer; PC, travel ing-v:ave 
configuration KD 2 PO 4 Pockels cell; C1-C3, 50-ohm coaxial cables; 
Rl, 107 ohm charging resistor; R2, 50 ohm load resistor; SG, high 
pressure, N 2 ~fiHed, laser-triggered spark gap; PD, photodiode 
detector; HV, high voltage supply. 
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speed of light. The right side is the temporal interpulse 
spacing in the pulse train (7 nsec here) . This one pulse 
with polarization now rotated back to its original state is 
then transmitted through polarizer P2 and constitutes the 
selected single pulse. Subsequent pulses on the train are 
still focused into the spark gap, but no voltage is present 
then, due to the long recharge time of cable Cl as guaranteed 
by resistor Rl. The cable C3 conveys the voltage pulse to 
the 50 ohm grounding resistor R2, and C3 is long enough to 
disallow voltage reflections at the C3-R2 connection from 
reaching PC while light is still passing through it. The fast 
(<0.5 nsec rise time) biplanar photodiode (ITT F4000) PD is 
used to monitor the entire pulse train, and conveniently dis- 
plays which pulse in the train was selected and excised by the 
shutter. A typical oscillograph of the laser j.ulse train minus 
the selected pulse is shown in Fig. II. 4a, along with the actual 
switched-out single pulse in Fig. II. 4b. It was found that the 
pressure in the spark gap (nominally 160 Ib./sq. in.) could 
be varied sufficiently to allow the selected pulse to be taken 
from any position in the pulse train. In normal operation the 
pulse in the train center was selected. However, the N 2 pressure 
could be used in some cases to make adjustments of the single 
pulse intensity level. (This technique would not be possible 
with a Ndrglass laser because pulses in different positions in 
the train have widely different duration and substructure.) 



Oscillograph of mode- locked 
pulse train after switch-out 
(2 nsec/div.). 



Oscillograph of single selected 
pulse (0.5 nsec/div.}. 


Fig. II. 4 

Temporal Waveforms 
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The jitter, or fluctuation, in the selection position of our 
shutter was typically 15-30 nsec, or 2 to 4 pulse positions. 
The actual Pockels cell used in our system is a product of 
Lasermetrics Corp. (Model 1071) , and has a cylindrical ring 
electrode (CRE) geometry, as indicated in Fig. 11,3. This 
configuration provides a more uniform transmission across 
the electrooptic crystal face at a given voltage than does 
the more conventional end-plate electrode geometry. In 
addition, the CRE geometry reduces the cell capacitance and 
yields a faster cell response time (in our case, ~1 nsec) , 
important for clean switch-out of a single pulse. 

D. Dual Amplifiers 

The final stages of our laser system are two Nd:YAG, 
single-pass, traveling wave amplifiers. These amplifiers 
provide a maximum power gain of about 100 for the single pulse 
which is selected by the Pockels cell. The amplifiers have 
operated at this excitation level for a period of almost three 
years without suffering any damage. 

The amplifier heads are different from the oscillator 
head in only one respect. The Nd:YAG rods (from the same 
supplier and of the same size as the oscillator rod) have 
faces cut, not at Brewster's angle, but at a small wedge angle 
of 6°. Because the amplifiers are single-pass devices, the 
loss from the face reflectivity is not important compared to 
the alignment ease of the wedge cut over the Brewster cut. 
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The faces are cut at 6° rather than 0° for two reasons. The 
small angle eliminates problems caused by back propagation of 
intense reflections into the Pockels cell/ etc. And the 6° 
wedge angle was calculated to be just large enough to suppress 
low-order parasitic inodes which would otherwise be able to 
oscillate inside the rod between the faces, thereby decreasing 
useable gain. 

The cooling system and power supply systems were described 
in previous paragraphs. Although it was not necessary for any 
experiments undertaken thusfar, a different arrangement of 
these amplifiers could be introduced to provide an additional 
gain factor of perhaps 8 to 10. Because the first amplifier 
is not at all saturated by the single pass of the single pulse, 
it could be used instead in dual-pass mode. The second am- 
plifier would remain in single-pass formation. This overall 
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triple-pass amplifier chain would then provide over 10 watts 
in a single 30 psec output pulse. 


E. Laser Pulse Characterization 

For a laser pulse which is both spatially and temporally 
Gaussian, the intensity waveform I{r,t) in a plane normal to 
the propagation direction is expressible as 




( 3 ) 


Here r is the radial coordinate of distance from the 
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propagation axis and t is the time variable. Thus the 
intensity waveform is specified entirely by the three param- 
eters p, and T, which are the peak intensity, the spatial 
radius, and temporal radius, respectively. This condition is 
true for focused, as well as for unfocused, Gaussian laser 
pulses. Therefore the task at hand in order to specify quan- 
titatively I(r,t) is to record I^, p, and t values for each 
laser pulse. The techniques for these three measurements 
which we have developed and employed in our experiments are 
described in the following paragraphs. Because the physical 
processes studied in our experiments (as discussed in later 
chapters) are highly nonlinear, even exponentially nonlinear, 
in tne intensity, precise knowledge of the pulse parameters 
is a solemn necessity. The unusual care that was exercised 
in the accurate determination of the individual pulse param- 
eters in our experiments was of fundamental importance for 
our observations . 

1. Measurement of Pulse Spatial Dimension, p 

The measurement of the unfocused pulse radial dimension 
offered the least difficulty of the three parameters. In 
initial experiments, measurement of the spatial profile was 
made with a multi-lens camera photographic technique^ which 
is independent of all film exposure nonlinearity. Photo- 
densitometric analysis of the data film yielded the parameter 
p. As one would expect, fluctuations in the spatial profile 
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p 

2 

and thus the area A (=11 p ) were small (unobservable at 1.06 ym) , 
as the amplifier stages were not used in the nonlinear ^ saturation 
range. Therefore, a measurement (or rather, the average of 
several measurements) of the pulse area needed to be performed 
only once for an experiment, rather than for every shot. This 
was fortunate, considering the time involved in the photographic 
recording and analysis ^ 

In later experiments, far greater accuracy in the measure- 
ment of p was made possible by the purchase of an electronic 
device, a (Reticon) self-scanning photodiode array. This 
monolithic integrated circuit linear array is composed of 256 

w 

individual Si photodiodes, each of area .001 in, x .001 in., 

I, spaced by .002 in. (50.8 ym) . Associated circuitry repetitively 

scans the 256 diodes (at 10 kHz) in order and processes the 
individual dicde light signals into a serial output format. 

This serial output signal is fed into an oscilloscope and 
provides a direct display of the distribution of light illu- 
minating the array over the scanning period. The array is in 
practice illuminated with a single pulse and an oscillograph 
^ is made. The pulse area is obtained immediately and to an 

accuracy of better than 105. The spectral response of the 
device is such that it is useful at 1.064 ym and at the second-, 
third-, and fourth-harmonic wavelengths. The photodiode array 
has been invaluable in another respect in experiments using 
the harmonic wavelengths, especially the third-harmonic. The 


,/ 
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array was used to phase-match the KH 2 P 0 ^ third-harmonic 
generation crystal to the particular angle which yielded the 
smoothest ultraviolet spatial waveform. 

2. Measurement of Pulse Temproal Duration, t 
a. Measurement of Average Duration, <t> 

A major series of experiments in our characterization 
of the laser pulse output concerned the ifteasurement of the 
average temporal duration, t. Specifically, t is the 1/e 
half duration of the intensity. Two techniques were used for 
these measurements: i) a second -harmonic autocorrelation 

5 

(SHAG) technique, and ii) a two-photon fluorescence (TPP) 

5 

technique. Fig. II. 5 illustrates the optical arrangements 
for the two techniques. In both techniques, the incident 
1.064 ym pulse is split into two equal parts which are then 
recombined to interact in a nonlinear medium- In the SHAG 
technique. Fig. 11.5a, an optical delay is imposed on one of 
the separated pulses by a movable mirror in one arm of the 
interferometer. The two pulses then travel collineaily through 
a nonlinear second -harmonic generation crystal (KH^PO^) . There 
the two infrared pulses produce 5321 S light with polarization 
orthogonal to the infrared light. In addition to a constant 
background signal due to the individual pulse harmonic gener- 
ation, there is a correlational contribution which is propor- 
tional to the temporal overlap of the two pulses. The amount 
of second-harmonic light produced ip measured by detector D. 
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0 . SECOND HARMONIC GENERATION AUTOCORRELATION TECHNIQUE 



b. TWO-PHOTON FLUORESCENCE AUTOCORRELATION TECHNIQUE 



c. AUTOCORRELATION DATA 


DELAY OR 
DISTANCE 


Fig. II. 5 Schematic diagram of two 
conventional experimental methods for 
measuring the duration of picosecond 
laser pulses ^ and an illustration of 
their autocorrelation output data. 
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If the detected harmonic light is plotted versus optical 
delay, a curve as shown in Fig. II, 5c is the result. From the 
width of the peak, one extracts t. In the TPA technique, no 
delay is necessary. The two pulses are directed ant ico 11 inear ly 
through a two-photon absorbing meditim, such as Rhodamine 6G, 
and the fluorescence spatial pattern is photographed by a 
camera. Photodensitometric analysis of the film again yields 
a curve as shown in Fig. II. 5c, and again t is extracted from 
the peak width. An important consideration of the contrast 
ratio between the peak height and the background height is also 
necessary in these experiments. Because both experiments require 
many tens or even a few hundred shots, each produces a measure 
of the average duration, <t>. The significance of that obser- 
vation will be mentioned below. The result of both experiments 
was that for our laser <t> = 18.0 ± 1.8 picoseconds and the 
full-width-at-half -intensity duration is 30.0 ± 3.0 picoseconds. 
This result has been confirmed by more recent, direct experiments 
using streak cameras. 

In Appendix I, we include a paper which was published as 
a result of our TPP investigations. In the work in that paper, 
account for the effects on the contrast ratio from one- and/or 
two-photon attenuation of the laser pulses in the TPP medium 
was included for the first time. These observations were 
important for the extraction of accurate <x> values from the 


data. 
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b. Measurement of Individual Pulse Duration, x 
The previous section described our measurements of 
the average duration <t>. For a laser that is mode*-locked 
with a saturable absorber, knowledge of <x> is insufficient 

■ a 

for many purposes because the single pulse duration, t, fluc- 
tuates substantially from shot to shot. Saturable absorber 
mode-loc'^ing is intrinsically a stochastic process. Letokhov 

7 

and others developed a successful model of the process by 
beginning with the realization that, prior to absorber bleach- 
ing, the radiation in the laser resonator is well characterized 
as a narrow-band Gaussian random process. This behavior results 
from the beating of the many modes which are brought above 
lasing threshold and which lack a deterministic phase relation. 

Q 

The well developed formalism of Gaussian random processes 
then immediately yields probabilistic expressions for the 
parameters of the noise radiation, such as the most probable 
noise spike duration, the spread therein, etc. It is the 
distribution of probability for the temporal width of the most 
intense noise spike, in the initial stages of pulse development, 
which results in a distribution of duration for the output 
pulses . 

Thus one must detect and account for this variation in 
duration from shot to shot if he is to arrive at accurate 
knowledge of the single pulse intensity distribution. Prior 
to our experiments, the only way in which this fluctuation 
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had been considered was by direct observation with streak 
photography. Only large laboratories have such equipment. 
Therefore we developed a nonlinear optical technique for the 
indirect monitoring of the pulse duration, which requires 
minimal equipment and has several advantages over direct streak 
camera photography for experiments such as ours. This technique, 
called the TA technique, is described in a published paper 
which is included in full as Appendix II. The tA technique 
involves the simple measurement of two energy signals with 
photodiodes for each pulse. One signal is the energy in 

a sample of the laser pulse, and the other is the energy 
at the second-harmonic wavelength which was obtained by sending 
a second sample through a KH2pO^ crystal. It is shown in 

r\ 

Appendix I that the ratio of is proportional to the 

individual pulse duration t and to the pulse area A. Fortunately, 
for our laser we could treat A as a constant in most cases. By 
recording the signals for each of the shots in an 

experiment and forming their average, one may renormalize that 
data into invididual pulse duration data, in units of pico- 
seconds. This tA technique was invaluable for obtaining 
accurate single pulse intensity data in all our subsequent 
experiments, and is mentioned further in later chapters. 

To make use of the tA technique itself to investigate 
the statistical nature of the pulse duration, an experiment 
was performed to measure the pulse duration probability function 
of our laser. This experiment is also described fully in a 



paper, presented as Appendix III. In that experiment the 
pulse duration probability function was found to be approx- 
imately 20 picoseconds wide, indicating pulses of about 20 
to 40 picoseconds are produced by our laser. This degree of 
fluctuation has been since verified directly with streak 

Q 

photography by workers at the Naval Research Laboratory - 
It is evident that such wide variation in one of the three 
pulse parameters (p, x, 1^) must be taken into account if 
accurate experiments are to be performed. 

3. Measurement of Individual Pulse Peak Intensity, I^ 

Present optical detectors do not have sufficient 
frequency response to monitor the instantaneous intensity 
of a picosecond laser pulse. They must be used in the 
integrating mode. Furthermore, the large surface area of our 
detectors performs an integration over the spatial coordinates 
of our laser pulse. Expressed more concisely, the signals we 
record are energy ( B ) signals according to the relation 
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where I(r,t) is the Gaussian intensity function for our pulses 
as written earlier in Eq. 3 of this chapter. Yet, intensity 
(I^; values, not energy values, are needed for our data evalu- 
ation. Therefore we are forced to measure p, x, and B for each 
pulse in order to arrive at I^^ for each pulse. The techniques 
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for measurement of p and t were discussed in the previous 
sections of this chapter. It remains only to outline the 
relatively simple procedure used to measure the energy of 
each laser pulse. In almost all of our experiments, biplanar 
photodiodes, which have characteristics of high linearity of 
response, stable sensitivity over long periods of time, and 
large photocathode area, were used as the optical detectors. 

They were generally fitted with narrow-band dielectric 
interference filters. These detectors were always calibrated 
with the same procedure, using the same two thermopiles as 
standards. The thermopiles are products of the Eppley Lab- 
oratory (Newport, R.I.), and were calibrated in absolute units 
by the manufacturers using National Bureau of Standards reference 
light sources . Our calibration procedure involved firing the 
laser at a non-damaging level at 7 pulses per second into the 
thermopile and recording the signals from the thermopile and 
the photodiode. The calibration figure of the thermopile was 
then used to determine the photodiode absolute sensitivity. 

The high repetition rate w;s necessary for two reasons: i) 

the thermopile is too insensitive to respond to single pulses; 
ii) the thermopile is calibrated with a continuous light source, 
meaning that our repetition period had to be less than the 
instrument response time by at least a factor of 10 to enable 
us to use the calibration figure. 

Therefore, the typical pulse intensity measurement involved 
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recording the pulse energy with a calibrated photodiode, and 
dividing it by it ^ P t, where p is the essentially constant 
spatial radius and x is the individual pulse duration obtained 
with the tA technique {Appendix II) . Our procedures yielded 
i-ntensity values with estimated accuracy in the unfocused case 
±20^ (±10^j in electric field), and +305 (±155 in electric 
field) in the focused case. These uncertainty figures were 
substantiated by the reproducibility of standard measurements 
performed at different times. 

F. Additional Applications of Laser System 

To indicate the utility and versatility of the laser 

system and characterisation procedures that were discussed in 

this chapter, in experiment;" other than those constituting the 

remaining chapters and principal theme of this thesis, papers 

concerning two different research areas are included here. 

Part A of Appendix IV contains a paper published on the subject 

of laser-induced photoemission from a metal surface caused by 

multiphoton absorption. In particular, four-photon photoemission 

was observed from a tungsten surface, using 1.064 ym light over 

8 2 10 2 

an intensity range of 5 x 10 W/cm to almost 10 W/cm . The 
photocurrent density varied over a corresponding range of 
slightly more than 4 orders of magnitude. The careful control 
of the picosecond laser pulses permitted the elimination of 
lower-order thermionic heating effects, a necessity for the 
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observation of such a high-order nonlinear process. 

In Part C of Appendix IV/ we include a paper published 
on the different subject of two-phot m absorption in semi- 
conductors. The two-photon absorption coefficient at 1.064 ym 
was measured in five semiconductors. The coefficients we found 
were among the smallest ever reported, and were in surprisingly 
good agreement with theoretical extimates. The reason for the 
systematic smallness of our coefficient values, as compared with 
those from previous experiments, is again the degree of charac- 
terization and control permitted by our experimental apparatus 
and techniques in determining the intensity distribution of our 
laser pulses. 


i 
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CHAPTER III. 


STUDY OF PICOSECOND LASER-INDUCED BREAKDOWN AT 1.064 ym 

This chapter is composed of three parts. Part A (pages 
III-2 to III-IO) describes the first systematic study of pico- 
second laser- induced breakdown, and the major results and 
conclusions of this chapter are presented there. Part A is 
a reproduction of a 1975 publication. In Part B (pages III-ll 
to II 1-12) the experimental techniques developed in Part A are 
used to study several new laser-related optical materials. 

Part B is also a reproduction of a published paper. Part C 
(pages III-13 to III-26’' describes an additional aspect of 
the 1.064 ym study, the breakdown morphology, which was not 
included in Parts A or B. Part C is in the form of a manu- 
script which is to be published in Optics Communications 
(1976) . 
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Dielectric-breakdown tlireshold and nonlineai^refractive-index measurements with picosecond 

laser pulses* 

W. L. Smith, J. H. Bechtel, and N. Bloembergen 
Gordon McKay Laboraiory, Han/ard UnivarsUy, Cambridge, Massachusetts 02138 
(Received 18 February 1975) 

Single pulses from a well-calibrated mode-locked YAGjNd laser have been used to study the dielectric- 
breakdown propel ties of 14 transparent solids. Results are presented for the breakdown threshold electric 
fields, accurate to un estimated 15%. The experimental technique permits simultaneous measurement of the 
nonlinear susceptibility (— <o), via the nonlinear refractive index. Susceptibility results with 
estimated accuracy of 50% are presented. 


I. INTRODUCTION 

It is well known that damage to transparent opti- 
cal components in high-power laser systems con- 
stitutes a limiting factor to their performance. In 
many practical cases the damage threshold is set 
by the presence of absorbing inclusions, or by the 
occurrence of self-focusing. It has, liowever, been 
established that an intrinsic damage threshold, 
characteristic of the built transparent material, can 
be defined. If the effects of inclusions have been 
eliminated, and if one has corrected for the effects 
of-self -focusing, the threshold is determined by the 
onset of avalanclie ionization, in a manner wliich is 
very similar to low-frequency electric breakdown,^ 
Since the conduction-electron density has to build 
up from the initial density of about 10® to 10^° elec- 
trons/cm® to a critical value of about 10‘® electrons/ 
cm® during the duration of the laser pulse, the 
damage threshold has a characteristic dependence 
on pulse duration. Most available experimental data 
were taken with pulses of nanoseconds or longer 
duration. In fact, heretofore, there was only one 
quantitative datum point taken with a picosecond 
pulse.® The reliability of that point has been ques- 
tioned on the grounds tliat, for picosecond pulses, 
impurities of submicroscopic size may play a domi- 
nant role and that no moi'phological study of the 
damaged region had been made. The purpose of 
this paper is to provide systematic picosecond data 
on the damage threshold in 14 transparent materi- 
als. The data obtained are characteristic of the 
material and are very reproducible. 

It is necessary to take into account the effects of 
seif -focusing in tins experiment, even though the 
maximum pulse power was always less than the 
critical power for self- focusing. This power 

is related to the intensity-dependent contribu- 
tion to the index of refraction where 

( 1 ) 

by 

12 


P„ = cAV32?r®H2 . (2) 

In the above equations, W(, is the linear refractive 
index, c is the speed of light, and A is the vacuum 
wavelength. Wagner et discuss tlie deforma- 
tion of a Gaussian light pulse fay the self-interaction 
represented by Eq. (1). Let u>i be tlie beam waist, 
i, e,, the radius at which the intensity drops to 1/e 
of its on-axis value, at low power levels. Then the 
low-power focal area irzij is reduced by a factor 
(1 -P/f’jf) due to self-action. The maximum inten- 
sity is consequently increased by a factor {1 -P/ 
Ppr)‘*. Zverev and Paslikov® and Fradin® have sug- 
gested how this effect may be used to determine P„, 
and consequently tin, by measuring tlie damage 
threshold as a function of focal length, utilizing dif- 
ferent lenses. 

Consider a low-power diffraction-limited Gauss- 
ian pulse with an initial l/t> intensity radius p. A 
lens, witliout spherical aberration and wilh focal 
lengOi/, will produce a focal spot ivith a radius (to 
the 1/e intensity point) of 

H'i=A//2jrp. (3) 

If the pulse total power is P and we define an area 
A = uicf, the focal plane intensity distribution is 
given by 

/{r)-(P/A)e-<‘^^'“i’\ (4) 

For liigh-power pulses where self -focusing tenden- 
cies cannot be ignored, the relationship between the 
maximum intensity and the power becomes 

^^^(P/AKl-P/P^)*! (5) 

or 

+p;J . (6) 

If one plots, therefore, the reciprocal measured 
power, necessary to produce the value of re- 
quired for breakdown, as a function of the recipro- 
cal focal area one should obtain a straight line. 
The slope of this straight line gives the brealtdown 
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rms electric field strength 

where is the free-space permittivity. The in- 
tercept with the vertical ajds determines JP„ and b^. 

One restriction implicit in the preceding discus- 
sion concerns the response time of the self-focus- 
ing mechanism. Various mechanisms for self- 
focusing exist, with response times from the or- 
der of 10"® sec for electrostrictlve and thermal 
self-focusing, to the order of 10'^® sec for self- 
focusing originating in the electronic hyperpolar- 
izability. For picosecond-pulse experiments in 
solids, the only mechanism fast enough to reach its 
steady-state value is the electronic hyperpolariz- 
ability. That mechanism follows our laser-pulse 
intensity essentially instantaneously, and so the 
assumption of a steady-state nonlinearity is a valid 
one for our experiment, hideed, picosecond pulses 
provide results free from ambiguity introduced by 
slower contributions to n^, an ambiguity which may 
hinder Q-switched laser-pulse experiments. 

The method described above has been used to 
measure the rms breakdown field strength Fj and 
the critical power for self-focusing Per several 
transparent solids of interest. Section II describes 
the experimental apparatus and procedure, includ- 
ing a technique for measuring the product of the 
laser-pulse duration and pulse area on a shot-to- 
shot basis in a simple manner. Results are pre- 
sented and discussed in Sec. HI, 

II. EXPERIMENTAL METHOD 
A. Laser system 

A mode-locked YAG: Nd laser system generated 
the picosecond pulses used in these experiments. 

The experimental arrai^ement is illustrated in Fig. 
1. A 6-mm-diam by 3-in-long Brewster-ended 
YAG: Nd rod supplies gain. Mode lo citing is 
achieved by a contacted dye cell which confines a 

0. 5-mm-thick sheet of flowing Kodak 9860 saturable 
absorber. TEMqo transverse profile is established 
with an intracavity aperture. Single-pulse switch- 
out is effected by a cylindrical-ring electrode potas- 
sium dideuterium phosphate Pockels cell, positioned 
between crossed polarizers, and pulsed by a nitro- 
gen-filled high pressure laser-triggered spark gap, 
A selected pulse then propagates through two single- 
pass YAG* Nd amplifiers for augmentation of ener- 
gy. A repetition rate of 20 pulses/ sec is available, 
and an amplified single-pulse pealt power exceeding 
300 MW has been measured. CRO-1 is a Tektronix 
519 oscilloscope used in conjunction with high-speed 
biplanar photodiode PD-1 (ITT F4000, S-1 re- 
sponse), for single-pulse menitoring. Two Gian 
polarizers form the variable attenuator A in Fig, 

1. The calibrated photodiode PD-3 (S-1 response) 
measures the energy incident on the lens. The sec- 



FIG. 1. Schematic diagram of the experimental ar- 
rangement. MLO, modo-loct:ed YAG : Nd oscillatory PC, 
Pockaia celly LTSG, laser-triggered spark gap; AMP, 
YAG :Nd amplifier; BS, beam splitter; F, filter; PD, 
biplanar photodiode; CRO, oscilloBOOpe; KDP, second- 
harmonic generation crystal; A, Double Glan-prism 
variable attenuator; L, focusing lens; S, dielectric sam- 
ple. 

Fig. III.l 

ond-harmonlc generation channel, ending with 
photodiode PD-2 (S-20 response), is used with PD-3 
to obtain single-pulse values of the product of the 
pulse duration and pulse area, as will be described 
later. CRO-2 is a Telctronix 555 dual-beam oscil- 
loscope. In Fig, 1, L represents the focusing lens 
and S the sample materials tested in the experi- 
ments. 


B, Pulse characteristics 
/. Spatid citaraciertstics 

Since precise measurement of threshold optical 
electric fields requires an accurate Icnowledge of 
both the spatial and temporal pulse shape and the 
pulse energy, great care was talten toward those 
measurements. 

At the actual site on the optical table where the 
focusing lens was later placed for the breakdown 
measurements, three scans across the beam were 
made with a 50 -pm pinhole to determine the trans- 
verse Intensity distribution. The results of one of 
the scans is shown in Fig, 2, In addition, pictures 
were taken of phase-matched second-harmonic light 
with a multilens camera® (MLC). The result of 
each of the above tests was that the intensity dis- 
tribution is of Gaussian form to at least the 1/c® 
points, and that the average l/e intensity radius 
(p) is 1,1 ±0,1 mm. 

Z Tetnporal clwractertsflcs 

In order to determine the temporal pulse prolile 
of the output, two tecluiiques^® were used: a sec- 
ond-harmonic -generation autocorrelation method 
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DISTANCE FROM PULSE CENTER (mm) 

FIG, 2, Intensity distribution of laser pulse at loca- 
tion of focusing lens. The solid Une Is a Gaussian fit to 
the data, yielding a 1/e intensity radius 1. 1 mm. 

Pig. III. 2 

(SHAG), and two-photon-fluorescence (TPF) photog- 
raphy. The two methods are formally equivalent, 
hoth providing 1 + 2G2 (t). In eitlier method, a single 
pulse Is split into two pulses of equal amplitude 
(any disparity in the two amplitudes degrades the 
contrast ratio). In the TPF arrangement, the tivo 
pulses are then redirected by mirrors to overlap 
in a precisely collinear fashion within a two-quan- 
tum-absorptivo material such as Rhodamine 6G, 
and the dye fluorescence is photographed. In the 
SHAC method the beam splitter is part of a Michel- 
son interferometer. The pulses are recombined 
collinearly, after one has been delayed by a vari- 
able amount, and directed through a potassium di- 
hydrogen phosphate (KDP) crystal. One plots sec- 
ond-harmonic output as a function of delay. Tlie 
result of one of the SHAC runs is shown in Fig. 3. 
Tile enlianced scatter of the points near the center 
of the graph is due to interference effects in the 
interferometer. The curve is a four-parameter 
least-squares fit to the function A'l +K^ e?£p[- 2(s 
- It shows a contrast ratio [equal to 

(Ai +Ks)/Kt] of 2,7:1. The results of TPF traces, 
each one the result of many laser shots, agree with 
the data of Fig. 3: the average full-ividth-at-half- 
maximum (FWHM) temporal duration of the laser 
output is 30.0 ±6.0 psec, implying an average 1/e 
half-duration (t) of 18. 0±3. 6 psec, 

3. Indh'idimUpiilse t, A meosnrements 

During the course of the laser calibration, a nar- 
row distribution of individual pulse durations was 
actually observed, the average of which is the val- 
ue stated above. In order to investigate the pulse- 
vddth statistics further, and to be able to obtain in- 
dividual T measurements, a simple pulse-duration- 
pulse-area monitor system was devised and added 
to the laser system. “ The rationale of this monitor 
system follows from the proven Gaussian spatial 
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and temporal distribution of the pulse intensity. Let 
5jr be the energy content of the sample of the funda- 
mental pulse that is measured by photodiode PD-3 
in Fig. 1, and (Sgi, be the energy content of the sec- 
ond-harmonic sample that is measured by PD-2. 
Then, for small harmonic conversion efficiencies, 
it can be show that 

( 8 ) 

where is a collection of constants for a given ex- 
periment. If, in the course of an experiment, 
enougli data points are taken so that the average of 
the ratio values given above is meaningful and rep- 
resentative, then that average value can be set 
equal to the product of the known, previously mea- 
sured average values, <t) (A>, thereby calibrating 
absolutely the individual pulse values of the i*atio 
given by Eq. (8). Using these individual-pulse val- 
ues of tA, one obtains individual-pulse intensity 
values from the easily measured pulse-energy val- 
ues. This simple scheme is useful in obtaining ac- 
curate data for processes wluch are highly nonlin- 
ear in the intensity, such as laser-induced break- 
down or multiphoton photoemission. 

C, FontI coiiskicralioiis 

As discussed in Sec, I, a plot of the reciprocal 
brealtdown power versus reciprocal local area 
yields the breakdown field and the nonlinear re- 
fractive index. For each material tested in tliese 
experiments, the breakdown power was measured 
for three different lenses. The focal lengths were 
0.5, 1.0, and 1.5 in. All were of proper shape to 
minimize spherical aberivition. In the absence of 
all aberration, the intensity distribution in the focal 
plane inside the test material would be the Fourier 
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FIG. 3. nesiiltB of measurement of the laser piitse 
averaae temporal duration, using a seeontl-harmonic- 
goneradon teehnlnue. The enhanced scatter of t!ie polnta 
in the middle of the data is due to interference effects in 
the Intcrferometor. The solid curve is a least- squares 
fit to the data, yielding a FWHM duration of RO, o psec 
and a contrast ratio o£ 2. 7;i. 
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transform of the distribution in the bade focal plane 
of the lens. The latter is of the form I(r) =/oexp 
- (r/pf, and so the focal-plane distribution would 
be f^(r)--^f5exp- where is i;iven in Eq. 

(3). Such performance would be diffraction limited. 
In reality, however, contributions to the focal radi- 
us due to spherical aberration by the lens and by 
the planar surface of tlie test material must be con- 
sidered. The radius that would result from the 
first of these two aberrations iCj is given by^® 

jcj =4«p®//® , (9) 

where k is a function of lens refractive index and 
has the value 0.069 for the 0.6- and 1. 5-ln. focal- 
length lenses (quartz, »(j = 1,4493), and 0,0455 for 
the 1 -in. focal-length lens (glass, «q = 1.6667), at 
1.06 pm. Equation (9) was derived for the case of 
a plane wave incident on the lens, but it is still use- 
ful here in a manner to be described below. The 
focal radius that would result frorn spherical aber- 
ration by the planar test material surface” ii)^ is 
given by 

iCa = D(>Jo - l)pV2)io/^ , (10) 

where B is the distance between the interface and 
the focal point and, in the experiments here, was 
always 2.0 mm. The effective focal radius «■ is 
related to the above constituents by 

m'=^(m| + h| + h|)'^® . (11) 

In Table I, the quantities introduced above are tabu- 
lated for the three lenses used here. Because of 
the shallow depth of focus D maintained here, the 
contribution to m’ from «’j ts seen to be less than 1% 
in all cases. The contribution from is also seen 
to be less than 1% for the /=l-in. and / = 1.5~in. 
lenses. Thus, in these two cases, the use of Eq. 

(9) to compute h’s is acceptable. For the / = 0. 5-in. 
lens, however, the spherical aberration contribu- 
tion to fc as calculated from Eq. (9) is nonnegligible. 
It was felt necessary to perform an experiment to 
determine the focal distribution for this lens. To 
that end a chromium-steel razor blade was mounted 
on a two-dimensional precision translator with 
0. 13-pm resolution. The front surface of the blade 
edge was then positioned precisely in the focal plane 
of the 0.5-in. focal-length lens, the focal plane 


Table III.l 

TABLE I, Focal parameters. 


Focal length / 
(in.) 

re, 

<fim) 

H'2 

(Mm) 

«'3 

(Mm) 

(Mm) 

A -■ 

(Mtn=) 

0.5 

1.95 

a 

0.23 

3.3» 

34,2 

1.0 

3.80 

0, 37 

0,03 

3. 94 

48.2 

l.s 

5,85 

0.25 

0.01 

5.85 

107,9 


709 

being identified as that requiring the least irradi- 
ance to produce a spark on the blade su’*face, T!;c 
blade was then swept in increments across the 
beam, and the fraction of tlie incident pulse energy 
(as measured by a reference detector) that was not 
blocked by the blade was measured as a function of 
position. Several shots were takeit at each posi- 
tion, and the irradiance was always maintained be- 
low tlie blade damage level. A simple calculation 
shows that the transmission function T(x) should 
he of the form 

T(x) = (irzc®)-^ r f dy 

= I [1 + erf (jf/ic)] . (1 2) 

The experimental data conformed well to a curve 
of the above form if ic was given the value 3. 3 pm. 
Microscopic examination of the blade edge after the 
experiment verified that it had not been disfigured. 
Therefore, the value of 3. 3 pm for the focal radius 
for the/=0. 5-in. lens was used in all data analysis 
in this report, and it appears in Table I. 

Finally, aberration due to decentration and tilt” 
of the lens must bo considered. In these experi- 
ments, such further contributions to tr were ren- 
dered negligible by tlie following alignment pro- 
cedure, which was used for every experiment. 
Firstly, in order to deal with tilt, the lenses were 
mounted in holders which allowed the attachment of 
a microscope slide In front of the lens in such a 
manner as to be parallel to the plane of the lens. 
Next the beam was apertured so as to allow only a 
small (0. 5-mni-diara) central section to strike the 
attached slide. The lens was then adjusted to be 
normal to the propagation direction to within better 
than 1,0“ by adjusting the holder so that the reflec- 
tion was collinear with the pulse propagation direc- 
tion. Secondly, decentration was made negligible 
by two operations. Before the lens was placed in 
the system, the positi'>n of the center of tlie laser 
pulse at a point past the lens site was located by a 
marked phosphorescent card. Then with the lens 
inserted, and corrected for tilt as discussed above, 
it was translated so that the center of the trans- 
mitted light was coincident with the position pre- 
viously located. Additionally, the reflection of the 
pulse from the curved, first face of the lens was 
observed with another viewing card and when align- 
ment was obtained, the reflection was symmetrical- 
ly distributed about the propagation direction. The 
lenses could be centered to within 0. 1 mm easily, 
making associated aberration negligible. The re- 
producibility of lens positioning afforded by the 
above simple techniques is necessary for accurate 
absolute breakdown field measurements over a 
large series of materials, 
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D. Protedure 

To begin each experiment, the sample was opti- 
cally polished and blown with dry nitrogen to re- 
move any residual abrasive particles from the sur- 
face. The sample was then clamped onto an x-y 
translator, and the entrance surface was aligned 
normal to the propagation direction. It was neces- 
sary to maintain a hot-air stream over the most 
hygroscopic of the materials to retain surface qual- 
ity. The amplifier gain and the variable attenuator 
were then adjusted so that the approximate brealt- 
down irradiance entered the material and so that 
sufficient signal was available for the rA measaiir- 
ment devices, as discussed earlier. Thereafter, 
the attenuator was left unchanged in ordt that the 
ratio values of Eq. f8) retain their usefi .. ,ss, as 
noted earlier. An> irther fine adjustment of the 
input irradiance was done with the amplifier gain 
control. The material was translated between 
pulses so tliat a new portion of the material was 
probed by each pulse, and the coordinates of each 
probed location were recorded. A fast photodiode/ 
Tektronix 519 oscilloscope combination {PD-1/ 
CRO-1 in Fig, 1) was monitored during each firing 
in order to eliminate data from pulses witii satel- 
lites. At each firing, the fundamental energy and 
the harmonic energy were recorded. Loss of en- 
ergy by reflection from the lens and material sur- 
face was taken into account in each experiment. 

The fundamental energy was measured with a photo- 
diode calibrated with an Eppley Laboratory tliermo- 
pile. Breakdown sparks were initially detected with 
well dark-adapted eyes operating in a dark room 
shielded from laser flashlamp liglit. Latei*, using 
the recorded coordinates of tl.a probed sites, the 
crystal was observed under a microscope to verify 
breakdown at any visually uncertain points. The 


n 

observed damage was confined to a eylindrically 
shaped volume typically 2 //m in diameter and 25- 
50pm in length, and therefore occupied a small 
fraction of the focal volume, just above threshold. 
Tlie reproducibility of the data when probing dif- 
ferent specimens and different volume elements 
within the same specimen indicate tiiat a property 
characteristic of the bulk material is observed. 

Only a few isolated cases of inclusions wore ob- 
served under the microscope, and they were elimi- 
nated, A further detailed account of tlie morphology 
of the damage will be given elsewhere. Typically, 

50 shots were used near the threshold region with 
each lens for each material. The data were then 
plotted as shown in Fig. 4 for KCl. The siiarpness 
of the threshold is typical of the behavior encoun- 
tered in these experiments. The uncertainty in the 
breakdown energy was generally less than 10%. 

The three encircled points were associated with 
pulses having unusually short tA values, and did 
cause breakdown even though their energy was be- 
low the tln*eshold energy evident in the figure for 
pulses with tA nearer to The damage ener- 

gy was determined in this manner in each case, and 
then (t> was used to convert the energy to power. 

A pulse energy Gf is related to the pulse power P 
by P = The reciprocal damage-power 

values were then plotted versus reciprocal focal 
area as discussed in Sec. I. An example, that of 
KCl, is shown in Fig. 5. A least-squares fit to the 
three points then yields the desired quantities— t!ie 
breakdown field E^, independent of self-focusing ef- 
fects, and the nonlinear index of refraction n^, via 
the critical power for self-focusing, 
m. RrsuLTs 

In this section results are presented from tests 
on 14 materials, AH samples were obtained from 
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Fig. III. 4 

FI(5. 4. KcBUlts of 
breakdown experiment 
In KCl for /'l. Ti-ln. lens. 
The .average pulse dura- 
tion was usetl to con- 
vert energy to power via 
P=S/k}^'Hr), as dis- 
cussed in See. Ill), The 
three circled points re- 
sulted from atypically 
short pulses. 
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commercial sources, and all except SiOg and ED-4, 
"vvere raonocrystolliiio. All results are for a [j;©- 
ometry ^vith the propagation direction along the 
[001] axis and electric vector along the [100] axis, 
except for the YAG : Nd material, in which the 
propagation direction was along a [111] axis. This 
sample orientation was chosen in order to measure 
the same quantities that are operational in YAG: Nd 
laser rods, 

A. Pulsi^duration dependence of the breakdown tlircsliold 
in NaCl 

Graphs such as Fig. 5 yield the breakdown- 
threshold electric field according to Eq. (6), as 
discussed previously. Field values so obtained 
are independent of all self-fo casing effects. The 
field values obtained by us are all absolute, and 
the experimental uncertainty, estimated by repro- 
ducibility and by calculation, is ±15%. 

In tills section we examine the pulse-duration de- 
pendence of the breakdown threshold in alltall 
halides probed by pulses of 1.06- jam radiation only. 
Measurements have been made in NaCl for four dif- 
ferent pulse durations, from 10 nsec to 15 psec, by 
previous investigators,^ Those data are reproduced 
in Fig. 6, along with the datum point for 30 psec 
measured by us. Tlie solid curve in that figure is 
the result of an approximate theoretical calculation 
discussed in Ref. 2. The previous investigators 
described the pulse duration dependence of the 


Fig, III. 5 



FIG, 5. Uesults of breakdown esporimcnta In KCl for 
three lenses. The solid line ia a loast-sciiiarcs fit to the 
data. Tile elope equals the reciproeal threshold break- 
down intensitv and the vertical-axis intercept equals tiie 
reciprocal critical power for self-focuslnK. 



rms OPTICAL ELECTRIC FIELD (MV/cm) 

FIG. 6. Functional relationship between the break- 
down threshold electric field strength and the pulse dura- 
tion. The solid curve is a semlemplrical prediction for 
NaCl, diBOUSsed In Ref, 2, 

Fig. Ill, 6 

threshold in terms of a field-denendent ionization 
rate inversely proportional to the pulse width. 
An exponential relation between the number of con- 
duction electrons and time was obtained 

N(f)=Noexp(jT* a[E(f')]dt') “NoA/(f) , (13) 

where M(f) describes the evolution of the initial 
number N(, (~ 10“ to 10‘® per cra“). As mentioned 
earlier, because Af(/) must reach -10‘° during the 
presence of the laser pulse, stronger fields are 
necessary as the pulse duration is decreased. 
Therefore the field dependence of the ionization 
rate, as well as the pulse duration dependence of 
the threshold, is displayed in Fig. 6. 

In addition to data for NaCl, tlie data now avail- 
able at 1.06 pm for five other alltali halides are 
illustrated in Fig. 6. The three Na halides display 
a similar relative dependence of threshold on pulse 
duration, which appears to be somewhat different 
for the K halides. In combination vdth dc measure- 
ments by others, the data of Fig. 8 make available 
the approximate field dependence of the ionisation 
rate in six materials. 
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FIG. 7. Experimental trend of the breakdown field 
through (a) the Nn halides and (b) the K halides. Shown 
for comparison arc data from previous studies; (x), 

Itef, 19; (A), Ref, 18; (•), dc data from Ref, 16, 10,6- 
fim data from Ref. 17. In the figure nsec denotes a 
pulse duration the order of nanoseconds, and psec, 
picoseconds. Each set of data Is normalized to the field 
for NaCl of that same data set. Error bars for the pres- 
ent work are shown displaced to the right. 

Fig. III. 7 


D. Vnrlntiun of the breakdown-threshold field through 
the Periodic Table 

We now consider the variation of the breakdown 
threshold among different substances for picosec- 
ond pulses at 1. 06- wavelength. The results of 
our measurements are presented in Table II. The 
threshold is observed to vary from a minimum of 
3.4 MV/cm to a maximum of 22.28 MV/cm among 


the materials tested. KDP possesses a threshold 
so high as to approach the level necessary to malce 
multiphoton ionization probable. Figure 7 displays 
the trend of our threshold measurements across the 
Periodic Table for the Na and K halides. For com- 
parison, previous data from dc, and (?-switched 
CO 2 , ruby, and YAG : Nd, laser-damage studies 
are presented. Note that all the data for a particu- 
lar combination of pulse duration and frequency are 
normalized to the NaCl threshold of that same data 
set. 

C. Tliird'ordcr electronic susceptibility 

The critical power for self-focusing, as given by 
Eqs. (2) and (6), is obtained in these experiments 
as the inverse of the intercept of the fitted line wltij 
the vertical axis, in graphs such as Fig. 5, The 
nonlinear index Hj and x*®* (- - w) are re- 

lated by 

Wg = 12irx‘®’ (- CO, w, w, - (u)/ho . (14) 

The experimental condition that tlie laser-pulse 
duration is much shorter than the electrostrictive 
response time sec) and that the materials 

investigated are ionic solids, indicate that the val- 
ues of P„, }t2, and X*®’ obtained in these experi- 
ments are tliose of electronic (E) origin only. Sev- 
eral techniques have previously been employed to 
measure among them third-harmonic genera- 

tion, three-wave mixing, intensity-depen- 
dent ellipse rotation, nonlinear birefringence, 
and time-resolved interferometry,’® 

The values obtained in the present experiments 
for the critical power the nonlinear refractive 
index, and the third-order electronic susceptibility 
XrJif (- u), w, - w) are presented in Sec. A of 

Table in. Also tabulated is the linear refractive 
index used to convert iin into x*”^> The experimen- 
tal accuracy of P„, and of x”*^» is estimated to be 
±50%, Section B of the table lists previously ob- 
tained values of x‘®*^ for comparison. The smallest 
value was registered by NaF, a factor of al- 
most 18 less than that of KBr. Agreement with 
previous values is substantial, except for LiF, 
given the rather large experimental uncertainties 


Table III. 2 

table II. rms broalcsown electric field strengths (Fj, uncertainty, ± ISTi;; 
field ratio uncertainty, ±10%). 



NaF 

NaCl 

NaBr 

KF 

KCl 

KBr 

KI 

Fj (MV/cm) 

10.77 

7,34 

3.67 

8.34 

5.86 

5.33 

5.87 


1.47 

1.00 

0.77 

t.l4 

0. 80 

0. 73 

0,80 




ED-4* 


Fused 




LIF 

Rbl 

glass 

YaG :Nd 

SiO, 

CaF; 

KDP 

Ft <MV/em) 

12.24 

3.40 

9. 90 

9.82 

11,08 

14.44 

22,28 

Ft/FjJ^':' 

1.67 

0. 46 

1.33 

1.34 

1.39 

1.97 

3,04 


“Owens-minole (Toledo, Ohio) undopotl loser glass. 
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TABLE III, ROBultQ for xfiVf w, - w). 


Material 

■^er 

<MW) 

A. 

«o 

Present worlc 

(10'*^ esu) 

110*'^ eau) i-50'J 

I 

±3007d 

B. (10- 

Experiment* 

11 m 

±20% 

ecu), other work 

Theory'^ 

Other I U 

NaBr 

0,111 

1.62 

96,15 

4,13 







NaCl 

0.16G 

1.532 

64.66 

2.63 

1.2 

1.7 



0.8 

2.0 

NaF 

1.13 

1.321 

D. 45 

0.33 

0.25 




0, 13 

0.3 

KI 

0. 035 

1,638 

112.35 

4.88 





3.1 

9.0 

KBr 

0.075 

1,544 

142. 31 

5.83 

2.8 

3.0 



1.3 

3.7 

KCl 

0.322 

1,479 

33.15 

1.30 

1.2 

1.9 



0.9 

1.0 

LiF 

0,433 

1.387 

24,37 

0, 90 

0.3 

0.2 

0.34 ±0,06 


0. 28 

0.3 

CaFj 

0.383 

1.429 

27.70 

1,05 


0. 35 

0.43 ±0.13 


0.4 

0.4 

fused SIO; 

0.733 

1.450 

13.63 

0,52 


0.7 

0,70 ±0,12 

0.39 ±0.03*''' 

1.4 

1.1 

KDP 

0,295 

1.494 

36.18 

1.43 







ED-4 

0.503 

1.550 

21.22 

0, 87 



1.07±0, 08 

0.71 ±0.05’'^ 



YAG : Kd' 

0,302 

1.823 

35.34 

1.71 



2,18±0.13 

(1,98 ±0,14*^ 
ll.53±0. 30^^ 


0. 62 


*1. Values from third-harmonic fieneratlon, VVan^f and Bardeen^ Ref, 20, II. Values from threo-ivavo mixing. Maker 
and Terluine, Rof. 22. III. Values from throe-wave mixing, Levenson and Bloembergen, Ref. 22. IV. Values from in- 
tensity-dependent ellipse rotation, Owyoung, Ref. 24. V, Values from tlme-rosolved Interferometry, Bliss etal., Ref, 
26. 

”l. Calculated from generalised Mlller^s, = esu, Wang, Ref. 27, II. Calculated from Wang's rule. 

Ref. 27. 

'All exporimmital values in this row, with one exception, are for pulse propagation along a |lll| direction, the usual 
laser rod configuration. In these cases, the measured susceptibility is +Xjjjf +2x}i|f). The value in the column 

designated IH Is for the single element Xitlf. however. The present work Investigated a YAG;Nd crystals sH other values 
in the row pertain to undoped VAG crystal. 
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listed. The only serious discrepancy occurs for 
LiF. The difference between our value for LiF and 
those obtained by others is a factor of approximate- 
ly 3. We tested two different samples of LiF, a 
total of three times, and our results were repro- 
ducible. Therefore, we know of no reason to sus- 
pect that the error in our value for LiF exceeds the 
stated uncertainty of ± 50%. The expected trend of 
the third-order susceptibility element through the 
periodic table is displayed for the Na and K halides. 

IV. SUMMARY 

The method of controlled dielectric breakdown®' '' 
has beeii demonstrated to be a useful way to mea- 
sure simultaneously the two intrinsic parameters 
which must be known to enable propagation of max- 
imum picosecond-light intensity through a material 
in a stable nondestructive manner. Those param- 
eters, the breakdown field strength and the elec- 
tronic nonlinear refractive index, have been mea- 
sured in several materials of scientific interest 
^vith 1.06-pm laser pulses of 30-psec average dura- 
tion (FWHM). The measurements were made with 


well-calibrated pulses and focal geometry, in a 
manner allowing elimination of any inclusion effect 
and measurement of intrinsic processes charac- 
teristic of tlie bull;. A simple technique was re- 
lated for the measurement of (individual picosec- 
ond pulse duration) -(pulse area) products. The 
values of presented here are in 

reasonable agreement with those of otlier investi- 
gators. The measured breakdown field strengtiis 
have been shown to be those of an avalanche ioniza- 
tion process. Further picosecond experiments at 
the second-harmonic frequency are in pr ''gr*'ss to 
investigate the onset and properties of dielectric 
brealidown driven by multiphoton ionization, 
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Laser^induoed breakdown and nonlinear refractive index 
measurements in phosphate glasses, lanthanum beryllate, 
and AI2O3 ' 

W. Lee Smith and J. H. Bechtel 

Gordoft McKay Laboratory, Harvard UnivErsity. Cambridge, Massachusetts 02138 
(Received 2 February 1976) 

We present results of measurements of the laser-induced dielectric breakdown threshold in two phosphate 
composition laser plasses, in lanthanum boryllateiNd, and in AljOj. The nonlinear refractive indices for 
these materials are also presented. The measurements were made with single 30-pscc 1.064-ftm pulses. The 
results are compared with those for other laser materials in the 1.06-ftm spectral region. 


FACS numbers: 79.20, Ds, 77.50. +p. 78.20.Dj, 42.65. ~k 

The properties of two recently introduced laser 
materials which operate near LOG pm have attracted 
considerable interest. The prospect of a smaller non- 
linear refractive Index but larger gain coefficient than 
offered fay familiar silicate glasses has created Inter- 
est in new phosphate composition glasses. ‘ Another new 
laser material, neodymium -doped lanthanum beryllate® 
(LasBcsOs :Nd) is of interest due to its biaxial mono- 
clinic (Cjh) structure which results In plane-polarized 
output radiation. We report here measurements of the 
threshold electric field for laser -induced damage and 
the nonlinear refractive index in these two new mate- 
rials, and also in a carefully grown sample of AlaOj. 
These new results will be compared to measurements 
from a similar study at 1.064 pm,® which included 
YAG :Nd and ED-4 laser glass. 

The experimental technique and apparatus used in 
this experiment have been described previously,® and 
will only be outlined here. A saturable absorber (Kodalc 
9860) mode-locked YAG rNd oscillator, constrained to 
the TEMq 5 mode, produced the picosecond, 1.064 -pm 
laser pulses. A fast Pockels cell/spark gap combina- 
tion excised one single pulse from the oscillator pulse 
train, and It was amplified as needed by two single - 
pass YAG ;Nd amplifiers. Subsequent detectors moni- 
tored the energy content and the temporal duration,® 
and identified pulses with satellites or bad switch-out 
structure. The pulses of this laser were determined 
previously to be of 30-psec average duration. Three 
precision lenses of the proper shape to minimize spher- 
ical aberration, with focal lengths of 1.0, 1.5, and 
3.0 in,, were used with each material. The calculated 
focal area A {defined as the area in the focal plane 
within which the intensity exceeds times the peak 
axial value) for the three lenses 'n this experiment was 
31.0, 50.8, and 180.5 pm®. These values include the 
small effects of spherical aberration by the lens and the 
planar sample surface as discussed in Ref. 3, 

The experimental procedure involved the passing of 
focused increasingly powerful pulses each through a 
new recorded position in the sample and noting the 
threshold power level P at and beyond which breakdown 
occurred. Breakdown was detected in by observa- 
tion of microscopic spark radiation with well dajk- 
adapted shielded eyes. Subsequently, these sightings 
were checked by 340-power microscopic examination, 
at \vhich time any inclusion effects would have been 
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spotted. A threshold power for breakdown was thus 
determined for each of the three lenses. The desired 
quantities, the threshold (rms) electric field for break- 
down, Ej, and the critical power for self-focusing. 

Par, were extracted using the operative equation® 

where is the peali intensity 
reached at the focal point. This equation describes the 
mild reduction of the pulse area due to self-focusing 
and assumes that the pulse Gaussian spatial profile it- 
self does not become distorted. It may be rewritten into 
a more useful form as 

( 1 ) 

This equation is valid if only minor self-focusing oc- 
curs, that is, for pulse Input power less than a critical 
power Per which is defined as 

P„ = c?l®/32/«2. (2) 

Here Hg is the nonlinear refractive index, \ is the opti- 
cal wavelength, and c is the vacuum speed of light. 
Noting the linear form of Eq. (1), the reciprocals of 
the three threshold input power levels are plotted ver- 
sus their respective reciprocal focal area. A linear 
least -squares fit to the three data points then yields 
the experimental values of and P^p, Eb Hg are 
immediately obtainable using Eq. (2) and the relation 
dcnai/Ho^oe), where hq is the linear refractive index 
and is the free-space permittivity. The real part of 
the electronic third -order nonlinear optical susceptibil- 
ity, “w), which is associated with jjg in 

this experiment may be obtained v.a the relation 

Hg= (12ir/«o)xJ®{i(-w,a),w, -w), (3) 

where the frequency degeneracy factor of 3 has been 
included in the factor 12. 

T^l3l.G XXX *4 


TABLE I. Brealidown (rms) olcctric flvlt! (UrDngths un- 
certainty, i LYi ; field ratio uncertainty, r 10^’), 


Material 

(MV/cm) 

PEI, :N(1 

ir,.r, 

UIG-.I :Nd 

irs.n 

LHG-fIzNd 

i:!.n 


12. 3 

ED-1 ^ 

!l,f) 

YAC 

3. ^ 


^DaCt from Ref. ;t. 
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table II. Mcafiurcd nonlinear refractive Itidleeo (f.s uncer- 
tainty, ± GO^'). 


Material 

do*"* esu) 

AljOj 

1,30 

LHG-5 :Nd 

1.30 

LUG-6 :Nd 

1,59 

BELtNd 

2,12 

ED-4‘" 

2.12 

YAG :Nd“ 

3.53 


*Data from Ref. .1. 


Table III. 5 

The results for the brealtdov/n rms electric field 
thresholds are presented in Table I. These values are 
independent of self-focusing effects, in that Eq. (1) 
corrects analytically for the present mild self-focusing, 
as just discussed. The breakdown threshold was sharp 
(^10% uncertainty in input power) for all the tested 
materials; inclusion effects were not observed. In this 
experiment the materials KCl and silica, examined in 
Ref. 3, were retested to check the accuracy of the new 
focal parameters. For both materials, the threshold 
values and critical power values agreed to better than 
10% with those obtained earlier. Therefore, one may 
compare the absolute values obtained in the present 
experiment with these for the various materials tested 
in Ref. 3. The tested phosphate coniposition glasses 
were types LHG-S and LHG-6 of Hoya Optics, Inc.^ 

The AlgOa sample was a carefully grown sample from 
Union Carbide Corp.'’ All the materials listed in Table 
1 except A1j03 are candidate materials for high-power 
near-l.OG-^im laser oscillators and/or amplifiers. 

The laser transition in LaaBejOs ;Nd (BEL :Nd) is at 
1. 070 or 1.079 pm, depending on orientation,® while 
in LHG-5 and LHG-6, the laser center wavelength is 
at 1.056 pm. The breakdown threshold values for the 
new materials show considerable Improvement over 
those of the two more commonplace 1,06 -pm materials, 
YAG :Nd and ED-4 glass. The Ej value of 15.6 MV/ejn 
registered by BEL ;Nd is 59% greater than that of 
YAG : Nd, which represents an Improvement factor of 
about 2. 5 in intensity. The phosphate composition 
glasses also show a marked improvement in damage 
threshold. 

Next we present the results for the nonlinear refrac- 
tive indices. Table II presents the values obtained in 
the present experiment and also lists the values for 
ED-4 laser glass and YAG :Nd measured by the identical 
technique in Hef. 3. The experimental uncertainty in 
these values is unfortunately inherently rather large, 
i 50%, due to the employed slope -Intercept method. The 
orientation of our AlaOj sample was such that the laser 
pulse propagated along the rhombohedrai axis of the 


crystal. The electric field direction within the normal 
plane was not determined. The orientation of the 
BEL :Nd sample was not determined; however, we ex- 
pect the orientational dispersion in not to exceed 20%. 
We now compare these measured values with those 
available from other methods. Recently, a value of 
1.48±0. 15xl0'‘® esu for Ha in ruby was obtained by a 
time -resolved interferometry technique. “ This com- 
pares well with our value of 1.30x10*^® esu. The two 
phosphate glasses gave Identical results, 1.50x10'*® 
esu, in our measurements, whereas measurements at 
Lawrence Livermore Laboratory’ have produced values 
of 1.16 ±0.12 and 1.01 ±0.10x10'*® esu for LHG-5 and 
LHG-6, respectively. The small nonlinear indices of 
these neiv materials have important implications with 
regard to the propagation and focusing behavior of 
terawatt laser pulses.* The m value of the other new 
sample, BEL :Nd, along with its relatively large dam- 
age threshold (Table I) and output polarisation charac- 
teristics, indicate interesting possibilities for new laser 
oscillators. For completeness we note that, while our 
results gave the same value of for BEL ; Nd and 

ED-4, other techniques offering better accuracy have 
found values for ED-4, or its Nd-doped version ED-2, 
of 1.52±0.3,“ 1.73±0.1,“and 1. 53±0. 15,“ in units of 
10'*® esu. For final comparison, other published values 
for undoped YAG in units of 10'*® esu are 3. 16± 0.6,“ 
4.09±0.3,“and 3.47±0.35.“ 
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port and critical review of the work. 


Tlesearcb supported In part hy the National Aeronautics anti 
Space Administration, by the Joint Services Electronics 
Program, and by the Advanced Hesearch Projects Agency, 

'Laser Focus 11 (No. Tj), 07 (1975); 11 (No. (5), 44 097.5); 11 
(No. 10), 40 (1975). 

®n.C. Morris, C. F. CUne, R.F. Begley, M. Dutoll, P.J, 
Harget, H.P. Jenssen, T.S. LaFrance, andB. Webb, Appl, 
Phys. Let-t. 27, 444 0975); sample supplied by Allied Chemi- 
cal Corp. , Morristown, N.J. 079C0. 

*W. Lee Smith, J.H. Bechtel, and N. Bloemhergon, Phys. 
Rev. D 12, 700 (1975). 

''Hoya Optics USA, Inc, Menlo Park, CA 94025. 

*Unioti Carbide Corp,, East Chicago, IN 40312. 

‘‘M.J. Moran, C.-Y, She, ondR.L. Carman, IEEE J. Quan- 
tum Electron. QE-11, 259 (1975), 

Milam and M.J. Weber (unpublished). 

"B.S, Bliss, D.R. Speck, andW.W. Simmons, Appl. Phys. 
Lett. 25, 72fl (1974). 

®A. Owyoung, IEEE J. Quantum Electron, QE-9, 1064 (1973). 


607 


Appl. Phys. Lett., Vol. 28, No. 10, 15 May 1976 


W.l_ Smith and J.H. Bechtel 


607 



ni-13 


PICOSECOND LASER -INDUCED DAMAGE MORPHOLOGY: 
SPATIALLY RESOLVED MICROSCOPIC PLASMA SITES* 


W. Lee Smith, J. H. Bechtel and N. Bloembergen 


Gordon McKay Laboratory 
Harvard University 

Cambridge, Massachusetts 0E138, U.S.A. 


The morphology just above threshold of damage sites caused by 
picosecond 1. 06 |xm laser pulses is shown to consist of a collection of 
micron-sized, spatially distinct vestiges of individual plasmas. From 
the observed site density, the density of electrons which may initiate 
breakdown can be inferred. The morphology is consistent with an 
avalanche ionization model, but not with absorbing inclusion damage. 


i 


ni-14 


Nt'Jv observations are presented here concerning the initiation and 
growth of regions of avalanche plasma inside initially transparen* solid 
dielectrics due to the passage of intense picosecond laser pulses. These 
observations provide experimental data for the :Q,rst time on the question 
of the operative density of electrons which may serve to initiate laser- 
induced breakdown at a given intensity level. Because of the short 
duration of the laser excitation used in this experiment, avalanche growth 
does not develop to the point of plasma microsite overlap and extensive 
disruption of the material. The damage morphology in the picosecond 
regime, induced by well controlled pulses just above the breakdown 
threshold, is shown to consist of a chain of individual, spatially -resolved, 
micron-sized plasma vestiges. They are limited to a narrow axial part 
of the focal region, the perimeter of which is well correlated with an 
avalanche ionization model. 

Previously we have reported thresholds for picosecond laser-induced 

1 2 

breakdown in eighteen transparent materials, ’ The present observations 

concern the morphology of the damage produced in those studies. Single, 

1. 06 [Jim pulses from a passively mode-locked, TEMqq, YAlG:Nd laser 

system were used to cause breakdown. The average full-width-at-half- 

maximum (FWHM) pulse duration was 30 psec and a nonlinear optic 

technique was used to infer the duration of each individual pulse. The 

4 

pulse duration is known to fluctuate in passively mode -locked lasers. 

This technique allows accurate single pulse intensity measurements. The 
influence of incipient seif-focusing was minimized by tightly focusing the 
laser pulses with an external lens and by always using incident pulse power 
less than the critical power for self-focusing. Corrections for incipient 
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self -focusing were taken into account by a technique using several lenses 

of different focal area. Any damage track, -which showed presence of 

inclusions upon subsequent microscopic examination, was eliminated. 

This occurred in roughly one percent of the shots. 

In figs, la, b, c we show photomicrographic data for picosecond 

laser -induced breakdown near threshold in a NaCl single crystal sample 

from Harshaw Chemical Co. These damage tracks were produced with a 

well-corrected 0.5 in. focal length lens, resulting in an intensity 1/e 

radius of 3. 3 pm. ^ Only one pulse traversed any single region in the 

cryst**.l. The tracks were photographed with a Leitz Ortholux microscope 

and blue light through a 0, 4 N. A. objective at 340x total magnification. 

The calculated resolution and the observed minimum detectable diameter 

is ~ 0. 7 pm. The scale lines in the figure indicate a distance of 10 

microns, and were transferred from a photograph made under identical 

magnification conditions of a 10 micron-spacing objective micrometer. 

Because the photographs were made under bright-field illumination, 

the damage microsites appear black in the photographs. That 

appearance does not necessarily suggest that the microsites are strong 

absorbers after creation. The same appearance would result for no 

absorption, but instead only scattering of light out of the collection solid 

angle of the microscopic objective. The microsites do scatter light 

strongly, as verified by observation with dark-field illumination. The 

morphology of fig. 1 is typical for all the crystalline materials that we 
1 E 

tested. ’ The morphology in the tested non- crystalline materials 
(silica and laser glasses) showed less distinct microsites, however. Such 
a difference is to be expected due to the lack of surrounding crystalline 
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a) 


b) 


c) 
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Fig. 1 Breakdown morphology in a NaCl sample produced by 
1.064 ym, nominally 30 psec laser pulses focused to an e"^ 
rms intensity radius of 3.3 ym. The scale markers denote 
a distance of 10 ym, and the pulses were incident from the 
left. The maximum rms electric field values reached at the 
focal point in Figs, la-c were 1, 20, and 52%, respectively, 
above the previously measured threshold. 
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medium. The crystalline backdrop provides a light scattering discontinuity 
between unmelted crystalline structure and melted material which re- 
solidifies presumably into a disordered state. 

Figs, la-c illustrate the damage morphology for three different 
input intensities. In fig. la, the maximum rms electric field reached at 
the focal point was 7. 4 ^ 1. 1 MV/cm (corresponding to 65. 6 kW), The 
breakdown threshold found previously^ for NaCl was 7. 3 i 1. 1 MV/cm, 
and therefore fig. la illustrates damage essentially at threshold. The 
damage is accordingly observed to consist of a small number of isolated 
scattering centers localized on the focal axis, each of which is the vestige 
of a laser-induced microplasma. Each microplasma began from a single, 
or perhaps a few, easily freed electrons along the focal axis during the 
passage of the laser pulse. Fig, lb Illustrates the damage in the same 
sample under the same focusing conditions due to an 86. 3 kW pulse with 
maximum focal point electric field of 8, 7 MV/cm, about 20% above 
threshold. A pulse of 127. 6 kW and 11. 1 MV/cm maximum field (52% 
above threshold) produced the damage shown in fig. Ic. 

The diameter t f the microsites in the figures is typically about 1 {xm. 

If we assume that the micrcplasma began with an electron located at point 
r' at the beginning of the pulse and model the growth of the plasma by a 
simple diffusion model, we may estimate the expected spatial extent of 
the light- scattering remains of such a microplasma to have a diameter 
fB Zjr - r'l of roughly 4 n/d^, Here T^ is the laser pulse duration 
(FWHM, 30 psec. ) and D is the electron diffusion coefficient, which is 
expected to have an average value for the high electric fields encountered 
here of roughly 1 to 10 cm /sec. These two numbers bracket the microsite 
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diameter by about 0. 2 and 0. 7 microns. This rough estimate is consistent 

with the observed diameter. It is evident that for laser pulses with 

duration a few times longer than the present 30 psec FWHM duration, the 

focal region would show continuous breakdown. Resolved microsites 

would not be observed Such continuous morphology was observed by 
6 7 

Yablonovitch, Fradin^^ , and others in damage studies with nano- 
second pulses. 

11 2 

Just at threshold (2. 2x10 W/cm maximum focal point intensity) 

the volumetric density of microsites from fig. la is approximately 

9-3 3 

4 X lO cm , corresponding to the observed linear density of 1, 5 x 10 cm 

11 -3 

The maximum density in fig. lb is about 2 x 10 cm , for a maximum 
focal point intensity level of 3. 1 x 10 W/cm . This density decreases 
slowly out to the extremities of the damage, where the maximum field 
attained was the threshold value and the density is approximately equal to 
that in fig, la. In fig, Ic one observes along the track the progression in 
density of microsites from approximately that of fig. la, at the extremities 
of the damage, through that of fig, lb, into a region where the damage 
along the focal axis has become continuous, except for short discontinuities 
caused or enhanced perhaps by diffractive shadowing. 

The damage sites just above threshold are statistically spaced and 
symmetrically distributed within their region of occurrence along the 
axis. There is initially no shadow effect due to absorption and scattering 
out of the incident beam. The light diffracts readily around the damage 
sites in the trailing edge of the pulse and the effective integrated intensity 
on axis is hardly diminished. At higher intensities the loss of incident 
intensity due to damage spots creates the asymmetry noticeable in fig, Ic. 
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The maximum intensity reached at the focal point in fig. Ic was 
11 2 

5, 0 X 10 W/cm , It is important to note that the electric field and 

intensity values discussed in this paper are maximum rms values reached 

only instantaneously as the temporally gaussian laser pulse passed through 

a given spatial point in the sample. 

8-15 

Much discussion has taken place in the recent literature con- 
cerning the number and origin of initiating electrons for the avalanche 
process, and on their influence in associated aspects such as the statistical 
nature of laser -induced breakdown. In the original work by Yablonovitch 
which established 10. 6 (xm optical breakdown as an avalanche process 
closely related to DC breakdown, the initial electron density in the 
avalanche relation 


N(t) 


N exp 
o 



Ti{E(t')j dt» 


( 1 ) 


had to be obtained essentially by order -of -magnitude guesswork. The 

resulting successes of the avalanche model have been due to the weak 

dependence of moat aspects on N^. However, statistical aspects of laser- 

13-15 

induced breakdown depend markedly on N^. The microsite density 

in figure 1 may be identified with the density N^. It is of course possible 
that more than one initiating electron leads to the production of one micro- 
site. Thus the derived value of N is approximate, and has the order of 
magnitude of 10 /cm , in agreement with past aummptions. The 

production mechanism of these electrons could be investigated in more 

12 

detail with our observation technique. It has been argued from electrical 

9 3 

conductivity properties that 10 free electrons per cm do not exist in 
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9 12 

insulating crystals, but that initiating eledrims are produced * at 1.06 gm 

promptly by two-photon ionisation of ground-state F centers which may 

J 1 - 3 

occur in conct^ntrations as high as 10 /cm , even in a clear crystal with 

-4 - 1 

linear absorption coefficient a ^ 10 cm at 1.06 pm. It is evident 
that further studies of this nature as a function of F-center concentration, 
temperature, ambient light, etc. , would be useful in fully understanding 
this aspect of laser -induced breakdown. Any color center densities must 
be implanted without concomitant introduction of metallic colloids which 
would confuse subsequent data interpretation. 

The damage microsites are clearly confined to a narrow region on 
the axis of the focal spot, where the laser intensity is maximum. The 
region grows both in length and in diameter, with the amount by which the 
maximum laser intensity exceeds the threshold value. The region has an 
ellipsoidal shape and is plotted as the dotted line in figure 2 for the case of 
ligure lb. It may be roughly interpreted as the volume element in which 
the laser intensity exceeds the threshold value. These data are testimony 
to the reproducibility of the laser pulses and their effects. 

The damage is due to heating of the lattice by resistive losses of the 
electrons in the plasma. Since the heat diffusion length during the pico- 
second pulse is much smaller than the electron diffusion length which 
determines the size of the microsite, the temperature rise in a microsite 
located at the point {r, z) may be calculated from 

2 t 

T(r,z,t) — - j T -o " N{r,z,t*)lE(r,z,t')l^ dt* , (2) 

mC(l 4 

where the laser electric field distribution is given by 


i 
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E / Z ( 

lE{r,z,t)[^ = 2 ^ ^ 2 ' expj -ry[p^(l-P{t)P“J)(l-z^2^p}] 

1+z 

(3) 

The electron density is given by eq. (1) where the ionization rate q may 

18 1 

be calculated from ^ 

log T] {E(r, z, t)} = 3. 58 log E(r, z,t) + 8, 62 , (4) 

for E between 1. 0 and 9. 0 in units of MV/cm. For in eq. (2) we use 

12 -3 

one electron within a microsite volume, or N =2x10 cm . In 

’ o 

expression (2) above, e and m are the free-electron charge and mass, 

respectively, C is the sample heat capacity per unit volume (1. 85 j/cm^*^C} 

for NaCl, o) is the laser angular frequency, and is the effective 

- 1 18 

electron collision time which we have taken here to be u) . In eq. (3) 
is equal to E(0, 0, 0) which is the maximum rms field reached at the 
focal point, occurring at t = 0. The e”^ intensity radius in the focal plane 
is p(= 3. 3 pm here), and z^^^ is the Rayleigh distance for the electric 
field, 2TTnp /K, or 98.5 microns here (where n is the linear refractive 
index). The e temporal intensity half- duration is denoted above by t, 
having values of 23. 2, 19* 4 and 27. 0 psec, respectively, for the pulses 

3 

which produced figs, la, b, and c, P(t) denotes the laser pulse power. 

1 

P^^ is the critical power for self -focusing, taken here to be 166 kW for 
NaCl. Using eqs. (1-4), we may compute ths three dimensional surface on 
which the maximum temperature reached just equals the melting temp- 
erature of NaCl (800*^C) and within which any volume element of microsite 
size would be melted, if it contained an initial electron. The result of 
this calculation is plotted as the solid line in fig. 2. The good agreement. 
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Fig. 2 comparison of the breakdown 
morphology perimeter of Fig. lb /in- 
dicated by the dashed lines, with 
an avalanche heating isotherm as 
computed with Eq. 2. The isotherm, 
the solid curve in the figure, indi- 
cates a temperature of 80Q ®C, the 
melting point of NaCl. The z axis 
is the focal axis. Note the differ- 
ent scale factors for the r and z 
axes. 
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especially as to the shape of the region with damage sites, is strong 
evidence that the avalanche ionization mechanism is responsible for 
the damage. 

The difference in size of the region with damage shown in figures la 

and lb, corresponding to a difference in input pulse power of only 32 

percent, cannot be explained on the basis of absorbing inclusions. These 

(invisible) inclusions would have to be uniformly distributed with a density 
11 -3 

of 10 cm , and would all have to be uniform in absorptivity and size. 

They would also have to vary in a systematic way to explain the observed 

characteristic threshold for each bulk material, ^ in particular the 

systematic compositional variation throughout the periodic table for the 

alkali halides. Since the specter of absorbing inclusions, influence of 

which was carefully assessed and eliminated in the present experiments, 

19 

was raised in a recent publication, the question of metallic Na inclusions 
in Nad is considered in further detail. The invisible inclusions are 
obviously smaller than the wavelength of blue light (< 0. 4 pm), but in order 
to melt a volume of this size to produce microscopically visible damage, 
they must have a diameter at least as large as about 0. 1 pm. This number 
was calculated by solving the transient heat conduction boundary value 

problem for the absorber-host system, using the Mie absorption cross 

20 11 3 

section and the incident intensity. If the sample contains 10 /cm of 

Na particles with a diameter larger than 0. 1 pm, one would of course 

find a linear attenuation coefficient many orders of magnitude larger than 

-4 - 1 

the observed attenuation coefficient of < 10 cm in Harshaw NaCl 
samples. 
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19 Z2 

Other experiments ’ also showed filamentary damage of up to 

1 mm track length and were carried out at power levels that allow a 

23 24 

moving self-focal spot ' to develop. 

In our experiments the effects of incipient self -focusing were 

1 25 

carefully corrected by the method proposed by Zverev and Pashlcov 

26 

and by Fradin. Our results add further strong evidence to that of 
1 Z 6 7 27 

earlier papers > » ’ J ' that the fundamental damage threshold is 
determined by the avalanche ionization process. 


The authors are indebted to Prof. Ij. Taylor of Harvard University 
for the use of microscopy equipment, and to Prof. F. Luty of the 
University of Utah and Prof, E. Yablonovitch of Harvard for helpful 
discussions. The skillful sample preparation of S. Maurici is greatly 
appreciated. 
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CHAPTER IV. 


PICOSECOND LASER-INDUCED BREAKDOWN AT 5321 8; 

OBSERVATION OF FREQUENCY-DEPENDENT BEHAVIOR* 

W. Lee Smith‘*‘|. J, H. Bechtel t N. Bloembergen 

Gordon McKay Laboratory, Harvard University 
Cambridge, Massachusetts 02138 

Laser- induced dielectric breakdown at 5321 8 has been 
investigated in six materials: KH 2 PO^, fused Si02, NaCl, CaF 2 # 

NaF, and LiF. The laser pulses were nominally 21 psec in dura- 
tion and were obtained by the frequency doubling of single 1.064 
Jim pulses from a mode-locked YAlG:Nd laser system. The fre- 
quency doubling was carried out in a temperature phase-matched 
CsH 2 AsO^ crystal, resulting in spatially smooth, reproducible, 
diffraction-limited, 5321 8 pulses. Strict attention is given 
to the complications introduced by self-focusing. The new thres- 
holds for breakdown at 5321 8 are compared with thresholds 
measured previously at 1.064 with the same laser system. 

This comparison, as a function of material band gap 3 ' q * charts 
experimentally for the first time the transition of the char- 
acter of laser-induced breakdown as it becomes strongly frequency- 
dependent. A frequency-dependent decrease in the threshold of 
~30S is observed for KH 2 pO^, for which three 5321 8 photons 
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exceed in energy. Frequency-dependent threshold increases 
of up to 445, varying smoothly with (fg, are observed for the 
other materials. These observations are discussed with regard 
to models for the processes involved in the breakdown, which are 
avalanche and multiphoton ionization. 


\ 
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A. Introduction 

In this paper we present experimental results which 
characterize for the first time the transition of laser-induced 
dielectric breakdown from the familiar DC limit into the regime 
where both avalanche ionization frequency dependence and multi- 
quantum effects are clearly observed. This study represents a 
systematic frequency extension of our previous 1.06 ym investi- 

1-3 

gations of picosecond dielectric breakdown in transparent solids. 

Laser-induced breakdown was reported in the literature early 
in the 1960's, However, the earnest study of DC breakdown has 

4 

taken place since the 1930 's. Indeed, it has been pointed out 
that laboratory observation of DC breakdown in glass was reported 
as early as 1799, by A. Van Marum in the first volume of Annalen 
der Physik . Van Marum 's efforts at producing high voltages were 
limited by breakdown through the walls of his storage jars. Be- 
ginning in 1931, experimental progress in DC breakdown threshold 

g 

measurement was achieved by von Hippel. Theoretical developments 
7 8 

by Frohlich beginning in 1937 and Seitz in 1947 proposed and 

described avalanche ionization as the basic physical mechanism 

involved in the DC case. In 1971 Yablonovitch performed laser- 
. 9 

induced breakdown experiments at 10.6 ym in which the similarity 
of the threshold between the infrared and earlier DC data, both 
in the trend with material composition and in absolute magnitudes, 
was compelling. Avalanche ionization was established as the 
intrinsic process whereby an infrared laser pulse is able to 
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induce breakdown in pure dielectric material. Further controlled 
laser experiments at 1.06 ym^^ again found breakdown thresholds 
not significantly different from those at DC and 10.6 ym, either 
in trend or magnitude. At the still shorter ruby wavelength of 
6943 S, a systematic increase over previous thresholds was 
reported. As will be discussed in a later section, the pattern 
of increase was unusual, however, and it now appears that the 
threshold increase in the low band gap materials may have been 
an artifact. Yet at least part of the large relative increase 
in the NaF threshold reported in that study was quite probably 
an observation of true behavior and therefore the first indication 
of frequency dependence in avalanche ionization breakdown processes. 
An increased resistance to surface damage in several materials 
at 6943 R over that at 1.064 ym was also reported^^ in 1973, 

Because that investigation involved rather different experimental 
conditions (surface breakdown, multiple pulse irradiation of 
the same sample volume) than the present work, we will not consider 
it further. 

At wavelengths shorter than 6943 S several experiments 

13“18 

have been performed to date. However, in each of these 

experiments self-focusing was present. No attempt was rade to 
extract the self-focusing influence in order to measure a break- 
down threshold. At best, such experiments yield a measure of 
the critical power for self-focusing. Further interpretive 
difficulties were caused by the use of multimode lasers. 
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entire mode-locJced pulse trains, or hajcmonic generation 

crystals allowing walk-off- Under such conditions reliable 

breakdown thresholds are rarely obtainable. 

Thus, up to electromagnetic angular frequencies w of 
15 

almost 3x10 Hz, the basic behavior of bulk dielectric break- 
down has been observed to be quite similar to the behavior at 
zero Hertz. The concise reasons for this steadfast behavior 
over such a remarkable frequency range are two-fold; i) The 
hot electron collision time (t^) in dielectrics is very fast 
(<10~^^ sec) . ii) In experiments up to 6943 S, the photon 
energies were “'ot a sufficiently sizable fraction of the 
material band gap,(^. In Fig. 1 we illustrate the schematic 

frequency dependence of the breakdown threshold Eg, as has 

19 

been anticipated for some time. The low frequency, pure 
avalanche region (solid segment) is the range of investigation 

9 

prior to this study. As discussed elsewhere, the relation 


£3(01) 


= E 


DC 
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„ V- 


1 + 


2 2 


( 1 ) 


has been used to predict the frequency dependence of laser- 
induced breakdown in the region of departure from the DC limit. 
In Pig. 1 the rolid segment and its dotted extension have the 
form of Eg. 1. The aforementioned increase in the NaP thres- 
hold cit 6943 S indicates that hT^^ ^'^ruby NaP, 

and roughly indicates the ht”^ position on the axis of Fig. 

o 
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Fig. 1 Schematic fregiiency dependence of the two intrinsic 
physical processes involved in laser- induced breakdown: 
avalanche ionization and multiphoton absorption ionization. 

denotes the optical band gap of the material and hT“^ 
the energy associated with the electron collision time f . . 
The solid curve is the region investigated prior to this‘s 
study. 'T’he present work investigates the region located by 
the arrow. 

Fig. IV. 1 
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The dashed curve outlines the frequency dependence of the 

20 

Keldysh tunneling mechanism. As the photon energy becomes 
sizable compared to the material band gap, this tunneling 
mechanism is described as multiphoton ionization. At some 
point in this transition the multiphoton ionization probability 
becomes large enough to effectively assist the breakdown process, 
still classifiable mainly as avalanche ionization. Several 
papers recently have considered aspects of this combination 

process. One iiranediate consequence of the multiphoton ioniza- 
tion assistance is the ready supply of initial electrons with 
which the breaJcdown process may proceed. For sufficiently lower 
frequencies, the productions of such electrons is a subtle 
matter. As the photon energy exceeds <^q/ 3, we expect the thres- 
hold to reverse its trend as the arrow indicates in Fig, 1. 

As the photon energy is increased beyond gross two-photon 

absorption (TPA) will occur. At such frequencies, breakdown 
experiments utilizing an external focusing lens will be com- 
plicated considerably by TPA energy depletion and distortion 
of the laser pulse in the medium before the lens focal region. 

One should note the manner in which the two curves in 
Fig, 1 will shift horizontally as the material parameters 
and hT^^ are changed and hw is held fixed. Because hi^^ is 
likely to be less sensitive to than is the threshold for 
strong multiphoton absorption, the most demonstrative effect 
on Fig, 1 from changing to a matarial with larger would be 
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the shift of the dashed curve to the right and creation of 
a larger hump in the composite threshold trajectory. There- 
fore a comparison of breakdown thresholds at two different 
laser frequencies, both in the region of the transitional 
hump, for a series of materials and band gap, would provide 
a test of the surmised behavior which is illustrated in Fig. 1. 
From our picture, one would expect from such a comparison to 

( 1)2 ^ ^ 1 , as a function of material band gap. We will present 
data later in thj.s paper which depicts just such a systematic 
increase. 

A recent review of dielectric breakdown induced by laser 

pulses is found in Ref. 19. Recent theoretical papers which 

will serve as well to introduce the reader to segments of the 

literature are found in Ref. 24 by Holway and Fradin and Ref. 

25 by Sparks and Duthler. Numerous papers concerning all 

aspects of laser-induced dielectric breakdown appear in National 

2 6 

Bureau of Standards Special Publications. The 1973 edition 

27 

of O' Dwyer’s book may be consulted for information on DC, 

as well as high frequency, breakdown phenomena. 

Section B outlines the theory necessary for the evaluation 

of our experiments. Section C relates the experimental details, 

and we present results in Section D. Comparison of our exper- 

mental results with theories on multiphoton contributions to 

4 ,21 

breakdown are given in Section E. 
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B. Theory 

The objective of our experiment is to obtain the 
high electric fields requisite for breakdown, while at 
the same time minimizing the input pulse power and thus 
the effects of self-focusing. This we do by tightly 
focusing the laser pulses into the sample with short focal 
length lenses. In the absence of any self-focusing, we 
know that a diffraction-limited light pulse with a Gaussian 
spatial intensity distribution of radius p j^I(r) = exp 
-(r/p)^ j , focused by a lens of focal length f, will produce 
an intensity distribution in the focal plane of 

If(r) = (P/A) . (2) 


Here P is the pulse peak power, defined as 


P 


2tt 


CA 

r 

dr r I(r) 

J 

0 


(3) 


2 

where A = ttw , and w = Af / 2'trp, if A is the light wavelength 
in air. In addition, we know that if the power P in the light 
pulse exceeds some fraction of the critical power value known 
as P^, defined as 


= cA^ / 32Ti^n2 


(4) 
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for a Gaussian pulse, the above focusing behavior is 

modified by self-focusing. The material parameter which 

determines is ii 2 f the nonlinear refractive index. If 

E is the rms field amplitude and n is the total refractive 
. . . 2 

index, TI 2 is defined via n = n^ + n 2 E +... The equations 

above may be altered to account for self-focusing with the 

help of the analytical "aberrationless " or "constant shape" 

28—3 0 

(CS) approximate theory if P is "less than" P^. In 

that case, Eq. 2 becomes 

2 

If(r) = (P/A) d-P/Pj^}"^ . (5) 

Depending on whether P^ is known for a given material of 
interest, Eq. 5 may be used in two different ways for break- 
down threshold measurements. 

If U 2 and thus P^^ are not known, then Eq. 5 may be 

used to measure both the breakdown threshold electric field 

31—33 

Eg and P^, simultaneously, as was done by others and 

1 — 2 

by ourselves in earlier 1.064 pm studies. In that situation, 

lenses of different focal area are used to produce break- 
down at the focpl point, each lens requiring a different 

input power P^. One may then plot the various reciprocal 

-1 -1 
input powers versus the respective focal areas A^^ . As 

long as all are much "less than" P^ (an unknown) , the points 

will define a straight line, as indicated by Eq. 5. Fig. 2 


EEPRODUcmiLiTy or the 

0RIG5NAL PAGE IS POOR 



CONSTANT SHAPE 
APPROXIMATION ' 


18 


‘ A • * > « I 


•BREAKDOWN 
CURVE.: 


mm 


Fig. 2 Illustration of the slope-intercept breakdown construction, as given 
by Eqs. 5 and 6, wherein one plots for a given lens the reciprocal pulse input 
power (p-l) versus reciprocal focal area (A-^) . The dashed curve is defined by 
the breakdoxvn properties of the material, and the solxd line is from, the constant 
shape approximation, and Po are the "paraxial'* and the "catastrophic '* critical 
powep, resp. See Sec. B for discussion of the two curves. is the minimum 

focal area obtainable with a given wavelength. 
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indicates these features schematically. There the slope of 
the solid line is 1^(0)"^ (and equals the reciprocal break- 
down intensity at threshold) and is obtained using the 

r 1 1/2 

relation Eg = 1^(0) / n^e^c , where n^ is the linear 

refractive index, is the free-space permittivity, and c 
is the vacuum speed of light. The intercept of the line with 
the P ^ axis yields If# however, TI 2 is known, then the 

threshold power measurement from only one lens suffices to 
determine Eg from Eq, 5. It should be stated that the slope- 
intercept, controlled breakdown method discussed here and 
previously is an intrinsically less accurate way to measure 
than some other available schemes. It is, neverthe- 
less, in the face of unknown n^ values, the only known way 
to obtain bulk breakdown threshold measurements in the pico- 
second pulse regime which are corrected for self-focusing - 

Aside from the paucity of accurate available n 2 values, 
another situation exists which thwarts the extension of 
picosecond breakdown measurements to shorter optical (UV) 
frequencies. All the above discussion was based on input 
pulse powers P being "less than" P^, the critical power for 
the onset of self -focusing by paraxial rays in an unfocused 
geometry. The problem lies in our inability to ascribe more 
precision to the phrase "less than" in the previous sentence. 
It is not known quantitatively how, in the focused geometry 
(where an external lens focuses the light pulse into the 
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material sample ) , the constant shape approximation breaks 

down as P approaches and, moreover, exceeds P^. Fig. 2 

illustrates our present understanding of this behavior. 

The cross-hatched area is the region of validity of the 

CS approximation. The dashed curve is the locus of (A 
-1 

P ) points defined by the breakdown threshold of a particular 

material. The straight line, which is asymtotic to the dashed 

curve in the low power limit, is defined by Eq. 5 and is the 

1-2 

construction used in our previous papers to measure Pj^. 

The dotted area in the figure is the region within which the 

self-focusing behavior has departed from the CS description. 

The exact functional form of this departure, and thus of the 

dashed breakdown curve in the dotted region, is not known. 

Therefore, the exact limiting power level (indicated by Pj^^ 

in Fig. 2) , beyond which one should not apply the constant- 

shape approximation to deduce breakdown thresholds, is not 

28 

known. It is indicated by careful numerical calculation 
that the power in a Gaussian pulse must exceed P^^ by a factor 
of 3.7 before a catastrophic self-focus (a singularity in 
axial intensity) will occur, in either externally focused 
or unfocused geometry. This intensity singularity manifests 
itself by breakdown and other processes, and represents the 
horizontal segment of the breakdown curve in Fig. 2. This 
multiple of P^^ leading to the singularity is denoted in the 
literature by p 2 * For pulse powers very near to but less 
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than P^f the maximum axial intensity may be expressed then 

2 8 

by Eq. 5 with P^ replaced by P 2 . And in fact, the maximum 
axial intensity reached at any input power level might be 
expressed by the relation 

^MAX " [^" . (6) 

Here ^(P) is presently unknown and must be the result of 
numerical calculation, but will behave monotonically between 
the following limits: 

<|){P) ->-1 for P << P^ , 

(f)(P) •> 3.7 as P "> p 2 » 

2 8 

Data published by Marburger indicates that efforts may 

be made to compute the function 0(P). Such new results, 

in a form applicable to any focusing geometry, would be of 

great use for experiments of many types, expecially with 

2 

higher laser frequencies considering the restricting X 
dependence of P^. Knov7ledge of (|)(P) would, for example, 
remove our present dependence on the CS approximation and 
allow more than 4 times more power to be used to produce 
calculable intensities for breakdown studies such as this 
one. In that case the segment of the breakdown curve 
within the dotted area and down to the P 2 ^ level would be 
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useable. Similar applications could be made to studies of 
absolute coefficients for UV two-photon absorption in solids 
and liquids, etc. The data presented in this paper will be 
tabulated so that the breakdown threshold values may be 
recalculated by the reader in the future, should data on 
the <ti(P) fi;i> .* ion appear. 

C. E:cperiment 

1 . Apparatus 

The experimental arrangement is shown in Fig. 3 and is 

in large part the same as used in Refs- I- 3. The YAlG:Nd 

laser oscillator is mode- locked by saturable absorber (Kodak 

9860). It produces a train of TEM^^ transverse mode, 1.064 

um pulses of 30 psec average duration. One of these pulses 

is selected for use by a KH 2 PO^ Pockels cell which is driven 

by a laser-triggered spa-k gap. The selected pulse is then 

monitored for proper switch-out behavior on a Tektronix 519 

oscilloscope (CRO 1) . Photodiodes P2 and P3 measure the energy 

content in two samples from the main pulse, one at 1.064 pm 

and the other converted to 5321 S by a KH 2 PO^ crystal. From 

these two signals we obtain the temporal duration of each 

40 

individual picosecond laser pulse. With this temporal 
characterization, the main part of the infrared pulse is 
converted to a 5321 S pulse by harmonic generation cell 
(HG2) containing a cesium dihydrogen arsenate crystal. The 
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crystal is maintained in an oven at the proper temperature 

for noncritical or 90° phase-matching. Two characteristics 

of the green pulse , crucial for this experiment/ accrue from 
o 41 

90 phase-matching. First/ no spatial distortion of the 
green pulse occurs by "walk-off /“ as occurs with angle-tuned 
phase-matching. The 5321 S pulse is cylindrically symmetric, 
without spatial noise, and Gaussian in shape. Fig. 4 is a 
tracing of an oscillogram from a Reticon linear photodiode 
arr-y (Reticon Corp., Sunnyvale, California) of a typical 
green pulse profile. The waveform was smoothed by no more 
than the width of the tracing pen. Such a distortion-free 
profile is required for accurate calculation of the electric 
field distribution of a laser pulse in the focal region of 
a lens. Second, this harmonic generator produces, for an 
unlimited length of time, green pulses with the same spatial 
profile. Such reproducible behavior, over a period of several 
hours and from day to day, is difficult, if not impossible, 
to obtain from angle phase-matched crystals. A serious attempt 
was made by us, prior to the present experiment, to accurately 
measure breakdown thresholds with 5321 2 . pulses generated in 
an angle-tuned KH 2 PC^ crystal. The effort was not satisfactory 
for the reasons outlined above. The same set of difficulties 
will thwart efforts to study breakdown at the third-harmonic 
frequency of the Nd:YAlG laser, unless materials are found to 
provide 90° phase-matching for 3547 2 generation. 
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Pig, 4 Oscillogram tracing of the spatial 
intensity profile of the 5321 A pulses pro- 
duced by the CSH 2 ASO 4 crystal. The oscil- 
logram was recorded with a commercial linear 
array of 256 photodiodes with 0.05 mm resolu- 
tion. 
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The 5321 8 pulse is separated from the infrared pulse by 

a prism, oriented to maintain the circular cross section of 

the green pulse. Photodiode P4 monitors the absolute energy 

content of the green pulse just before the focusing lens L. 

Photodiode P4 was calibrated with two Eppley Laboratory 

(Newport, Rhode Island) thermopiles, themselves calibrated 

against National Bureau of Standards sources. By using the 

energy content signal from P4 and the individual pulse 

duration which we derive from P2 and P3 signals, we obtain 

the peak power in each individual pulse. The measured full 

angle of divergence of the harmonic pulses is less than 0.9 

mrad. By measuring the cross-sectional area of the pulses at 

the lens position with the linear photodiode array, and using 

the design parameters of the lens, we may calculate the area 

in the focal plane of the lens.^ The measurements which we 

report here were made with a 1.0 inch focal length corrected 

2 

lens, with a calculated focal area (defined as ttw , where w 

-1 2 
is the e radius of intensity) of 14.8 ym . This figure 

includes contributions from spherical aberration of the lens 

and the planar sample surface,^ both of which amounted to a 

calculated correction to the area of less than 2%. This 

calculated focal area was checked experimentally by a scan- 

ning aperture technique. To within the limiting accuracy of 

the technique, ±20%, the experimental and calculated areas 

agreed . 
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We have employed an average full-width-at-half-intensity 

duration of 21 psec for our 5321 S pulses. This duration is 

shorter than the 30 psec parent pulse duration by 1/ /2, as 

given by the theory for the ideal frequency doubling of 

Gaussian pulses in the low conversion efficiency regime. 

Final individual pulse durations were calculated using this 

40 

average duration value ^ the tA technique, and a computer 
program to account for the temporal stretching of the pulse 
due to saturation of the conversion efficiency. This last 
correction is discussed a.n Section D.l. 

2. Procedure 

The procedure for each measurement began with the 
alignment of a sample (dimensions typically 6 mm x 20 mm X 
30 mm) in the x-y translating holder so that the pulse 
propagation direction was normal to the sample face. Then 
the lens itself was aligned with a reproducible procedure, 
discussed previously, which was made easier by the visibility 
of the light in the present experiment. Next the sample was 
checked for unacceptable surface polishing roughness, a check 
which was not possible with invisible 1.064 urn light. Shots 
were fired through the sample, below the breakdown threshold, 
and the transmitted light was observed on a card. If the 
transmitted light pattern showed a contribution due to dif- 
fraction from polishing scratches, the sample was either 
repolished until no surface scratch diffraction was observed, 



rv-21 


or the sample was deleted from the experiment if sufficient 
smoothness was unobtainable. The polishing was performed 
by an experienced technician using synthetic lap material 
and alumina abrasive particles down to 0.05 ym diameter. 

Several materials were unfortunately elimiated due to this 
problem - the relatively soft, small band-gap materials 
including NaBr, KBr, and KI. It now appears that such 
surface distortion may have been the cause of the hereto- 
fore unexplained rise in the breakdown thresholds, relative 
to NaCl, of KBr, KI at 1.064 ^ and 0.694 over their values 

at DC and 10,6 ym. The absence of this problem in another 

in 9 

study ^ at 1.06 ym is not understood. At 10.6 yra, however, 

such submicron roughness would have almost no distorting effect. 

After a proper sample and the lens were aligned, a series 
of 50-100 shots were fired into sites with recorded location 
in the sample, one shot per site. The pulse power was increased 
throughout the series, by altering the Pockels cell voltage 
which changes the position of switch-out from the oscillator 
train. In contrast to Nd: glass mode-locked lasers, Nd:YAG 
laser pulses do not vary in duration as a function of pulse 
position in the train,. For each shot, signals from detectors 
P2, P3, and P4 were recorded and visual observation of micro- 
spark formation was used to indicate whether breakdown occured 
on each shot. Such visual sightings were not used to determine 
the breakdown threshold, however. If one wishes to meaningfully 
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compare thresholds measured at different light wavelengths, 
then the threshold observation mechanism cannot itself be 
wave length -dependent. The ocular observation of spark 
production (with its attendant scattering of subsequently 
arriving laser light) would be particularly subjective and 
wavelength “dependent between 1.064 ym and 5321 S. So the 
threshold here was defined by the same criterion as used in 
our previous 1.064 ym study, namely as the minimiim electric 
field value, reached at the focal point and at the peak of 
the pulse, which produced microscopically observable, per- 
manent damage in the material. The microscopic observation 
of the matrix of sites was performed after completion of 

the laser irradiation, using observation conditions identical 

3 

to those in our 1.064 ym study : Leitz Ortholux microscope, 

340 power total magnification, 0.4 NA long working distance 

objective, blue light field illumination, '-7000 2 resolution. 

Although the visual sightings of spark formation were not 

appropriate for identification of breakdown thresholds, the 

visual thresholds obtained with the 1.0 inch focal length 

lens, and also with 1.5 in. and 3.0 in. focal length lenses, 

were recorded for each material. The calculated focal areas 

2 2 

of the two latter lenses were 32,9 ym and 131.1 ym , respec- 
tively. These visual thresholds were used to investigate 
self-focvjirg behavior for pulses having powers between 
and Discussion of that behavior is found later in 
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Section D. 2. 

At least three different experimental runs were made on 
each tested material with each lens to ascertain the thres- 
hold reproducibility. 

D. Results 

1. Dielectric Breakdown Threshold Data 

In this section we will present and discuss new results 
for dielectric breakdown thresholds at 5321 S, induced by our 
nominally 21 psec pulses. 

Several general characteristics of this damage study were 

observed to be the same as those encountered in Refs. 1-3. 

. . 23 

The thresholds were again noted to be sharply defined. That 

is, our resolution of the thresholds was limited by our exper- 
imental uncertainty and reproducibility, which was observed 
to be ±155S. Again, only rare evidence of the presence or 
influence of inclusions was observed under the microscope 
CilS of the sites) . Again, the results did not vary from 
site to site within a sample, or among different samples of 
the same material, by amounts larger than our experimental 
uncertainty. Fig. 5 exhibits the morphology in two of the 
sample materials, KH^PO^ and LiF. These breakdown sites were 
photographed with blue light on a Leitz microscope at 250 
power magnification and ~0,8 ym resolution. The focal spot 
radius to the 1/e point of the electric field (/2 w) is 
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(b) 


Fig. 5 Breakdown morphology in (a) KH 2 PO 4 and (b) LiF, 
induced by 0.532 yn picosecond laser pulses focused to 
an e"^ rms electric field focal radius of 3.1 ym. The 
scale markers denote a distance of 10 ym, and the pulses 
were incident from the left. The maximum rms electric 
fields reached at the focal point for photographs (a) and 
(b) were 14% and 40% above the respective thresholds for 
the two materials. The mottled background appearance in 
(b) is immaterial. 


Fig. IV. 5 
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3.1 14m/ and the Rayleigh distance for the electric field 
(given by nn^w /A ) is approximately 85 ym. The scale 
markers in the figure represent a 10 ym distance. The laser 
pulses which caused the damage of Fig. 5a and 5b had max- 
imixn focal point' electric fields of approximately 14*o and 
405 above the respective thresholds for the two materials. 

The typical microscopically observed damage morphology of 

3 

Fig. 5 is similar to that reported earlier at 1.064 ym- 
Except for rare impurity sightings / damage always occurred 
within a needle shaped volume along the propagation axis at 
the focal point, the dimensions of which were only ~1 ym 
by -^10 ym at threshold. The perimeter of this damage was 
nonlinearly correlated with the incident electric field 
strength/ again as encountered at 1-064 ym. Within the 
damage perimeter occurred nucleation sites of ~1 ym diameter, 
at a density of 10 cm at threshold. Ko major differences 
were observed, between 1.064 ym and 0.532 ym morphology, which 
indicates that the final stages of the damage process are 
still due to avalanche ionization. A major difference in 
morphology would indicate massive or dominant influence from 
multiphoton ionization processes. 

Table 1 lists several material parameters for the six 
tested dielectrics: KH^PO^, Si02 (fused), NaCl, CaF 2 , NaF, 

and LiF. All the crystalline (monocrystalline) samples were 
cut so that the laser pulses propagated along the [001] axis. 
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TABLE IV. 1 

Test Material Parameters 


Material 

(eV) 


KH^PO^ 

6.92 

1.5316 

Si02 

7,8 

1.4608 

NaCl 

8.7 

1.5491 

CaF2 

10,0 

1.4354 

NaF 

11,0 

1.3269 

LiF 

12.0 

1.3937 


"2 

Cio'^^esu) (kwl 

I ■ ■■ ' III 

1.00 ± 0.30^^^ 267 

1.00 ± 0.08 281 

4.20 ± 0.84^'^^ 64 

1.13 ± 0.34^*^^ 237 

0.94 ± 0.47^®^ 283 

0.92 ± 0.16^^^ 


(kW) 

988 

1040 

238 

877 

1047 


a . Pef • 39 

b. Ref. 35 

c. Ref. 34 

d. Ref. 36 

e. Chapter III, Part A of this thesis. 


290 


1073 
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with the electric vector along the [100] axis. The KH 2 PO^ 

samples were purchased from Cleveland Crystal Co. (Cleveland, 

Ohio ) , while the other crystalline materials were purchased 

from Harshaw Chemical Co. (Solon, Ohio) . The fused Si02 

samples are Suprasil I, purchased from Amersil, Inc. (Hillside, 

N. J.). Column 2 of Table I lists band gap (^^) values, all 

of which (except that for KH 2 P 0 ^) were obtained from standard 

reference books or manufacturer's data. For KH2P0^, we were 

unable to find either or optical absorption data in the 

o 43 

ultraviolet beyond 2100 A, which indicates a disturbing 

lack of knowledge of impurity absorption behavior of this 

important material in the UV. An ultraviolet spectrogram 

44 

of our own samples indicates strong absorption beginning 

at about 1800 2 ., or 6.92 eV. For convenience, column 3 

lists linear refractive index data, which is needed to 

convert nonlinear refractive index data (column 4) to third 

( 3 ' 

order nonlinear optical susceptibility (x ') data, for 
comparison with data published elsewhere |^n 2 = 127in~^ • 

. For KH„P0. and CaF,, ordinary ray n 
values are appropriate. In coliiran 4, we list the most accu- 
rate n 2 values available from the literature for our materi- 
als, according to their quoted uncertainties which range 
between ± 8 % and ±50%. The last two columns tabulate the 
critical powers and ? 2 , which differ by a factor of 3.7 
as discussed earlier. 
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In Table II we present the new breakdown threshold 

results at 3321 S along with comparison data at 1.064 ym. 

The second column lists values of input power P at breakdown 

threshold. The ratio of these power thresholds to the 

respective critical power values from Table I are listed 

in the third column. The aforementioned focal area of 14.8 
2 

pm was corrected slightly (~95) for each material to account 

_ , 45 

tor saturatxon of the conversion efficiency of 1,064 ym 
light to the second harmonic. These power dependent correction 
factors were computed by numerical solution of the harmonic 
generation equations including pump beam depletion. The 
actual focal area values used are listed in column 4 of Table 
II. This spatial broadening (the profile remains Gaussian 
to a good approximation) due to conversion saturation also 
lengthens the pulse temporally. All power values in this 
paper have been corrected by the appropriate small factors 
(^5^i) , obtained from the same numerical computation. 

The values of P, P^, A, and n^ from the tables suffice 
to determine using Eq. 5. The new values so obtained 
from our studies for the breakdown rms electric field strength, 
at 5321 S with nominally 21 psec pulses, are denoted E^(21 
psec) and listed in column 5 of Table il. The superscript g 
denotes the green wavelength. The thresholds range from a 
low of 12.4 MV/cm for NaCl to a high of 26.5 MV/cm for LiF. 

The discussion of the uncertainty in the E^(21 psec) 
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Dielectric Breakdown Threshold Results 


• * 

Material P ^/^l ^ e|{ 21 psec) Eg^(30 psec) ^ Eg^(21 psec) 



(kW) 



(MV/cm) 

(MV/cra) 

[MV/cm) 


KH2P04 

151.2 

0.57 

16.0 

23.4 

22.3 

33.6 

0.70 

Si02 

129.0 

0.46 

15.9 

19,0 

11.7 

17,6 

l.OS 

NaCl 

38.4 

0.60 

15.1 

12,4 

7.3 

11.0 

1.13 

CaF^ 

146.0 

0.62 

15.9 

25.2 

14.4 

21.7 

1.16 

NaF 

126.4 

0.45 

15.8 

19.4 

10.8 

16.3 

1,19 

LiP 

171.3 

0.59 

16.1 

26.5 

12.2 

18.4 

1.44 


^ Table 2, Chapter III of this thesis. 

^ The quantity A is defined as Eg (21 psec) / Eg^^(21 psec). 
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values of Table II is complicated by the unavailability of 
<{){P) data, our experimental conditions allow us to know the 
/p/a factor in the expression for Eg (Eq. 6) to an estimated 
±155. The absence of 4>(P) data has no effect on the repro- 
ducibility of our thresholds, of course, and our experimentally 
observed reproducibility limit of '^±155 reflects that fact. 

The absence of ^{P) data does, however, compound the uncer- 
tainty in the absolute numbers which are listed in Table II 
for Eg (21 psec) . In coliimn three of Table II, the P/Pj^ 
values range from 0.45 to 0.62, and it is plausible that they 
are too large to permit the CS approximation to yield high 
accuracy for the breakdown threshold electric fields from 
Eq. 5. For that reason we will later consider the full range 
of possible modification of our results, which in the worst 
case would be obtained by using denominator of Eg. 

5 in place of P^. For the present, however, we will continue 
to discuss Eg values obtained with P^ in Eq. 5 (4>(P) = 1) in 
the CS approximation. 

Our greatest interest lies in determining what effect 
the doubling of the light frequency has had on these thres- 
holds, by comparing them with the results of our previous study 
1 

at 1.064 urn. Column 6 of Table II lists the thresholds from 

T V 

that infrared study [denoted Eg (30 psec)]. As indicated by 
the notation, the ir data was taken with nominally 30 psec 
duration pulses, and it is desirable to separate out the pulse 
duration dependent changes in these two sets of data so that 
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we may observe changes due to light frequency, or photon 
energy, only. In order to convert the Eg (30 psec) data 
to Eg (21 psec) data, we require the pulse duration depen- 
dence in the picosecond range of the breakdown threshold for 

the six materials tested. Unfortunately such data is known 
1 

only for NaCl and, to a lesser degree, for five other 
alkali-halides. We are restricted, therefore, to an approx- 
imate treatment for removal of the pulse duration dependence, 
based on NaCl data only. From Fig. 6 of Ref. 1, one obtains 
that the ratio of Eg (21 psec) to Eg (30 psec) for NaCl is 
about 11.0 / 7.3 = 1.51. We have used this number to convert 
the entries of column 6, Table II, to those of column 7. It 
is, of course, unlikely that the threshold dependence on pulse 
duration would be exactly the same for the dissimilar col- 
lection of tested materials. However, over such a small range 
of pulse duration (30 psec to 21 psec) , the use of a single 
renomalization factor is reasonable. In that case, the change 
in breakdown thresholds due to frequency alone, denoted A, 
is the ratio Eg ( 21 psec) / Eg (21 psec) , and is listed in 
the final column of Table II. As discussed in Sec. A, the 
systematic variation of A with optical band gap should indicate 
the presence, if any, of frequency-dependent or quantum effects. 
Such systematic variation has never before been seen in experi- 
ments on laser-induced breakdown. Pig. 6 illustrates this 
point. In that illustration we have ordered the alkali-halides 
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Fig. 6 Illustration of the absence of multiphoton effects in breakdown threshold 
data for wavelengths down to 1.064 pm. The ratio of thresholds at the Nd frequency'^ 
to that at zero frequency® is plotted versus alkali-halide band gap. No systematic 
pattern is evident, in contrast with the higher frequency data of Fig, 7, 
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according to band gap. Plotted along the vertical axis is 
the ratio of breakdown thresholds from two previous experi- 
ments, one at 1.06 pim with '^5.0 nsec laser pulses [Eg^(5 

g 

nsec)] and the other at DC [E^ (DC)]. No systematic 

behavior of this ratio is observed, and the same is true 

for any other pair of laser-induced breakdown experiments 

reported thus far. The basis for this frequency-independent 

. 19 

behavior up to 1.06 iim has been discussed previously and 
is due largely to the twr jonditions mentioned in Sec. A con- 
cerning the brevity of the hot-electron collision time and the 
size of the band gap in these solids . Even the smallest band 
gap of 6.2 eV for the tested alkali-halides in the previous 
1,06 ym experiments would require at least six photons for the 
occurrence of direct multi-photon transitions, and such tran- 
sitions are vanishingly small even at breakdown intensity 
levels. At the ruby laser frequency (1.79 eV photon energy), 
four photons are needed to directly cross a 6-2 eV band gap 
and the probability is still small. At a 1,79 eV or 1.17 eV 
photon energy, it is probable that the only important role in 
the breakdown process in these insulators played by multiphoton 

transitions is the supply of avalanche initiating electrons by 

21 

the two-photon ionization of defect centers. As mentioned 
earlier, from our own work in the present and previous^ studies, 
it appears now that the unexplained increases in the thresholds 
(relative to that of NaCl) encountered by Fradin at 6943 2 
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and by us at 1.064 |im in soft, low band gap (^8 eV) alkali- 
halides may have been an artifact due to the undetected sample 
surface roughness mentioned earlier. Fradin also reported a 
large increase in the breakdown threshold of NaP over that of 
NaCl at 6943 S, and it is likely that most of that increase is 
a true observation, indicating the condition ^ 1 as discussed 
in Sec. A. This observation indeed indicated the departure of 
the breakdown process from the DC limit at 6943 S. Up to the 
present 5321 S. experiment, the NaF observation of Fradin and 
Bass was the only experimental datum point available for 
checking the predictions of frequency dependence siimraariaed 
in Pig . 1 . 

Now let us discuss the frequency-dependent threshold 
changes encountered in this study at 5321 2 (2.33 eV) . We 
plot in Fig. 7 the values of A listed in Table II versus 
material band gap Because of the specific uncertainty 

introduced by the lack of ij)(P} data, we have plotted in 
addition the A values obtained by using in each case in 
Eq. 5, rather than Due to the closeness of the 

values for the six materials (see Table I ) , it is probable 
that the true data points, were ij)(P) known, would all lie the 
same relative distance between the circle points and the 
solid points. That is, the adjustment of the circle points 
to their true location in Fig. 7 due to full knowledge of 
0(P) would probably result in only a small, nearly constant 



o-^(P) = 1.0 
• -^(P) = 3.7 




KH2PO4 SiOg NaCI CaF2 NaF LiF 


Fig. 7 Illustration of the band gap dependent behavior of the threshold. change. A, 
observed in the present experiment. Here A is the ratio E§(21 psec) / Eg^(21 psec) 
from Table II, plotted versus band gap (^q) for the six test materials. For each 
material, the open circle point and the dotted point result from the use of 4>(P)-1 
and (J)(P)=3.7 in Eq. 6, resp., to determine E^(21 psec). The discussion of the trend 
of this data is found in Sec, D. 
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lowering of tlie circle points. (Nevertheless, all information 
necessary to make such an adjustment in the future is given 
in this paper.) 

The variation of A, the frequency-dependent change in the 

5321 S and 1.064 ym breakdown thresholds, with material band 

gap in Fig. 7 is unmistakable. A A ordinate of 1.0 indicates 

no change in threshold on doubling the light frequency. The 

exact position of 1.0 on the vertical axis is not known due 

i^T 

to the uncertainty introduced by the Eg ( 21 psec) renormal- 
ization. Nevertheless, the data of Pig. 7 indicates that the 
band gap value for a A of unity is approximately 7.5 eV. The 
energy of 3 photons of 5321 S light is 7.0 eV. Thus, for 
materials with < 7 eV, the influence of multiphoton tran- 
sitions dominates over the competing effect of the light fre- 
quency exceeding the avalanche collision time of hot electrons: 
the result is A < 1. For KH2P0^, we obtain a A of 0.70 which 
is a frequency-dependent decrease of 30S. For materials with 
> 7 eV, 4 or higher-order multiphoton transitions (neglecting 
defect populations) are less effective against the cot increase 
at 5321 8, and A values exceeding unity are the result. A 
maximum A value of 1,44 is registered for LiF, a material 
(<^_ ss 12.0 eV) unaffected by interband multiphoton transitions 
at 5321 S. 

20 

It is perhaps instructive to consider the Keldysh model 
for an indication of the relative strengths of multiphoton 
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transition effects in our materials. The Keldysh model 
allows one to estimate for a particular field strength, the 
transition rate per unit of volume and time for the promotion 
of electrons from the valence band directly to the conduction 
band via the lowest-order , energetically sufficient, multi- 
photon transition process. We have used the free electron 
mass, in place of the unknown effective mass, the band gap 
values listed in Table I, and the breakdown threshold field 
strengths from Table II [E^(21 psec)l to calculate the 
Keldysh rates. Here, because we wish only to consider the 
rough variation of this rate for our six materials, such a 
single field strength estimation will suffice in lieu of a 

detailed computation. The Keldysh rate for KH 2 pO^, at 23.4 

28 3 

MV/cm with ho) = 2.33 eV, is large (~7 x 10 / cm -sec) due 

to the fact that 3hw slightly exceeds the band gap. As we 
go to SiO^ at a field strength of 19.0 MV/cm the rate drops 
by a factor of over 660 as the lowest-order transition then 
requires four photons. The further successive decreases in 
the rate as one proceeds through NaCl, CaF^, NaF, and LiF are less 
pronounced: 4>^,2,21,20, respectively. The trend of these 

numbers is reflected in the A data of Fig. 7. The much 
greater strength of the three-photon transition for KH 2 pO^, 
compared to the higher-order transitions for the other materials, 
is consistent with the significantly lower A value for that 
material from our experiment. The gradual increase in the 
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A value for > 7 eV is consistent not only with the gradual 
weakening of the higher-order transition probabilities, but 
also with the "low" frequency avalanche ionization. Joule 
heating model embodied in Eq. 1, Thus one sees that the data 
of Pig. 7 illustrates the continuous evolution of the character 
of the breakdown process as a function of material band gap 
from the "pure" avalanche ionization. Joule heating region of 
large into the regime where the influence of multiphoton 
transitions is observable. 

2. Self-focusing Behavior 

The experiments reported here provide qualitative 
verification of two aspects of self -focusing behavior in 
the tightly focused geometry, which to our knowledge were 
not investigated before. For the self-focusing experiments, 
we used the same experimental arrangement as in Fig. 3, 

However, three different lenses (focal lengths 1.0, 1.5, 

3.0 in.) were used rather than one as in the previous section. 
Pulses of increasing power were passed with each lens into 
the sample materials, one shot per site, and the threshold 
power for the production of a visible, light scattering 
plasma at the focal point was recorded. Note that this 
criterion for plasma scattering (PS) is not identical to 
the criterion for threshold breakdown as discussed in Sec. 

C.2. Great care was necessary to remove extraneous scattered 
light so that only the focal voliame was observed. Considering 
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the sensitivity of the eye at 5321 S to scattered photons, 
one would not be surprised that the PS criterion required 
less power to achieve than the breakdovm criterion of permanent 
microscopic damage. For the 1.0 in. focal length lens which 
was used for both experiments, the ratio of the PS threshold 
power to the breakdown threshold power varied with material 
between about 1.0 and 0.5. Nevertheless, the experiments 
showed that the PS criterion was sufficient to indicate 
reproducibly the attainment of a well-defined electric field 
in the sample. Because it is less time consuming to implement, 
the PS criterion enabled several materials to be examined. 

By plotting the PS power thresholds versus focal area for 

the three lenses, as discussed earlier and illustrated in 

Fig. 2, we investigated the trajectory of the breakdown curve 

in the dotted region of that figure. Referring to Fig. 2, 

consider the analysis of a data set containing points lying 

in the power range P^^ to P 2 . A straight line fit to those 

points will yield a vertical axis intercept below P^^ and, 

—1 

perhaps, even below ?2 , if the construction in the dotted 
region of Fig. 2 is correct. Our data for a typical material, 
Si02, is shown in Table III. The calculated focal areas, 
corrected for saturation as mentioned earlier, are listed in 
column 2, and PS power thresholds are listed in column 3. 
Columns 4 and 5 indicate that the power values almost span 
the P^ to ?2 region. If one plots this data as in Fig. 2, 
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TABLE IV. 3 

Self-focusing Data for Si02 


f 

(in. ) 

A 

(ym^) 

P 

(kW) 

P/Fl 

S/P2 

1.0 

15,9 

138.9 

0.49 

0.13 

1.5 

36.6 

260.2 

0.93 

0.25 

3.0 

155.7 

709.1 

2.52 

0.68 


^int 
^int ^ ^2 


= 1138 kW 
= 1.09 
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and fits a straight line to the three points, the intercept 

“1 

of the line v:ith the P a3:is yields a value of 1130 hW. 

As stated in Table III, this quantity indeed G::ceeds P 2 , 
verifying qualitatively the construction of the dashed breaJc- 
down curve as drawn in Pig. 2, A great need exists for accurate 
numerical computational results of <{){P) with which the data 
of Table III could be further tested. 

The second aspect of self-focusing behavior which has 

been qualitatively verified in our experiments concerns the 

exact input power value, in the tightly focused geometry, 

which will result in a singularity in axial intensity. In 

numerical computations, such a singularity is recognised as 

a rapid increase in intensity of orders of magnitude with 

only a small incremental increase in input power. In the 

experimental situation, the singularity is restricted to 

the breakdown intensity. As listed in Table III, in Si02 

we were able to focus up to 0.68 of P 2 before initiating 

a sizable plasma density. Of the six materials we tested, 

P/P 2 for the 3 in. lens was highest for 3102* The other 

materials registered P/P 2 between 0.32 and 0.57 for the 3 in. 

lens. All these P/P 2 values correspond to P/P 2 ^ values exceeding 

unity. Therefore, we have verified that in the tight externally 

28 

focused geometry, as discussed by Marburger, pov/er of 
is insufficient to produce a self-focus singularity. We v/ere 
not able with our lenses to verify that the full Pp is so 
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required, but we did determine that at least ■~0.7 of 
is required. This data could be extended further, of course, 
with longer focal length precision lenses. 

E. Discussion 

In Section D we interpreted the measured, frequency- 
dependent threshold changes (Pig. 7) in terms of the two- 
component picture of laser-induced breakdown which is 
schematised in Pig, 1, The interesting behavior there is 
provided by competition between the frequency trends of 
avalanche ionisation and multiphoton absorption. One would 
like to be able to directly compare the experimental data of 
Pig. 7 with existing theoretical predictions. At present, 
however, none have been presented for our experimental 
conditions of wavelength and pulse duration. This situation 
is due to the unavailability of experimental multiphoton 
absorption coefficients at 1.064 ym and 5321 S. This absence 
is unfortunate because the onJy available avalanche ionization 
rate data for laser-induced breakdown is for 1.064 ym wavelength 
in alkali-halides. 

Multiphoton ionisation coefficients have been measured 

(by photoconductivity detection) at 6943 S and 3472 S in NaCl, 

4 8 

KCl, and KI, however, by Catalano ^ and Dneprovskii 

49 21 

et Braunlich ^ have used that data, along with 

the 1.064 ym avalanche ionisation rate data (neglecting its 
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frequency dispersion) to investigate numerically the influence 
of multiphoton absoi'ption on breakdown in NaCl at 6943 S and 
3472 S (30 nsec FWHM pulse duration) . They calculated, as 
a function of light intensity, the conduction electron density 
generated by pure avalanche ionization, by avalanche ionization 
with five-photon interband absorption, and by avalanche with 
five-photon absorption plus an P-center population. Their 
results showed that multiphoton absorption, both interband and 
via P-center levels, is more efficient at producing initial 
electrons at low field strengths than is avalanche ionization. 
They also showed that the avalanche ionization efficiency 
dominates the plasma generation for electric field strengths 
over about 1 MV/cm, regardless of the presence of the other 

3 

multiphoton mechanisms. Our morphological data, at 1.064 pm 
and 5321 S (Pig. 5) , indicates that the last statement is 
indeed true. At both of those wavelengths we observed the 
same morphology, even in KH^PO^ at 5321 S where three-photon 
transitions are possible. At both wavelengths the morphology 
consisted of spatially distinct plasma sites separated irreg- 
ularly by short (-^1 urn) undamaged spaces. If strong multi- 
photon absorption were dominant in the breakdown process, one 
would expect to see, not spatially distinct plasma vestiges, 
but spatially continuous evidence of damage along the focal 
axis even at breakdown threshold. As Fig. 5 demonstrates, 
that is not the case. At 5321 S, in materials V7ith ^ 7 eV, 
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the character of intrinsic picosecond laser- induced break- 
down is dominated by the familiar avalanche ionisation 

21 19 

mechanism. This situation is not expected '' to persist 
for wavelengths much shorter than 5321 S. in our sample 
materials of band gap 7 to 12 eV. 

Another recent theoretical analysis has been made by 

A 25 

Sparks. ' He included electron-photon-phonon transitions 
as a new mechanism for sustaining an avalanche, and he also 
included multiphoton absorption for initial electron generation. 
His theoretical threshold values for 10.6 ym, 1.06 um, and 
6943 S laser-induced breakdown fall within the uncertainty 
limits of the experimental values. His calculations predicted 
a decrease in the NaCl threshold to be observed for wavelengths 
shorter than about 1.0 \xm {for nanosecond pulses). Our pico- 
second pulse data do not follow that prediction. They do not 
show the predicted ~505 decrease of the threshold between 
1.064 ym and 5321 S. If NaCl behaves in a way similar to 
KH 2 P 0 ^, our data suggest that a threshold decrease may be 
expected for wavelengths shorter than about 3hc/^^ 4300 S. 

Several other points raised in Refs. 4 and 25 merit discussion. 
It was stated in these papers that the difference between 

Q 

early DC thresholds and those measured by Yablonovitch was 
attributed to avalanche ionization frequency dependence, as 

9 

embodied in Eq. 1, That assertion is incorrect. Yablonovitch 
warned against such assignment of the systematic difference that 
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he found; and indicated that the more probable origin of 

. 19 

that factor concerns details of the DC e^perrments. Others 

have restated that point. Concerning the avalanche ionization 

frequency dependence as given by Eq. 1, both the present work 

11 

and an earlier experiment indicate that is of the order 
2 X lO" ” sec; whereas the calculations of Sparks predict 
that T_ should be larger roughly by the factor 10. A final 
point concerns the initial density of electrons that may 
serve to launch an avalanche. It is stated by Sparks that 
because of insulator electrical conductivity and photocon- 
ductivity values, initial electron densities should be less 
than -^lO^ cm"”^ and that densities as large as 10^^ cm ^ 
are unreasonable. Yet direct experimental evidence of such 
large densities has been found in picosecond laser- induced 
breakdown morphology. The resolution of this point lies 
in the fact that the photoconductivity measurements, to which 
Sparks alludes, must involve light intensities too weak to 
permit the ionization of color center traps by multiphoton 

absorption, and photon energy too small to allow single-photon 

11 -3 

ionization. Color center densities as low as 10 cm are 
difficult to detect by linear absorption spectrophotometry. 

It is therefore possible to have a clear NaCl sample which 
satisfies the insulator criterion (dc conductivity a ~ 

ohm""^ cm”^) , yet simultaneously possesses cm ^ 

F-centers which may be ionized by 1.064 ym tvjo-photon 
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absorption to produce breakdown morphology as observed in 
Ref . 3 . 
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CHAPTER V. 


STUDY OF PICOSECOND LASER- INDUCED BREAKDOWN AT 3547 S 

A. Introduction 

This chapter relates our effort to accurately measure 
breakdown thresholds in a series of materials using 3547 S 
picosecond laser pulses. In order to investigate the widest 
possible range of material band gap, the series of materials 
included as many as possible of the materials tested at 1.064 
urn and 5321 S (Chapters III and IV of this dissertation) : LiF, 

NaF, CaF^f SiO^, and KH 2 FO^ (KDP) . These materials cover a 
band gap range of about 12.0 to 6.9 eV and, with a 3547 2 photon 
energy of 3.50 eV, present the opportunity to study dielectric 
breakdown assisted by multiphoton transitions of various order, 
all the way down to two-photon absorption in KH 2 PO^. Results 
from this study have provided qualitative verification and 
extension of the documentation of frequency dependent breakdown 
behavior, as first described in Chapter IV, 

Section B of this chapter describes the experimental 
equipment, procedure, and difficulties. The results are pre- 
sented in Section C and discussed in Section D. 

B. Experiment 

1 . Apparatus 

The components of the experimental apparatus are 
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illustrated in Fig. V.l. The optical input to that series 
of components is the single, 1.064 pm, nominally 30 psec, 
laser pulse from the laser system described in Chapter II 
of this dissertation. The characterization and frequency 
conversion of this infrared pulse is begun with the components 
within the dashed line. The potassium dihydrogen phosphate 
(KDP) crystal HGl and vacuum photodiode detectors D1 and D2 
yield the two energy signals required for the pulse duration 
measurement by the ta technique, as described fully in Appendix 
II. The infrared pulse is next frequency-doubled in a second- 
harmonic generation cell HG2 containing a cesium dihydrogen 
arsenate (CDA) crystal. This crystal is electronically main- 
tained at a temperature of 42 ®C, which is the temperature 
that noncritically phase-matches the oriented CDA crystal. The 
benefits of noncritical phase-matching^ are the absence of 
deleterious walk-off effects due to double refraction, a 
relatively large acceptance angle of input beam divergence, 
and stability of Gaussian spatial profile for an unlimited 
period of time. These properties, not obtainable with the 
more familiar angle-tuned KDP doubling crystal, are essential 
for the production of precisely known, optical electric field 
distributions. The conversion efficiency for second-harmonic 
generation is monitored with the filtered detector d 3. The 
limit to the amount of 3547 2 light that we may produce with 
this equipment is set by the damage threshold of the CDA 
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Pig. V.l Diagram of the experimental apparatus. Notation: 
D1-D4, biplanar photodiode detectors; P, appropriate filters; 
HGl, KHjPO. frequency doubling crystal; HG2, CsH^AsO^ frequency 
doubling crystal; FMC, KH2p0^ crystal for 3547 generation; 
p, prism; LPA, Reticon linear photodiode array; C, curtain; 

L, Si02 lens; S, dielectric sample. 
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matorial. At the conservative maximum incident 1.064 ym 
irradiance level we have used :i 10 W/cm with 100 MW 
single pulse power) / the power conversion efficiency was '^SOvi, 

The 1.064 jim and 5321 8 pulses exit collinearly from 
the CDA cell into the final frequency mixing crystal (FMC) . 

In this crystal the orthogonally polarized infrared and second- 
harmonic pulses are mixed to produce light at the third-harmonic 
wavelength of 3547 8. This frequency tripling crystal is made 
of KDP and is phase-matched by angle-tuning. Unfortunately, 
no material is available with the proper temperature behavior 
of refractive indices (ordinary and extraordinary) to allow 
temperature- tuned, 90° phase-matched frequency summing at 
these two input wavelengths. The maximum power we have obtained 
at 3547 8, under the conservative limiting conditions mentioned 
above, is about 8 MW for a power conversion efficiency of 
from 1.064 ym to 3547 8. Conversion efficiencies for frequency 
tripling of up to 30^ have been obtained by others^ using longer 
crystals. Further detailed characteristics of the 3547 8 pulse 
will be related in the following paragraphs. Spatial separation 
of the three output pulses from the tripling crystal is performed 
by quartz prism P. The UV pulse is then sampled by two quartz 
beam splitters, before passing through a light-shielding curtain 
and into the damage chamber. The first of the two samples is 

4 

directed onto a commercial Reticon self-scanning linear array 
of 256 miniature photodiodes (LPA) . The diodes are spaced by 
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0.05 Him# and are incorporated onto an integrated circuit chip. 
This device, coupled to an oscilloscope, displays the spatial 
energy profile of the impinging light pulse. With it, we are 
able to measure the area of the 3547 S pulse and check its 
spatial profile quality. Further mention of its use will be 
made below. The second sample of the main pulse enters a 
calibrated photodiode D4, in order to provide a measure of the 
UV pulse energy just before it enters the lens. The charac- 
terised UV pulse finally is tightly focused by a quartz lens 
(L) into the dielectric sample (S) under investigation. The 
sample is held in a three-dimensional translation holder so 
that the pulse may be directed through any volume element of 
the sample. 

In Chapter II of this work, it was made clear that three 
pulse parameters must be carefully measured to allow the 
accurate determination of the electric field distribution of 
the optical pulse. The sufficiency of these three parameters - 
the spatial profile radius p, the temporal half-duration t, 
and the total pulse energy - is based on the assumption that 
the pulse spatial and temporal functional form is unchanging. 

We will now consider the degree to which the characteristics 
of our 3547 S pulses allow us to determine the field distribution 
in this experiment. 

The third-harmonic pulse energy measurement is 

straightforward, as it was in the previous 1.064 pm and 5321 S 
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experiments. The detector D4 in Fig. V.l was calibrated 

5 

again with tv 7 o Eppley Laboratory thermopiles, which were 
themselves calibrated with National Bureau of Standards spectral 
sources. Again, the detector area is large enough and the 
detector-oscilloscope response time is long enough that true 
spatial and temporal integration is performed over the pulse 
waveform by the detector. Our calibration method yields pulse 
energy measurements vjith an estimated accuracy of better than 
±105. 

The temporal waveform of the 3547 S pulse is closely 
Gaussian for small conversion efficiencies in both the doubling 
and mixing crystals. In that case, the average duration of the 
UV pulse will he l/i/3 of the 30 psec 1.064 ym parent pulse 
duration (or 17 psec) . This l/v/3 degree of pulse shortening 
has been roughly confirmed experimentally with streak camera 

3 

photography. As the conversion efficiency increases, pref- 
erential saturation of the doubling and mixing processes at 
the temporal peak region of the pulse results in an effective 
lengthening of the duration. The saturation results simul- 
taneously in an effective broadening of the pulse area as the 
pulse becomes more square in profile. Fortunately, our experi- 
mental procedure - specifically the tA technique (Appendix II) - 
allows us to incorporate these broadening factors into our data 
analysis. In the low conversion regime, it is easily shown 
that the ratio of the fundamental (f) and third- 

)F THi^ 
POOE 
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harmonic (TH) pulse energies is directly proportional to the 
third-harmonic pulse area A and half -duration t and is inde- 
pendent of the pulse intensity. For date in the low conversion 
regime, we may normalize the individual pulse^c^^ / signals 
to a product of a half-duration t of 17 psec and the actual 
pulse area as measured with the photodiode array. Then for 
the remaining data points in the higher power range, the 
signals will become progressively larger and will properly 
take into account the broadening of the duration and area- 

We have discussed the energy and pulse duration measure- 
ments and the minor saturation broadening of the pulse area, 

2 

irp . It remains to consider the distortions in the pulse 

spatial profile which are due to walk-off and self-focusing 

in the KDP mixing crystal. Pig. V.2 illustrates the typical 

spatial distribution of energy that we obtained in a 3547 8 

7 

pulse. Walk-off is the term which describes the decrease in 
overlap, proportional to the distance traveled inside the non- 
linear crystal, between the nonlinear polarisation wave and the 
generated light wave. Walk-off is caused by crystal birefrigenc 
and it results in a cylindrically asymmetric, distorted spatial 
waveform for the 3547 8 pulse. The poor spatial profile exhib- 
ited in Fig. V.2a is the result in the plane most severely 
affected. Fig. V.2b shows the distortion which occurs in the 
plane normal to that of Fig. V.2a. One observes in these 
photographs the departure of our spatial profile from a 
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(a)‘ 



(b) 

Fig. V.2 Oscillographs of the spatial 
energy distribution of the 3547 8 opti- 
cal pulse. Figures (a) and (b) are 
taken in the orthogonal planes of most 
and least influence from walk-off, resp. 
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Gaussian form. 

Nevertheless, the input area of the UV pulse of Fig. V.2 

must be estimated by fitting to it a Gaussian curve, as shown 

in Fig. V.3. The input profile radius p (0.81 mm) obtained 

from the dotted Gaussian curves of Fig. V.3 yields (with Eg. 3 

of Chapter III, Part A) a focal area A {=itw^) of 5.6 pm^. If 

the laser profiles of Pig. V.3 fit the dotted Gaussian curves 

perfectly, the uncertainty in A would be less than ±20^, as was 

the case in the experiments at 1.064 ym (see Pig. 2, Part A of 

Chapter III). The actual fit of the profile in Fig. V.2b is 

acceptable, but because of the modulation evident in Fig. V.2a, 

the fit there is poor and leads to additional uncertainty in 

the focal plane area. Since the lens essentially performs 

a Fourier transformation, the abrupt transverse intensity 

variations (high spatial-frequency modulation) of the input 

pulse result in a deconcentrafcion of the energy distribution 

in the focal plane which is not taken into account by the 

Gaussian beam propagation equations used above to arrive at 

2 

the area of 5.6 pm . Rough estimation indicates that the 
actual effective area might be as much as a factor two larger 
than the 5.6 ym value obtained by Gaussian approximation. 

If the laser spatial profile were constant from pulse 
to pulse over the period of time necessary for an experiment 
to be performed, the scanning aperture technique (Chapter III, 
Part A) could be used to measure the effective area in the 
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Ca) 



(b) 

Fig. V.3 Illustration of the approx- 
imate fit of a Gaussian curve to the 
laser pulse spatial profile. The pro- 
file radius p is used to compute the 
area of the pulse in the focal plane 
of a lens. 




focal plane. That avenue of relief was blocked in this 
experiment by the last distortions which must be mentioned 
and which are both time- and power-dependent. Over a period 
of a few (~3 to 5) minutes ^ minute movements in the optical 
table resulted in variations in the KDP phase-matching angle/ 
and the large-scale fringe structure in Fig. V.2 wandered 
correspondingly . 

At high power levels self-focusing of the 3547 S pulse 
in the KDP crystal and in the subsequent quartz elements sig- 
nificantly enhances the noise structure produced by the phase- 
matching conditions. Such self-focusing/ which grows from 
small-scale spatial hot spots on a pulse/ is known as "'small- 

O 

scale self-focusing" . This latter spatial instability is 
to be distinguished from the "whole-beam" or "large-scale" 
self-focusing which was discussed in the breakdown experiments 
of Chapters III and IV. Enhancement of the hot spots on the 
pulse by se_ f-focusing was easily observed by removing the 
sample (see Fig. V.l) and displaying the diverging 3547 R 
light on a fluorescent card (or bleached paper) . The strong 
power dependence of the enhancement was clearly observable. 

Hot spot enhancement by self-focusing was observed to occur 
for input power levels over roughly 150 kW/ and was present 
in the breakdown tests of all materials except CaF 2 . This 
effect results in different/ uncalculable reductions in the 
effective focal area for materials which undergo breakdown 
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at different power levels. The influence of this effect is 

in opposite direction to the enlargement of the focal area, 

mentioned earlier, caused oy the Fourier transform of high 

spatial-frequency components on the input pulse. Therefore 

the true effective focal area for CaF^ is probably nearer the 
2 2 

11,2 ym estimation than the 5.6 pm value, but for materials 

with increasingly higher threshold power, the probable effective 

2 

area is reduced to values closer to 5.6 pm and perhaps less 
by hot spot enhancement. This material dependent effective 
area makes difficult the extraction and substantiation of even 
relative breakdown thresholds. 

The conclusion of this section is that the best 3547 S 
pulse that could be produced was decidedly less suitable for 
absolute electric field measurements than the 1.064 pm and 
5321 S pulses used in previous experiments. Aside from any 
consideration of whole-beam self-focusing, the spatial degra- 
dations mentioned above impose an uncertainty on the focal 
point power per unit area of roughly 100^. Repercussions of 
this uncertainty will be discussed in Section C of this chapter. 
2 . Procedure 

The procedure for this experiment is identical to that 
of the 5321 S experiment as related in Chapter IV (pages iv- 
20 to IV-23) and the reader is referred to those pages. One 
of the six materials studied at 5321 S, NaCl, was eliminated 
from the 3547 S experiment because of distortion by surface 
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polishing scratches. Due to the magnitude of the breakdown 
input power relative to the self -focusing critical power, data 
in this experiment was talcen only with our shortest focal length 
(0.5 in.) lens. 

C. Results 

The data from this experiment are listed in Table V.l. 

The five tested materials are KH 2 PO 4 (KDP) , fused Si 02 / CaF 2 / 

NaP, and LiF. These five materials were studied at the two 
previous frequencies of 1.064 um and 5321 S. The linear re- 
fractive indices n at 3547 2 are listed in column 2. Columns 

o 

3 and 4 tabulate the "paraxial" critical power and the 

"catastrophic" critical power P 2 , calculated using the nonlinear 

refractive index (n 2 ) values listed in Table 1 of Chapter IV 

(page IV-26) • These quantities were introduced and discussed 

2 ? 

at length in Chapter IV. The quantity equals cl / 32 it n 2 , 

where c is the vacuum speed of light and X is the light wave- 

2 

length in air. ?2 is equal to 3.7 P^. Because of the X 
dependence of and P 2 / the 3547 2 critical power values in 
Table V.l are smaller than the corresponding entries in Table 1 
of Chapter IV by a factor of 9/4. In column 5 the measured values 
of the pulse input power at the breakdown threshold are listed, 
and columns 6 and 7 tabulate the fr ictionaJ quantities 
and P/P 2 * 

As derived in Chapter IV (see Eq, 6 ) the breakdown 
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3547 S Breakdown Data 


Material 

^0 

^1 

(kW) 

^2 

(kW) 

P 

(kW) 

P/Pl 

P/P2 

t{)(P) 

„uv 

^B 

(MV/cm) 


KH 2 PO 4 

1.4860 

119 

439 

242 

2.03 

0.53 

2.5“3.7 

35-77 

1.5-3. 3 

Si 02 

1.4762 

125 

462 

170 

1.36 

0.37 

2.0-3.0 

25-49 

1.3-2 . 6 

CaF 2 

1.4463 

105 

390 

73 

0.69 

0.19 

1. 0-1.5 

18-33 

0.7-1. 3 

NaF 

1.3359 

126 

465 

271 

2.15 

0.58 

2. 5-3. 7 

40-99 

2. 1-5.1 

LiF 

1.4024 

129 

477 

245 

1.90 

0.51 

2. 5-3.7 

35-70 

1.3—2 . 6 
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threshold electric field Eg may be obtained from the threshold 
power p with the relation 

-.1/2 


E 


B 


n^E„c A 
o o 


1 - 




( 1 ) 


All the quanuo-ties in Eq. 1 have been defined in this chapter 
except the pemittivity of free space, and <}>(P), the self- 
focusing function introduced in Chapter V. To recapitulate, 
the known limits of (Jj(p) are <J)(P) 1 for P << P^ and <|){P) 3,7 

for P ^ P 2 * Between these two limits ij)(P) is presently unknown. 

The obstacles to the interpretation and extraction of a 
breakdown threshold electric field from the data of Table 1 are 
two-fold. One of them has already been mentioned and it is 
the substantially non-Gaussian and unavoidably poor spatial 
character of the 3547 S pulse due to the details of the final 
frequency mixing process and self-focusing. The specific result 
of the propagation of that distortion into the focal plane is 
that the intensity distribution there is irregular with an 
effective area A known only to within a factor 2. The probable 
effective area varies with input power, as well. 

The other obstacle is the magnitude of the breakdown thres- 
hold power P relative to P^^ and P 2 * As the entries in columns 
6 and 7 show, the P values fit neither of the two limits mentioned 
above. The CaP 2 input power more closely fits the P << P^ limit, 
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while the other input powers more closely approach the P ~ Pg 

limit. It is therefore necessary to consider different ranges 

of the (j) (P) function for the different materials, in order to 

deduce realistic upper and lower bounds for Eg from Eg. 1. 

Column 8 of Table V.l lists the considered ranges of In 

comparison to the large uncertainties introduced by the area A 
2 

(5.6 to 11.2 ]jm ) and by t}i{P) in Eg. 1, the small (±10^5) uncer- 
tainty in the input power P is immaterial and therefore neglected. 

The next to the last colxmin of Table V.l lists the 3547 S 
breakdown threshold (E^^) results obtained in this experiment- 

xp 

U.V * 

The stated range in Eg results from the uncertainties in A 
and (j) (P) . The lower and upper limits of Eg^ for each material 
were obtained with the use of the upper and lower, respectively, 
limits of A and ({)(P). The range of Eg^ for a given material in 
column 9 constitutes an uncertainty of ~+35& about the average 
of the upper and lower limits. The reason the data were not 
displayed in that manner is that the power- (and thus material-) 
dependent hot spot behavior points to different probable location 
of the true threshold within that ±35& range for the different 
materials. For example, for CaF^ the true threshold is closer 
to 18 MV/cm than 33 MV/cm, due to the small factor and the 

resulting absence of hot spot enhancement. But for LiF, the 
true Eg^ may be nearer 70 MV/cm than 35 MV/cm, 

In order to find whether interesting frequency dependent 
behavior is resolvable in this data, in Fig. V.4 we plot the 
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ratio / Ef (or rather the range thereof) versus material 
band gap. The 5321 2 thresholds (Eg) are taken from Table II 
of Chapter IV. Here we ignore the common effect of the small 
pulse duration difference (nominally 17 psec at 3547 2 versus 
21 psec at 5321 2 ) on the threshold ratio. 

The relative and absolute positions of the data bars 
for CaF^ and LiP seem reasonable, in view of the similar data 
of Chapter IV (Fig. 7) . We also observe an absolute increase 
in the NaF threshold, as expected. Yet the reason for relative 
magnitudes of the WaF and LiP thresholds is not known, unless 
it is the simple manifestation of the large experimental uncer- 
tainties. The observation that the Eg^ / Eg values for KH2P0^ 
and Si02 in Fig. V.4 are as large or larger than that for 
LiP can be explained, however. In Fig. 7 of Chapter IV, the 
interpretation was advanced that the marginal onset of three- 
photon absorption (3PA) lowered the KH2P0^ threshold at 5321 2 
relative to that at 1,064 ym. If we apply that observation 
to the present 3547 2 experiment, we conclude that the SPA 
might serve to lower the thresholds for materials with band 
gap less than 10,5 eV. CaF 2 has a band gap of 10 eV and indeed 
in Fig, V.4 one observes that the threshold ratio bar for Cap 2 
allows the possibility that the 3547 2 threshold is less than 
the 5321 2 threshold. For Si02/ 3PA would be considerably 
stronger due to its band gap value of about 8.0 eV. It is 
expected that strong 3PA caused beam depletion along the focal 



KH2P0^ Si02 


Cap 2 NaF LiP 


Fig. V.4 Illustration of the band gap 
dependent behavior of the threshold ratlo^ 
observed in the present experiment, and 

Eg are the 3547 § and 5321 S thresholds, resp 
as tabulated in Table V.l. The origins of 
the uncertainty in the data are discussed in 
the text. 
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axis in the SiO^ sample and thus caused the value 

exceeding unity in Pig, V.4. A similar explanation may be 
advanced for the even larger / E„ values observed for 
KH2P0^. For that material/ two -photon absorption is energet- 
ii.ally allowed at 3547 S, and beam depletion before the focal 
point would be more severe, 

D. Discussion 

We see that the data of Pig. V.4 is in large both explicable 
and reasonable in light of the earlier data obtained at the wave- 
lengths of 1.064 ym and 5321 S and displayed in Fig. 7 of Chapter 
IV. Although the uncertainties introduced by pulse distortion 
did not allow the determination of precise breakdown thresholds 
at 3547 S/ the frequency-dependent behavior observed first at 
5321 S was verified and extended qualitatively. Larger thres- 
holds were found at 3547 2., than at 5321 2 or 1.064 lim/ for 
materials which do not permit two- or three-photon interband 
absorption to occur at 3547 2 . For one material/ CaF 2 / which 
supports weak three-photon absorption at 3547 2 , a decrease in 
the 3547 2 threshold/ compared to its 5321 2 threshold, was 
found. This threshold behavior at 3547 2 in Cap^ corresponds 
energetically to that found in KH 2 pO^ as the light frequency 
was doubled to 5321 2 in the prior investigation. At 3547 2 / 
for materials (KH 2 PO^ and Si 02 ) with optical band gap smaller 
than that of CaF 2 , an increase in the threshold input power 
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was also observed. However, those increases were interpret- 
able as arising from depletion of the laser pulse, by strong 
two- or three-photon absorption, as it passed through the 
sample toward the focal region. 

The indications of strong multiphoton absorption in 
KH^PO^ and Si 02 suggest that if the breakdown experiment had 
been conducted on the sample surface rather than in the bulk, 
large decreases in the thresholds might have been observed at 
3547 S, The study of breakdown on surfaces would alleviate 
two difficulties which complicated the study described in 
this chapter - self-focusing in the sample and multiphoton 
absorption pulse depletion in the sample. Because of the 
wavelength dependence of these two factors, surface break- 
down study will be increasingly attractive for wavelengths 
below 3547 S. The use of 2661 S fourth-harmonic light, 
generated entirely with 90° phase-matched materials (CDA and 
ADP) , for such a study would remove walk-off distortion, as 
well. The benefits of surface breakdown study in the ultra- 
violet spectral region will be offset somewhat by other well- 
known problems associated with surface breakdown. Those 

9 

problems originate with surface impurities and imperfections. 
Because they generally introduce threshold electric field 
uncertainties of - \t a factor 2, the surface quality problems 
will require serious attention and solution if accurate and 
reproducible thresholds are to be obtained. 


i 


V-21 

The main conclusion of this chapter is that the break- 
down threshold does increase with light frequency , as predicted 
by the avalanche ionization model, up to the frequency where 
three-photon absorption may take place. Beyond that frequency, 
the assistance of the essentially instantaneous multiphoton 
absorption transitions dominates the avalanche ionization fre- 
quency trend, with the result that the breakdown threshold is 
reduced. It can now be said that the experimental study of 

laser-induced breakdown has been extended to the quantiam regime, 

10 

and that the predicted behavior of the breakdown threshold 


has been observed. 
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pulse duration measurements* 
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Analytical and numerical results are obtained for two-photon fluorescence in an attenuating medium. Coses 
of one- and/or two-photon attenuation arc considered. These results are interpreted with respect to the 
measurement of ultrashort laser pulses. 

FACS numbers: 42.65.D, 42.60.N 


The method of two-photon fluorescence (TPF) has 
been widely used to measure the duration of ultrashort 
light pulses produced by mode-locked lasers.^ In gen- 
eral, if two oppositely directed pulses of equal duration 
overlap each other in a two -photon absorbing medium, 
then the spatial distribution of the fluorescence decay 
may be used to infer the duration of the pulse. A typical 
experimental arrangement for measuring the pulse 
length is exhibited in Fig, 1. 

For simplicity we will assume that the laser pulse 
has a uniform spatial profile and that there is total 
spatial overlap of the two oppositely directed pulses 
when they intersect each other in the center of the cell. 

If the total time integrated fluorescence at a position z 
is denoted by S{z), then 

S(z)cc r ^l{t-z/v)dt+ + s/v)dt 

+ ij Iid- z/v)l 2 {t+ z/v)dt, (1) 

where Ii and are the irradiances of the two pulses. 

The position 2=0 corresponds to the position of com- 
plete overlap of the two oppositely directed pulses, and 
V is the group velocity of the pulse in the fluorescent 
medium. 

F or the case of equal intensities and no attenuation of 
the pulses due to either scattering or one- or two-photon 
absorption, Eq. (1) reduces to 

S(2)oc 1 + 2G‘®'(22 /w)[G<">(0)]-‘, (2) 

where 


M 



FIG. 1. Schematic diagram of the TPF experiment, M, mir- 
ror; BS, beam splitter. 


Gf^Hx)=^S”iit-xmdt. 

mCO 

This result has been obtained previously^ and implies a 
contrast ratio /f=S(0)/S(<«) = 3, We wish to note, how- 
ever, that contrast ratios less than 3 will be obtained if 
the pulses have unequal irradiances. Moreover, the 
spatial distribution of the TPF decay will also be dif- 
ferent if there is pulse attenuation in the TPF cell. 

For the case of TPF in a cell of length 2d and with 
significant beam attenuation via one-photon absorption 
or scattering, Eq. (1) becomes 

5 ( 2 ) a: (ia)^^®Tj^o{exp[- 2fli(d + a)] + 6^ exp[- 2a{d - z)] 

+ i6exp[-Zad-2^/(r^v^)]}. (3) 

Here ^10 is the irradiance at 2 = -d and Is the ir- ' 
radiance at 2 =d. The one-photon attenuation coefficient 
is O', T is the 1/e pulse half -width of the assumed tem- 
porally Gaussian pulse, and & is equal to 

If /io=f 2 o= All the result becomes 

S(2)cc2(-^ir)^/®ri5exp(- 2 a d) 

x{cosh(2a2) + 2 exp[~22V(T“w“)]}‘. (4) 

The results of Eq. (3) are presented In Fig, 2 for 
several values of 8 and a. 



FIG. 2. Plot of two-photon fluorescence intensity S ftn arbi- 
trary units) vorsuQ 2. The total TPF cell length Is 10,0 cm. 
The pulse duration full width at half-mmtlmum le 30 psec. 
Curve &) 0 = 0, 0=1.0; curve (2) a = 0, 0 = 0. 75; curve (3) a 
= 0.13 cm"', 0=1.0; and curve (4) a = 0,13cm’', 0 = 0.75, 
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FIG. 3. Plot of hvo-photon fluorescenco Intensity S (In arbi- 
trary units) versus z. The total TPF cell length is 10, 0 cm, 
and the maximum trradtance Is 3 GW/era?. The pulse duration 
full width at half-maximum la 30 psec. Curve ft) a =0, 0 
= 1.0, fl = 0; curve (2) o: = 0.13 cm-*, 0=1.0, d=0; and curve 
(3) a = 0.13cm*‘, 0 = 1.0, d = 8.4xio-‘® cmM. 



FIG. 4. Daneltomoter trace of two-photon fluorescence pro- 
duced optical density versus s. This trace results from a sln- 
i;le laser pulse extracted from a Nd : YAG mode-locked laser. 
The TPF medium is a lO"® solution of Rhodamtne GO In ethanol 
placed In a 10-cm-long dye cell. 


If both one -photon attenuation and two -photon attenua- 
tion are significant, the results become somewhat more 
Complex. For the case of a spatially uniform but tem- 
porally Gaussian pulse the irradiances are given by the 
expressions 


Itit ~ z/v) = /jo exp ^ expt- a{d+ s)] 


x{l-exp[— aW + x 




and 


Lit + z/v) = Isa eJT exp[- a{d - z)] 

x(l -exp[- 0 !{d- 2 )]yj . 


Here /3 is the two-photon attenuation coefficient and is 
related to the two -photon absorption cross section 
the number of two-photon absorbing molecules per unit 
volume B, and the photon energy Jiw by the expression 


If the two-photon absorption cross section has units 
cm*®, and tiu is expressed in Joules, then fi Is expressed 
in cm/W. 

No simple analytic expression may be obtained from 
Eq. (1) in this case, but numerical results are presented 
in Fig. 3 for several examples. The numerical exam- 
ple for and j3?a 0 corresponds to a 10*® M solution 
of Hhodamine 6G in ethanol for a laser wavelength of 


1.06 pm. The incident laser pulse is assumed to be 
temporally Gaussian with a pulse duration of 30 psec 
full width at half-maximum; the maximum irradiance is 
3 GW/ cm®. For this example a= 0, 13 cm*^ and 13= 8. 4 
xl0*>® cm/W." 

From Figs. 2 and 3 we observe that if there is 
significant attenuation of the pulses in passing through 
the TPF cell by one-photon absorption and/or two- 
photon absorption, then the spatial dependence of the 
TPF profile may change significantly from the results 
predicted assuming no depletion of the oppositely 
directed laser pulses. It is evident from these figures 
that one must use care in interpreting contrast ratios 
derived from these figures. 

Finally, we note that we have performed several 
experiments using single pulses produced by a mode- 
locked Nd :YAG laser to observe TPF, These experi- 
ments clearly show that there is significant attenuation 
of the laser pulse as it passes through the dye cell. A 
typical densitometer trace is exhibited in Fig. 4. The 
asymmetry of the trace is produced by the inequality of 
the intensities of the two oppositely directed pulses in 
this experiment. Our experimental results are in rea- 
sonable quantitative agreement with the theoretical de- 
scription presented here. 
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Wc describe here a simple nonilnear optic technique for the measurement of the duration of individual 
picosecond pulses. The accuracy and relative simplicity of the technique increase with the number of pulses 
measured. An experimental test of the basis of the technique is described. 

PACS numbers: 42.60.N, 42.65. 


1. INTRODUCTION 

Altliough a decade has passed since tiie generation 
of the first mode-locked laser pulses, the measur - 
ment of the duration of individual picosecond laser 
pulses is still accomplished only with difficuUw. The 
only direct method of measurement is streaJ: tmera 

photography. ‘ The techn' gy of streak photogj aphy has 
advanced to the point that subpicosecond resolution has 
been achieved, ^ and clearly strealc photography is the 
optimal method of measurement when the experiment 
involves a limited number of pulses. However, the cost 
of a streal: camera may be prohibitive. In addition, 
the associated photographic development and photoden- 
sitometric analysis of the data is time consuming. In 
tlaat respect, a streak camera is not ideal as a simple 
and rapid monitor when many pulses are involved. 

Single-pulse two-photon fluorescence^ (TPF) is the 
most widely used method of measuring temporal dura- 
tion, The advantages of this technique are experimental 
simplicity and economy. Several factors may hamper 
the straightforward Interpretation of the data, however, 
as attested by the literature concerning the attainability 
of a contrast ratio of three.®’'* Another such factor, 
that of nonlinear film response, will be mentioned here 
below. The actual temporal profile of the pulse cannot 
be reconstructed from the TPF trace. Indeed, one must 
assume a specific temporal shape in order to extract 
the pulse duration. The accuracy of that duration, there- 
fore, depends on the degree to which the actual tem- 
poral profile conforms to the assumed profile. In tliis 
work with a Nd: YAG mode- locked laser we have as- 
sumed a Gaussian temporal profile. This assumption 
has been verified by others. ** The TPF technique also 
involves photographic development and photodensitome- 
tric analysis for accurate measurements. Therefore, 
TPF measurement is also not ideal for use when many 
pulses are involved. It should be mentioned, however, 
that the use of a self-scanning photodiode array as the 
detector in a TPF arrangement® eliminates some of the 
above disadvantages of the TPF method. 

We will describe in Sec. n a simple and economical 
technique which is a useful alternative to the above 
methods for many experiments, especially when large 
numbers of pulses are used. It is capable of measuring 
individual pulse durations relative to the average pulse 
duration of the laser output. The proper characteristic 
average pulse duration must be determined accurately 
by other means, ® The technique presented should be 


useful for measuring single-pulse durations as short 
as 10'*® sec or less, provided the average pulse dura- 
tion of such a short-pulse generator can be measured 
by other means. Section in describes an experimental 
test of the theory of the technique, the analysis and re- 
sults of which are presented in Sec, IV. 


M. THEORY OF THE tA TECHNIQUE 


We begin the description of this technique by consid- 
ering a laser pulse having Gaussian spatial and temporal 
distribution of intensity: 

Hr, t) = /(, exp[- (r/p)®] exp[- (1) 

The energy content in a sample of this pulse will be 
proportional to the quantity given by 

- LZ /o°* Hr, i) 27rr drdt = ir* (2) 


where A = ffp® is the pulse area and t is the 1/c half- 
width of the intensity distribution. Such an energy mea- 
surement is performed by the beam splitter BS-1, 
photodiode PD-1, and the oscilloscope in Fig, 1. If a. 
second sample of such a pulse is converted to second- 
harmonic light, as by beam splitter BS-2 and the phase- 
matched second-harmonic generation crystal SHG-1 
in Fig, 1, then for small conversion efficiencies the 
enei'gy content of the pulse will be proportional to the 
qujjitity cfah given by 


*|™®VeiBILirY OWHli 



FIG, 1. Apparatus for pulse duration measurements with the 
tA technique. Notation: BS, wedged beam splitters; F, 
appropriate filter; SHG, phase-matched second-harmonic 
generation crystal j PD, btplanar photodiode; CRO, dual- 
beam oscilloscope. 
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The quantity is a collection of constants, for a given 
experiment, involvii^ the second- order nonlinear 
susceptibility, length of the harmonic-generation crys- 
tal, etc. In Eq, (3), the quantity A' = irp'* is the har- 
monic-pulse transverse area and t' is the harmonic 
pulse 1/e half-width of tlie intensity. For small con- 
version efficiencies, we know that and 

p = 2*^*p'. Therefore, the ratio can be written 

e}/u=Q'TA, ( 4 ) 

where q* is a second proportionality constant equal to 

Tills ratio is observed to be linearly propor- 
tional to the product rA, a fact that was noted by Glenn 
and Brienza’ in 1967 and employed by Shank and Ippen® 
in conjunction with a mode-locked cw dye laser. The 
ratio value is obtained simply from two oscilloscope 
signals or, alternatively, it could be obtained from 
digital equipment as a single numerical readout. In an 
experiment, therefore, if the value of the ratio 
is recorded for each individual pulse, one needs only 
to normalize this set of numbers in order to obtain in- 
dividual pulse values of tA, The normalization is easily 
acliieved if two conditions are satisfied: (i) The average 
value of the pulse duration (t) and the pulse area (A) 
must be known. These are quantities that are routinely 
measured in the characterization of a laser source, 

(ii) Enough shots must be monitored such that the aver- 
age of their recorded ratio numbers from Eq, (4) be- 
comes meaningful and representative. During tliese 
shots, Uie conations determining q‘ (the phase- match- 
ing angle of the harmonic generator, etc. ) must not 
be allowed to change. If these two conditions are satis- 
fied, then Hie average of the recorded ratio numbers, 
Identified as <tA), may be set equal to the product of 
the known, previously measured average values, 

The individual pulse ratio numbers are thereby 
calibrated absolutely in terms of this average. The 
quantity tA is just that needed to convert the easily 
measured individual pulse energy into individual pulse 
intensity. The latter quantity fluctuates from shot to 
shot in all mode- locked lasers, and heretofore has been 
difficult to measure. Usually, the average value (r) and 
(A) are used for all pulses, neglecting entirely their 
fluctuant behavior. If, to a reasonable approximation, 
the pulse-to-pulse fluctuatlonB in A can be neglected, 
or If a self- scanning photodiode array is used to accu- 
rately determine A for each pulse, then individual pulse 
durations are easily obtained. It should be noted that the 
accuracy of the technique, in addition to its relative 
simplicity over the alternaUve methods, Increases with 
the number of pulses utilized. 

One other point should be noted. It concerns the ef- 
fect of variation in the pulse temporal profile, and not 
just in the characteristic duration of that form, on the 
accuracy of the technique. It is not required that the 
laser pulses have Gaussian temporal and/or spatial 
functional form. A temporal and/or spatial form other 
than Gaussian will simply result in a different constant 
value of q' in Eq. (4). The values of <5//^'st, will still 
be proportional to the characteristic duration and area, 
So long as the proper characteristic average duration 
and area of the pulse have been determined by other 
means,® and are used in the normalization procedure. 


the technique presented here wUl function. However, 
if the functional form itself, and not just its charac- 
teristic area and/or duration, changes from shot to 
shot, then q' in Eq, (4) will no longer be constant and 
the strict proportionality there is lost. The accuracy 
of the technique will therefore be adversely affected. 

An example in which this limitation would be serious 
is that of a Nd : glass mode- locked laser without suf- 
ficient taandmdth reduction to remove noiselike 
substructure. 

III. EXPERIMENT 

We describe next an experiment that was conducted 
to test the validity of the linear correlation of <f//^"8n 
with tA as given by Eq. (4). The expe.'imental arrange- 
ment is schematically illustrated in Fig, 2, The 
Nd: YAG oscillator is mode locked by a 0. 5- mm- thick 
sheet of flowing Kodak 9860 dye. Transverse oscillation 
is constrained by an intracavity aperture to the TEMp# 
mode. The single-pulse switch-out mechanism is a 
standard Pockels cell pulsed by a laser- triggered 
spark gap. The components within the dashed line in 
Fig, 2 are tlmse of Fig. 1 and they were used to record 
ratio values of Eq. (4) for each pulse. Photodiode PD-3 
is a high-speed biplanar device used with Tektronix 
519 oscilloscope CRO-2 to monitor each pulse and in- 
sure that only data from single pulses without satellite 
structure is employed. Beam splitter 4 directs a frac- 
tion of the pulse into a phase-matched second-harmonic- 
generation crystal SHG-2. The harmonic output strikes 
a MgO diffuser D and is photographed by camera C-1, 
This multilens camera® (MLC) is used to measure the 
harmonic pulse area A', from which the fundamental 
pulse area A is calculated. Kodak Royal-X pan film was 
found to be sensitive enough to yield good exposure 
from a single pulse. The remainder of the pulse is di- 



FIG. 2, Schomatio diagram of the experimental arrangement. 
Notation; MLO, Nd : YAG mode- locked oa dilator j SPS, olnglc- 
pulse selector emi loying a laser -triggered spark gap and 
Pockels cell; AMP, NdiYAG amplifier; BS, beam splitter; 

F, appropriate filter; SUG, phase-matched 10)P second- 
harmonld generation crystal ; PD, biplanar photodiode; CHO-l, 
Tektronix 555 dual-beam oscilloscope; CRO-2, TeWronlx 519 
osclllogcopo; D, MgO diffuaer block; C-1, multltens camera; 
C-3, camera; TPF, two-photon fluoreecence coll coatatnUig 
rhodamlne GG in ethanol. The componenta within the dashed 
line are those of Fig. 1. 
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but not in Uie TPF technique. Indeed, high-speed films 
used to photograph single-pulse TPF traces offer cor- 
respondingly lower latitude. 



TRACE TRACE 

(a) 



t[og*gj,)-Ol3 Iag*tih 


(bl 

FIG. 3, (a) Schematic Illustration of the TPF photodenaltom- 
eter output from one alioti (b) Hurter-Drlffield presentatloQ 
of the data from (a). The letter d denotes optical density and 
log denotes the common logarithm. The subscripts h, a, 
and h denote the unattenuated trace, attenuated trace, peali 
level, and background level, respectively. 


reeled into a TPF arrangement The TPF medium is 
rhodamine 6 G, and Kodak Royal- X Pan Film was used 
in camera C-2 to record the TPF traces. The TPF cell 
housing employed two mirrors to direct ttvo images of 
the TPF trace, one from either side of the cell, into 
the camera. However, one of these images was attenu- 
ated by a neutral density filter placed in its path inside 
the housing. The attenuated trace was used to check 
linearity of film response as will be described below. 

The experiment was performed by firing a series of 
single pulses and recording various signals. The two 
traces from CRO-1 yield values of The MLC 

data yields another independent measure of the pulse 
area A, and tlie TPF data yields an Independent value 
of T, for each pulse. Therefore, the product of A and t, 
from the MLC and TPF data, should be linearly corre- 
lated with tile value of frotu CRO-1, according 

to Eq. (4). 

IV. ANALYSIS AND RESULTS 

In order to determine with what accuracy the ratio 
detects variations in t and A, the MLC and TPF 
data must be accurately analyzed. The problem of film 
response nonlinearity Is absent in the MLC technique 


The problem of film response linearity was taken 
into account in our experiment by recording two images 
of each TPF trace, one of which was attenuated with a 
neutral density filter of known optical density. An ex- 
ample of the photodensitometer trace output from one 
shot is shown schematically In Fig, 3(a). The exposure 
e from the TPF has tlie spatial dependence 

6 ( 2 ) = foil + Q exp{“ 2{sw/ ?t)*]}, (5) 

assuming the incident pulse shape of Eq. (1). Here tlie 
z axis is tlie axis of propagation in the TPF cell, having 
its origin at the point where the two pulses fully over- 
lap. In Eq, (5), Q is equal to the contrast ratio minus 
1. 0, n is the refractive index of the TPF medium, and 
c is the vacuum speed of light. Figure 3(b) displays 
the Hurter-Driffield information obtained from Fig, 

3(a), The exposure values of the two peaks of Fig. 3(a) 
are related by the known filter attenuation factor, in 
our case 50*??. Thus, if log^g is the common logarithm 
of the exposure of the point corresponding to the density 
then the abscissa of the point corresponding to 
is log(egp) - 0. 3. The abscissae log(«gp) and log{£u^) - 0.3 
result from tlie background densities of the unattenuated 
and attenuated traces, respectively, in identical fashion. 
Therefore, only a relation between log(£^p) and log(£„s) 
is needed to plot all four points as shown in Fig, 3(b), 

In our analysis, use was made of the fact that 
was very small, i, e. , 0 , l(d„j, - dpb) or less, on all the 
traces from the experiment. This was due to fortuitous 
exposure conditions. Therefore, tlie film response was 
talten as linear over this small region, a very good 
approximation An effective slope y was computed as 
the average of the slope from the pair of background 
points and from the pair of points from the pealis of tlie 
curves. The points logt^ub) and log(e„p) - 0, 3 were then 
connected via the relation 

rfub = rfap + VLO . 3 + 10g(£gj,/£up) 1. (6) 

This process established tlie positions of tlie four points 
as shown in Fig. 3(b), A quadratic curve was then least- 
squares fit to the four points, yielding the solid curve 
in Fig. 3(b), The proper optical density dj/j that cor- 
rectly corresponds to tiie haU-intensity point of the 
fluorescence above the background was tlien found from 
tlie relation logei /2 = iog(egp + €gb) - 0,3 and the curve oi 
Fig, 3(b). The full v/idth IF of the TPF densitometer 
trace at the di/j level tiien yields tlie temporal 1 /c in- 
tensity duration via t= IKH/c( 21 n 2 )^'^®, where c is speed 
of light and n is the refractive Index of the TPF medium. 


TABLE I. Teat of the linear correlation with tA. 


Pulse 

No. 

ipacc) 


de.-, 

V, 

t'b 

1 

15.1 

1.00 

l.OD 

1.00 

* a.i 

It 

15.4 

1.00 

1.03 

0.&7 

0,0 

.1 

Ki.u 

1.14 

1.20 

o.ni 

-0.2 

4 

IH.I 

1.30 

1.51 

0,00 

-7.2 

5 

20.7 

1.H2 

1.71 

1.05 

- o,;i 


*thc iuU-wldtli— at— lialf- intensity duration ii> erjuaS to 
^No' raalizcd to J. 00. 


1067 


J. Appl. Phvn., Vot. 47. No. 3, March 1976 


EEPRODUCIBILrrY OF TiF W.L. smith and J.H. Scchtel 
ORICSNAL PAGE IS POOH 


1067 



1 


AH- 5 

In the above manner, “was deduced from TPF 
for each of the five laser pulses and the values are 
listed in Table I. The area Aui.q was found from the 
technique of Hef. 9 for each pulse. In the third column 
of Table I, the normalized values of the product of T^pp 
andytiiLc are listed. The normalized values of the ratio 
(5//(5a„ are listed in column 4, The ratio of the values 
of columns 3 and 4 is listed in column 5. The degree to 
which the fifth column ratio elements are mutually 
equal indicates the accuracy of tlie tA technique em- 
bodied by Eq, (4). The sixth column lists the variation 
of Uie entries in column 5 about their average. There- 
fore, the correlation of the quantity 6j/6a\i. 
product of tlie pulse duration t and area A is demon- 
strated to be linear within about ± 10%, 

The advantages of the tA technique outlined above 
are obvious: (i) minimal expenditure of pulse energy 
in the rA measurement itself (as opposed to 1. OG-jim 
single-pulse TPP),‘ (ii) freedom from error introduced 
by film nonlinearity, a problem often encountered when 
photographic film is used with picosecond pulses, e. g. , 
in TPF photography; (lii) low cost and simplicity of ap- 
paratus; (iv) rapidi^ of data acquisition and analysis, 
and ease of incorporation of digital methods. These ad- 
vantages suggest the attractiveness ot the technique as 
a permanent on-line monitor. The simple scheme de- 
scribed above is useful in obtaining accurate data for 
processes which are highly nonUneax in the intensity, 


such as laser-induced breakdown*® and electron emis- 
sion,** multiphoton ionization, etc. 
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PULSE WIDTH FLUCTUATIONS IN MODE-LOCKED Nd:YAG LASERS * 

J.H. BECHTEL and W.L SMITH 
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Received IS August 1975 

Using a nonlinear optical technique we have measured pulse width fluctuations in a passively mode-locked Nd; YAG 
laser. 


An important aspect of mode-locked laser pulses is 
the pulse duration. We have measured the pulse with 
fluctuations by a technique that is based on observing 
both the total pulse energy and the pulse eneigy 
created by second harmonic generation in a nonlinear 
crystal [1-3], 

Consider a laser pulse in the form of a spatially and 
temporally Gaussian pulse. In this case the irradiance 
is given by I{r, t) - /,„ exp(-(>-/ti)')exp(-(f/T)2). The 
total pulse energy Ef is given by the expression 

Ef = 

t m 

where A ~ nd^. For small second harmonic conversion 
efficiencies the total pulse energy created by 
phased matched second harmonic generation is given 
by 

where ^ is a constant that depends on the nonlinear 
susceptibility the length of the harmonic generation 
crystal, and other parameters that do not change dur- 
ing the experiment. Here t = 7/2*/^ and .4' = A/2, 
Thus the ratio becomes 

If we assume that both the shape of the pulse and 
the pulse area do not cliange from shot-to-shot tlien 
our measurement of /Tf/E'si, is simply proportional to 
the pulse duration. 

A schematic diagram of the experiment is exhibited 
in fig. 1 . The laser used for this experiment was a 

* Supported in part by the Joint Services Electronics Program 
and the National Aeronautics and Space Administration. 



Fig. 1. Sehemutic diagram of the experiment. Nd;YAG MLO, 
neodymium doped yittrium aluminum gartict modc-lotkcd 
oscillator; PC, Pockcls cell; LTSG, laser triggered spark gap; 
AMP, neodymium doped yittrium aluminum garnet laser am- 
plifiers: PD, photodiode; CRO, Tektronix 519 oscilloscope; 

F, filter; KDP, potassium dihydrogen phosphate second har- 
monic generation crystal; INT, integrator; EA, electronic 
amplifier; ADC, analogue to digital converter; SC, scaler; 

Tl’, teletype printer. 

flashlamp pumped Nd:YAG laser rod in the Brewster- 
Brewster configuration. The laser was passively mode- 
locked by a 0.5 mm thick flowing dye (Kodak 9860) 
in contact with the high reflectivity rear mirror. The 
output was restricted to the fundamental i. ansversc 
mode, TEMoo,by a 2 mm diameter aperture within 
the oscillator cavity. A single pulse was selected from 
the train of mode-locked pulses by a cylindrical-ring 
electrode Pockels cell. This pulse was subsequently 
amplified by two Nd:YAG amplifiers. 

A photodiode, PDl , monitored the laser pulse with 
a Tektronix 519 oscilloscope to insure that only a 
single pulse was selected from the train of pulses. Cali- 
brated photodiodes PD3 and PD2 measured respective- 
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PULSE DURATtON (ptsc) 


Fig. 2. Histogram showing pulse width variation for passively 
mode-locked Nd:YAG laser. The pulse duration is the full 
tvidth at half maximum (FWHM). 

ly the total pulse energy at the fundamental frequency 
and the total pulse energy at the second harmonic fre- 
quency, These photodiode signals were monitored by 
a digital data acquisition system and displayed on a 
teletype printer. The average pulse duration full width 
at half maximum (FWHM) corresponding to the meas- 
ured average value of was set equal to 30 psec. 
Tills value was determined for our laser by a second 
harmonic autocorrelation technique employing a 
Michelson interferometer [4] . Fig. 2 exhibits a histo- 
gram of the total number of pulses observed with 
pulse duration (FWHM) between the endpoints of the 
horizontal axis, e.g., 88 observed pulses between 29 
and 31 psec (FWHM). 

The statistical nature of passive mode-locking is 
well known [Sj, The variance of the fluctuations ob- 


served here is larger tlian that predicted by others on 
that basis. From ref, [5] one obtains 
where m is the number of locked longitudinal modes. 
For m = 100, StIt = 0.1 . A larger variation would be 
expected, however, if a more realistic model of the 
saturable absorber [6] were employed. In addition, 
natural nonuniformities in the flowing dye will in- 
crease the variation in the pulse width. 

In conclusion we note that substantial pulse dura- 
tion fluctuations can occur with passive mode-locking 
in Nd:YAG lasers. Similar pulse duration fluctuations 
have been observed in the initial pulses from the pulse 
train of passively mode-locked Nd:glass lasers [7], 
These fluctuations were, however, observed with a 
streak camera. 

We wish to acknowledge the support and encour- 
agement of Professor N. Bloembergen. 
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FOUR-PHOTON PHOTOEMISSION FROM TUNGSTEN ''' 
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Four-photon pliotoemission is observed from tungsten using carefully controUed and calibrated single pulses produced 
by a passivel,- mode-locked Nd:YAG laser system. The results arc shown to be In quantitative agreement with a fQu^plloton 
generalization of the Fowler and DuBridge theories of photoemission. 


Since the advent of the laser nuiltiple-photon photo- 
emission has been observed from several metals. Two- 
photon photoemission from sodium has been observed 
by Teich and Wolga by the use of gallium arsenide 
laser [I] . Both two-photon and three-photon photo- 
emission have been observed from gold by the use of 
a ruby laser and the second harmonic radiation pro- 
duced by a ruby device [2] . The observation of 
higher order multiple-photon photoemission has now 
been facQitated by the use of higli power, ultrashort 
pulse, mode-locked lasers. These lasers can produce 
the very hi^ i*radiances necessary for the observation 
of multiple-photon photoemission, and yet because of 
the short pulse duration, catastrophic heating effects 
may no+ occur. Preliminary results of electron emis- 
sion produced by a mode-locked Nd:glass laser have 
been preseiiiedin ref. [3]. Unfortunately, these results 
were inconclusive because the authors of ref. [3] used 
the entire pulse train of the laser, and they found an 
asymmetry in the electron emission depending upon 
whether it was produced by the leading or the trailtag 
portions of the train of pulses. They attributed this 
asymmetry to the chirping of the pulses [4] or to sat- 
ellite pulse formation. We report here the observation 
of four-photon photoemission from tungsten by the 
use of carefully controlled and calibrated single, am- 

* Research supported in part by the National Aeronautics and 
Space Administration (under contract NGR'22-007-117) 
and by tiie Joint Services Electronics Frogicm (under con- 
tract N00014-67-A-0298-0006). 
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Fig. 1. Schematic diagram of tlie o;.peiimental arrangement. 
MLO, mode-locked oscHiator; PC , Pockels ceil; LTSG, laser- 
triggered spark gap: AMP, Nd:YAG amplifier; PD, photodiode; 
CRO, oscilloscope; KDP, potassium dihydrogen phosphate 
second harmonic generation crystal; F, filter; HV, hi^i voltage 
power supply; VT, electron emission vacuum tube. 

pUfled pulses produced by a passively mode-locked 
Nd:YAG laser followed by two Nd:YAG amplifiers. 

Fig. 1 exhibits a block diagram of the experimental 
apparatus. The laser oscillator used for this experiment 
is a 6 mm diameter by 3 inch long Nd:YAG rod. The 
laser is simultaneously Q-switcIied and mode-locked 
by a flowing saturable absorber, Kodak 9860, placed 
in contact with the loigh reflectivity mirror of the 
oscillator cavity. The laser oscillator output is restric- 
ted to the TEMqo mode by a 2 mm diameter aperture 
placed witliin the cavity. The output of the oscillator 
consists of a train of subnanosecond pulses separated 
by approximately 8 nsec. From tills train a single 
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pulse is selected for amplification by a spark-gap acti- 
vated, cylindrical-ring electrode Pockels cell. After 
this amplification with two single pass Nd:YAGampU- 
fiers we obtam a single, Iii^ power, subnanosecond 
pulse with a wavelength of 1 .064 pm. Tlie spatial 
profile of this pulse was measured by transverse scans 
of a 50 pm diameter pinhole. Tite spatial distribution 
of the laser irradiance at the site of the experiment 
was found to be gaussian with a radius at the 1 /e 
point of the irradiance equal to 1, 1 mm. 

Tire determination of the pulse duration of a mode- 
locked NdiYAG laser imposes some problems because 
the pulse width is less than the response time of the 
fastest available oscflloscopes. In order to determine 
the average pulse length we have used both a two- 
photon fluorescence technique [5] and a second har- 
monic nonlinear correlation method [ 6 ] . Both of 
these methods utilized many selected pulses, and the 
results of both measurements give an average pulse 
duration jf 30 psec full width at half maximum 
(fwhm) for tlie laser used in tliis experiment. 

To determine the incident peak irradiance of an 
individual pulse it is necessary to know the pulse 
energy, the pulse duration, the beam area, and the 
pulse sliape. We shall assume that the functional form 
of the pulse shape does not change from shot to shot, 
and that it is both spatially and temporally gaussian. 
The energy of an individual pulse was determined by 
the response of photodiode 3 (see fig. 1) which is cali- 
brated with an Eppley tliermopile. The product of the 
pulse duration and the beam area may also be deter- 
mined fci an individual pulse by a second harmonic 
generation method. For a pulse irradiance of the 
form 

I(r, 0 = /m exp [- (rld)^] exp [“(f/r)^] , 


the pulse energy, Ef, is given by the expression A’f = 
^3/2 For small conversion efficiencies the total 

energy at the second harmonic wavelength produced 
by this pulse is given by the expression 




m' 


Here >s the effective nonlinear susceptibility. From 
this we note that the ratio EflE^ is proportional to 
the product of the beam area and the pulse duration. 
Therefore by equating the average value of E^jE^^ for 
a series of laser pulses to the known product of the 
average beam area and the average pulse duration we 


can determine the product of the pulse width and the 
beam area for an individual pulse. Since we know the 
total pulse energy, we can calculate tlie maximum ir^ 
radiance, for each pulse. Moreover, if it can be 
assumed that the beam area does not change from 
shot to shot, then this technique yields single pulse 
temporal widths in a simple manner. This detailed 
knowi'.'dge is essential for an accurate quantitative 
mea'surement of a highh' nonlinear process. 

During the talcing of data for tlie laser-induced elec- 
tion emission the incident laser pulse was monitored 
with a fast biplanar photodiode (ITT F4000) and a 
Tektronix 519 oscilloscope (oscilliscope !). Tliis was 
done to insure that the oscillator did not have more 
than one pulse in the cavity during tlie mode-locking 
and that the Pockels cell was selecting only one pulse 
for amplification. The energy of tlie 1. 064 fim laser 
pulse was measured for each data point by the com- 
bination of oscilloscope 2 and photodiode 3, The 
second harmonic radiation was produced in a phase- 
matched KDP crystal, and the energy at 0.532 /urn 
was measured for each data point by the combination 
of oscilloscope 2 and photodiode 2. 

The electron emission was detected by measuring 
the total charge emitted from the tungsten ribbon 
target. This emission was recorded photographically 
with oscilloscope 3. The tungsten ribbon was 0.003 
inches thick, 0.04 inches wide and approximately 
1 inch long. Tlie pressure in the vacuum tube con- 
taining the tungsten ribbon was measured as 
1 X 10"'^ tori. The ribbon could be outgassed by pas- 
sing a current of approximately 10 A tlirou^i it, and 
immediately before tlie taking of data very lii^ ir- 
radiar.ces produced by the laser were used to further 
clean the area from which the emission would occur. 
The voltage between the cathode and tire anode was 
6000 V. 

In fig. 2 we plot the effective current density versus 
the maximum incident irradiance. By effective current 
density we mean the total emitted charge per pulse 
divided by the product of tlie laser area (ird^) and the 
pulse duration (fwhm). 

In order to explain the observed electron emission 
we shall use a theory that can be considered as a gen- 
eralization of the photoemission theories of Fowler 
[71 and DuBridge [8,9] . The four-photon current 
density is given by the expression 
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Fig. 2. Plot of ttic effective cutient density [total emitted 
cliargo divided by tiie product of nd* and the pulse duration 
(fwhm)] versus tlie peak incident irradiance. Tlie theoretical 
curve is obtained from a four-photon generalization of the 
theories of Fowler and DuBridge. 

Here e is the electron charge, R is the surface reflectiv- 
ity, A is the theoretical Richardson constant, hv is the 
energy of aa incident pfiuton, ^ is the work function, 
k is Boltzmann’s constant, /is the incident irradiance, 
ris the absolute temperature of the surface, and is 
an empirically determined constant. The function 
Fix) is tlie Fowler function and is tabulated in ref. [ 8 ] , 
nie total charge emhted per pulse is theoretically pre- 
dicted by a numerical integration of over the spatial 
and temporal dimensions of the incident pulse. To per- 
form this integration we assume that the incident irra- 
diance is spatially and temporally gaussian. Using this 
form for the incident irradiance we can determine the 


spatial and the temporal dependence of the surface 
temperature by solving the heat conduction equation. 
The details of this temperature calculation will appear 
elsewhere. We note, however, that the result of this 
calculation gives a maximum surface temperature 
change of approximately 700 K for a 1 GW/cm^ maxi- 
mum irradiance, 30 psec (fwhm) pulse impinging on a 
tungsten target. 

At tire higliest irradiances used in this experiment 
the maximum temperature is estimated to be near 
3000 K. The question may tlierefore be raised, whether 
excited contributions to the emitted current become 
important. Quite generally the current density may be 
written as 


«=0 

with 


'n » 




The term^o 's independent of the incident li^t in- 
tensity and represents pure tliermionic emission. At 
a maximum irradiance of 4 GW/cm^ we calculate 
that the tliermionic current density is less than 1 A/crn^. 
This is more than three orders of magnitude less tlian 
the measured effective current density for a maximum 
irradiaf'Ce of 4 GW/cm^, and tluis thermionic emission 
is negliijible. 

It is possible to have contributions from terms ana- 
logous to J 4 involving one, two, or three photons [ 10 ] , 
given by , J 2 and 73 . These terms represent photo- 
electric emission processes from the thermally excited 
high-energy tail of the Fermi- Dirac electron distribu- 
tion. For these processes tlie argument of tlie Fowler 
function is negative, as it is for/g- Their contribution 
is very small even for the hipest temperatures achieved 
here. Moreover, if there were any contribution from 
these lower order processes, then the slope of the ex- 
perimental curve would be at least a factor of two 
larger than the observed slope. 

The result of tlie theoretical calculation of the ef- 
fective four-photon current density is exhibited also 
in fig. 2. Tlie constant 04 is chosen to be 5 X 
cm^-sec^/C'^, and the pulse duration is assumed to be 
30 psec (fwhm). Moreover, we also assume that the 
reflectivity at 1.064 jum is 0.6, the work function is 
4.52 eV, and the photon energy is 1.165 eV. Wlien the 
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theoretical curve and the experimental data points are 
plotted on a log-log scale as in fig, 2, then tite slope of 
each should be nearly four if the dominant electron 
emission process is four-photon pbotoemission. The 
thermally induced contributions would increase the 
slope to values higher than four, at the large irradiances. 

In addition, since different pulse widths produce 
different temperature changes for a given peak itra- 
diance, some fluctuations in the effective current den- 
sity ate expected for a given peak irradiance. 

We conclude that careful control and careful calibra- 
tion of individual picosecond pulses have permitted 
the unambiguous verification of a four-photon photo- 
electric process. For longer pulses, more heating occurs 
for a given level of irradiance. The contributions from 
thermally induced processes become usually significant, 
and for longer pulses it is difficult to identify unambig- 
uously one high-order multiphoton process. 
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Two-photon absorption in semiconductors with picosecond laser pulses’'' 
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Single puUes from a well-calibrated modc-loclscd YAlGiNd laser have been used to measure the two-photon 
absorption, coefficient at 1.064 >im in several semiconductors. The materials studied are four direct-gap 
semiconductors, QaAs, CdTc, ZnTe, and CdSe, and one indirect-gap scraiconduclor, QaF. The results for the 
direct-gap semiconductors are interpreted witli respect to the imaginary part of the third-order nonlinear 
susceptibility anisotropy of the two-photon absorption coefficient is ob.scrvcd in GaAs. 


I. INTRODUCTION 

With the development of quantum mechanics in 
the late nineteen-twenties it became clear that 
when electromagnetic radiation interacts with mat- 
ter, processes involving more than one photon can 
occur. One such interaction is two-photon absorp- 
tion (TP A) and the basic theory for this process 
was formulated by Maria Goppert-Mayer in 1931.* 
The experimental observation of two-photon absorp- 
tion had to await the development of the laser, how- 
ever, and the first experimental observation was 
reported by ICalaer and Garrett in 1961.® 

Several techniques are used to study two-photon 
absorption. One technique involves measuring the 
transmission of the two-photon absorbing material 
for different incident ir radiances. Another in- 
volves measuring the luminescence from a materi- 
al pumped by two-photon excitation, and a third 
involves measuring the change in the conductivity 
of a material under two-photon excitation. 

There are several reasons for studying two-pho- 
ton absorption in various materials. Not only can 
two-photon absorption measurements yield new in- 
formation on the electronic states in a solid, but 
also TP A can provide a novel method of exciting 
large volumes in semiconducting and insulating 
materials. Such excitation has been used previous- 
ly to pump semiconductor lasers. ® Moreover, 
two-photon absorption can provide a damaging 
mechanism in materials that are nearly transpar- 
ent to low- Intensity radiation. This will constitute 
a very serious materials problem to the develop- 
ment of high-power ultraviolet lasers. Finally, 
two-photon absorption techniques in both liquids 
and solids have provided useful tools in ultra- 
short-pulse laser diagnostics.^ 

For comprehensive reviews involving many as- 
pects of two-photon absorption and spectroscopy we 
refer the reader to Refs. 5 and 6. 

There have been several previous studies of 
TPA in semiconductors, but there are very large 
differences in the measured values of the TPA co- 
efficients, For example, in GaAs the measured 
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values range from 0,02-5.6 cm/MW at a laser 
wavelength of 1.06 fim. Values of the TPA coeffi- 
cient will be presented here for GaAs, CdTe, 

ZnTe, CdSe, and GaP. 

Since two-photon absorption is a nonlinear opti- 
cal phenomenon, the TPA coefficient can be ex- 
pressed in terms of nonlinear optical susceptibili- 
ties, In Sec. H we outline a theory of TPA and re- 
late this theory to nonlinear optics, hi Sec. Ill the 
experimental apparatus and procedure are dis- 
cussed, In Sec. IV the results of the TPA mea- 
surements are presented. Finally in Sec. V we 
discuss the results of our experiments and their 
relation to other experiments and theory. 

11. THEORY 

In semiconductors it is possible to have TPA 
from a single light pulse with angular frequency w 
if twice the photon energy, 2 ^Tu 3, exceeds the semi- 
conductor energy gap Eg, The number of two -pho- 
ton transitions per unit volume per unit time IV‘®* 
is given by second-order perturbation theory as 

V - 2 “' . «> 

where 

1/ V" Pfj-oPtfq 

Here n is the index of refraction of the semicon- 
ductor, r is the vacuum speed of light, >u is tlie 
electron mass, e is tlie electron charge, h is 
Planck’s constant divided by 2a, and I is the power 
per unit area of the laser radiation. Tlie sums 
are to be extended over aU final states (or bands) 
denoted by the index / and all intermediate states 
denoted by the index t. The ground state is denoted 
by the subscript and the transition rale should 
be averaged over this index if there is more than 
one possible initial state. 

The unit vector a is directed along the electric 
field of the laser, and p/, and p„ are matrix ele- 
ments of the momentum operator between respec- 
tively the final slate and tlie intermediate state, 
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and the intermediate state and the ground state. 

The angular frequencies W/, and Wj, correspond to 
the angular frequency differences between the final 
and ground states and the intermediate and ground 
states, respectively. These di^erences, in gener- 
al, depend on the wave vector k. If we express 
tlie transition rate in terms of the peak electric 
field amplitude E instead of the irradiance I the 
transition rate per unit volume becomes 





- 2ft w) . 


( 2 ) 


In addition to the perturbation-theory calcula- 
tions of that involve matrix elements between 
tlie Initial and final unperturbed states, it is also 
possible to calculate the transition rate by a method 
that includes the effect of the electromagnetic field 
on the Bloch wave functions at the beginning of the 
calculation. Such a theory has been described by 
Keldysh. In this model the band distortion due to 
the eiecU'omagnetic field is included and the re- 
sults are applicable to both the optical and the dc 
limits. Moreover, for higher-order multiple- 
photon processes the Keldysh model offers both a 
simple and convenient method of estimating transi- 
tion rates. 

One method of measuring TP A in semiconductors 
is to measure the attenuatlonofabeamorpulseprop- 
agating through a two -photon absorbing medium. 

If the pulse irradiance is /, tlien the pulse is at- 
tenuated according to the expression 


(iz 




(3) 


Here a is the one-photon absorption coefficient and 
d is the two-photon absorption coefficient. If we 
neglect multiple reflections within the medium, the 
transmitted irradiance is given by the expression 


l{r, t, t) =• 


1 + da - fi) Ar, 0, mi - 


(4) 


Here I is the thickness and R is the reflectivity of 
the medium. 

The TPA coefficient d is related to the two-pho- 
ton transition rate per unit volume by the expres- 
sion 


0 = . (5) 

An alternative formulation of the TPA problem 
is to start with Maxwell’s wave equation and treat 
TPA as a nonlinear source polarization. For a 
monochromatic beam the electric field amplitude 
E(r, tij) is determined by the solution to this wave 
equation given by 

VxvxE(r, w)-((//c®)e .E(f, w) 


= (4ntoVo®) . (6) 

Here ? is the linear dielectric tensor and P°*®(r, a. ) 
is the nonlinear source dipole moment per unit 
volume. 

For TPA the relevant nonlinear source polariza- 
tion term is cubic in the electric field anti is pure- 
ly imaginary in character. In centrosymmetrlc 
media the relation is 

P^'*(r, tj) w, - a#) 

. (7) 

Here Xi “J, «.>, w - w) is the imaginary part of 
tlie complex third-order nonlinear susceptibility. 

In non-centro-symmetric media for electric field 
directions such that second-harmonic generation 
is not symmetry forbidden there are contributions 
to not only from the imaginary part of 
x“'{— w, ai, w, - w) but also from terms propor- 
tional to a.', 2w, w) 2w, w, w}. For a de- 
tailed discussion of these effects see Ref, 8. 

We will be interested in determining the TPA 
coefficient d for laser electric field directions 
that preclude the generation of second-harmonic 
radiation. In this case the total nonlinear source 
polarization is given bv Eq. (71. By using 
the slowly varying amplitude approximation to 
solve Eq, (6) we can express 0 in terms of 

a-’, a?, - u). The explicit relation between 
d and x/3i{'(- a', u>, (ij - , 1 )) depends both on the crys- 
tal class and the laser electric field direction. 

For example, in a crystal of symmetry (e.g., 
GaAs, CdTe, or ZnTe) with the laser electric field 
polarized along the [OOl] direction we find 

d = (327f®u/n^r®) [3xiiir(- w, a', co, - w)l . (8) 

Here w is the laser angular frequency, c is the 
vacuum speed of light, and n is the index of refrac- 
tion at the laser frequency. Note that we are using 
the Maker and Terhune® convention for the third- 
order nonlinear susceptibility. In this convention 
the effective nonlinear susceptibility is 
3xiii['(- w;co, u, - w) because there are three dis- 
tinct permutations of the positive and negative fre- 
quency arguments. This is the same convention 
that has recently been used by Levenson and 
Bloembergen'^'’ in their study of the resonant be- 
havior of third-order nonlinear susceptibilities 
by three-wave mixing. 

The results that we have heretofore discussed 
are for pure two-photon transitions. These theo- 
retical results are appropriate for the interpreta- 
tion of data for which the direct energy gap of the 
semiconductor is less than the combined energy of 
two photons, It is also possible to have TPA in an 
indirect gap semiconductor by the simultaneous 
emission or absorption of a phonon. Although we 
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report here TP A from GaP, an indirect gap mate- 
rial, we shall not attempt to present here the the- 
ory associated with this tliird-order process. For 
a discussion of the theoretical aspects of TP A 
from indirect gap materials the reader is referred 
to Refs. 11 and 12, 

There are, however, several theoretical aspects 
of TP A from indirect gap materials that are note- 
wortliy. Firstly, the TP A coefficient 3 is expected 
to be a strong function of the temperature because 
of the phonon participation in the TP A process. 
Secondly, for the laser intensities and frequency 
used in this experiment the TPA in indirect gap 
materials is less than that in direct gap materials. 

Until now we have tacitly assumed that absorp- 
tion by two-photon created excess carriers is 
negligible. This assumption is not always valid. 
Usually the dominant absorption by excess carriers 
comes from free holes; therefore, we can write by 
analogy with Eq. (3) 

= , (fl) 


where «(r, t) is the number of two-photon excited 
free holes per unit volume, and tj), is the free hole 
absorption cross section. The excess hole density 
«(?,/) is determined by the relation 


8»(r, /) 
a/ 


=^D^V^n{r,t) 


>t(r, t) t) 

/a ’ 


( 10 ) 


where is the hole diffusion coefficient and if, is 
the hole recombination lifetime. In general, one 
must solve the coupled Eqs, (9) and (XO). For 
picosecond pulses both diffusion and recombination 
may be neglected; therefore, the hole density is 
expressed by the relation 

«(?,/) . (U) 


Note also that we can solve Eq. (9) for the picosec- 
ond laser pulse if we assume that all absorption 
mechanisms are small. In this case the perturba- 
tion solution for the irradiance is given by the ex- 
pression 

lir, 2 , #) = (1 - 0 , t) - 0 , 1) - 0 , 0 

o-„3r/nCnO./){l-R)2 

XU-/Z) 

xjf' (12) 

We have assumed that the incident irradiance is 
/a(r, 0, /) and that the irradiated medium occupies 
the half-space 2 ^ 0 . 

Finally we note that fox picosecond pulses the 
attenuation due to the excited hole absorption is 
comparable to the direct two-photon absorption if 
the irradiance is equal to a critical irradiance 


given by the expression 

where t is the l/e pulse intensity half duration. 

The corresponding expression for pulses longer 
than the hole lifetime is given by a similar expres- 
sion in which the pulse duration is replaced by the 
hole lifetime. We conclude, therefore, that in 
order to obtain pulse attenuation due to pure two- 
photon absorption one should use the shortest pos- 
sible pulse duration. 

Ill, EXPERIMENTAL APPARATUS AND TECHNIQUES 
A. Laser system 

A passively mode-locked YAIG: Kd (neodymium- 
doped yttrium aluminum garnet) laser system pro- 
duced the picosecond pulses used in these eiqjeri- 
raents. A schematic diagram of the experimental 
apparatus is exhibited In Fig. 1. The laser oscil- 
lator is a xenon flashlamp pumped 6-mm-diameter 
by 7.62-cm-long YAIG: Nd rod In the Brewster- 
Brewster configuration. The oscillator is simul- 
taneously Q-switched and mode-locked by a flow- 
ing saturable absorber (Kodak 9860) in contact with 
the tigh reflectivity mirror of the oscillator cav- 
ity. The output of the oscillator is restricted to 
the TEMaa mode by a 2 mm diameter aperture. 

The output of the oscillator consists of a train of 
picosecond pulses separated by approximately 8 
nsec. From this train of pulses a single pulse is 
selected by a cylindrical-ring electrode potas- 
sium-dideuterium-phosphate Pockels cell. This 
Pockels cell is positioned between crossed Gian 
prisms, and is activated by a laser-triggered 
spark gap. The selected laser pulse is subse- 
quently amplified by two YAIG; Nd laser ampli- 
fiers. 



S 


FIG. 1. Schematic diagram of the experimental appa- 
ratus. YAlOiNtlMLO, neodymium doped YAIG mode-lncl:ed 
laser; PC, Pockels cell; LTSG, laser triggered spark 
gap; AMP, neodymium doped YAIG laser amplllior; PD, 
phofudlodc; CRo, oscilloscppe, F, filter; KDP, second- 
harmonic generation crystal; s, sample. 













I 


Arv-10 

35ie J. H. BECHTEL AND W. L. SMITH 13 


A high speed biplanar photodiode PD-1 (ITT 
F4000) is used in conjunction with CRO-1, a Tek- 
tronix 519 oscilloscope, to insure that only a sin- 
gle pulse is selected from the oscillator train of 
pulses. The calibrated photodiode PD- 3 (S-1 re- 
sponse) measures the energy incident on the two- 
photon absorbing medium, and photodiode PD-4 
(S-1 response) measures the energy transmission 
of this medium. A phase matched KOP (potassium 
dihydrogen phospiiate) second-harmonic generation 
crystal and photodiode PD-2 (s-20 response) are 
used together witii PD- 3 to measure the individual 
pulse maximum irradiance. This technique is dis- 
cussed below. 

B. Puke characteristics 

For an accurate quantitative measurement of 
many nonlinear optical properties of materials it 
is necessary to know the spatial and the temporal 
form of the pulse incident on the material being 
studied. The spatial profile of the laser pulse at 
the site of the hvo-photon absorbing medium was 
determined by several independent techniques. 

In one experiment we employed transverse scans 
of a 50 jum diameter pinhole coupled to a photo- 
diode. The results of one of these scans are shown 
in Fig. 2. In another experiment we directly mon- 
itored with an oscilloscope tiie output of a linear 
array of 256 silicon photodiodes that were irradiated 
by an attenuated laser pulse. The results of each 
of these experiments indicated that the spatial form 
of the incident irradiance was Gaussian, and that 
the 1/e intensity radius was 1.1 ±0.1 mm. 

In order to determine the temporal pulse width 
two techniques were utilized— a second-harmonic 
autocorrelation (SHAC) method and two-photon 
fluorescence (TPF) photography. In the TPF 



DISTANCE FROM PULSE- CENTER (mml 

FIG. 2. Irradiance distribution at the sample. Solid 
line is a Gaussian fit to the data yielding a 1/c irradiance 
radius of 1. 1 mm. 



FIG. 3. Results of measurement of the average temporal 
duration of the laser pulse using a second-harmonic auto- 
correlation technique. Solid curve Is a least-squares fit 
to tlic data yielding a FWHM duration of 30 iisec and a 
contrast ratio of 2. 7 to 1. 

technique the initial pulse is split into two pulses 
by a beamsplitter, and these pulses are redirected 
by mirrors to overlap within a two-photon absorb- 
ing material such as Rhc-damine 6G. The dye fluo- 
rescence is photo ..aphed, and the pulse duration 
is determined by the spatial distribution of the ex- 
posure produced by the fluorescence track. In the 
second-iiarmonic autooirrelation method the pulse 
is split into two pulses by a Michelson interferom- 
eter. Tiiese two pulses are collinearly recom- 
bined after one pulse has been delayed by a vari- 
able amount. The recombined pulses are subse- 
quently directed through a phase matched second- 
harmonic generation crystal. One plots the sec- 
ond-harmonic output as a function of the delay, 
and the result of such a plot is exhibited in Fig. 3. 
The curve is a four parameter least- squares fit 
to the function, /Ci +K-> e?^[- 2{s - iro/2)®/c^T^], 

This curve gives a contrast ratio, [Ki + 
of 2. 7 to 1. The results of both the TPF photog- 
raphy and the SHAC method give the average tem- 
poral duration full width at half maximum of 30 
±6 psec, implying an average 1/e halfwldth (t) of 
18±3.6 psec. 

Although single-pulse TPF photographs are pos- 
sible, the SHAC method is based on using many 
laser shots. We have found, however, that there 
are pulse -to-pulse fluctuations in the pulse dura- 
tion. A convenient way to measure the pulse dura- 
tion, and thus the pulse irradiance for individual 
pulses, is to employ a nonlinear technique.^® This 
method was first suggested by Glenn and Brienza*® 
and involves measuring both the incident pulse en- 
ergy at tlie fundamental frequency and the total 
energy created at the second-harmonic fre- 
quency from d portion of the incident pulse. For 
small second-liarmonic conversion efficiencies in 
the phase matched second-harmonic generation 
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crystal the relation ^i/iSah is proportional to the 
product of the Incident pulse duration and the inci- 
dent pulse area. If in an experiment a sufficiently 
large number of data points are taken so that the 
ratio is meaningful, tlien this average value 

con be equated to the product of the known values 
of die average pulse duration and the average pulse 
area. Thus we can measure the individual pulse 
maximum irradiance on a shot-to-shot basis, 

C. Procedure 

The technique of obtaining the two-photon ab- 
sorption coefficient 0 is to measure for each laser 
pulse the incident pulse energy, the transmitted 
pulse energy, and the pulse energy created at the 
second-harmonic frequency. The measurement of 
the pulse energy at the laser frequency in conjunc- 
tion with the pulse energy at the second-harmonic 
frequency allows us to know the maximum pulse 
irradiance for each laser pulse as described in the 
previous section. The reciprocal of the energy 
transmission is then plotted versus the maximum 
incident irradiance. The incident irradiance was 
changed by either changing the gain of the amplifiers 
or by changing the position of the pulse selected 
from the oscillator train. 

iV. EXPERIMENTAL RESULTS 


In order to determine p, the two-photon absorp- 
tion coeffi:ient, we measured the energy transmis- 
sion for different irradiances. The irradiance exit- 
ing a sample of Ihiclcness I for a spatially and tem- 
porally Gaussian laser pulse is given by the ex- 
pression 


Hr, I,/) 

“ 1 +^ 4(1 - R ) (1 - ?-“')/« • ' 


The energy transmission coefficient T is thus 
T - jTrrrdrCUrJ, t)dt 

„ tt e*“' 1-R 

X 1“ In^l +1 4(1 - R)(l - e -“') dx . (14) 

Expression (14) represents the basic equation used 
to derive p from the e:q)erlmental measurement of 
T (or T‘^) vs 4. Of course Eq. (14) is only ap- 
propriate when the absorption by two-photon ex- 
cited free holes in negligible. The limitations thus 
required will be considered separately for each "f 
the materials studied. 

The TPA coefficient p was measured for three 
samples of single-crystal n-type GaAs with the 
laser electric field along the [OOl) direction of this 



FIG. 4. Plot of reciprocal cnei’isy transmission coef- 
ficient vs niaxlmum incident Irradiance for GaAs. Sam- 
ple thickness Is 0.22 cm, and the laser electric field 
was polarized along the [0011 direction. Theoretical 
curve is given by Ect, (14) with q = 0. 7 cm’' and fJ= 0. 03 
cm/MW. 


zinc-blende material. In aH experiments unless 
otherwise indicated tlie samples were Irradiated 
at normal incidence along tiie [llO] direction. 
Sample I was Te doped wiUi a carrier concentra- 
tion of 2xlo“ cm“"; the thicloiess was 0.09 cm. 

The experimental results indicate that 41=0.030 
±0,005 cm/MW with a one-photon absorption coef- 
ficient a =8.0 cm*'. Samples 2 and 3 were both 
cut from the same boule and boUi had a carrier 
concentration of 5x10'“ cm’“. Sample 2 had a 
length of 0, 22 cm, and j3 was measured to be 
0.030 ±0.005 cm/MW with or = 0.7 cm*'. The re- 
sults for this sample are plotted in Fig. 4, From 
this figure we note that for incident irradiances 
above 300 MW/cm“ there is a small but significant 
departure of the experimental data points from the 
theoretical curve. Similar effects were observed 
in other samples. This departure can be attributed 
to the free hole absorption, and thus there Is an 
enhanced attenuation of the pulse at higli intensi- 
ties, AnoUier possible explanation is beam distor- 
tion due to the real part of the third-order nonlin- 
ear susceptibility. This will be considered in 
Sec. V. 

Sample 3 had a length of 1.20 cm; the one-photon 
absorption coefficient was the same as Sample 2, 
but the TPA coefficient was found to be 0.023 
±0.005 cm/MW. The reason that the measured 
value of 43 is smaller for the longer sample prob- 
ably results from the reflection at the exit surface. 
The reflected wave suffers a n phase change that 
enhances the total irradiance at and near the exit 
surface. Since the model that we are using doesn’t 
include the effects of multiple reflections, we can 
expect the effective TPA coefficients to be slightly 
smaller for longer crystals. 
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TABLE I. Various experimental values of the hvo- 
plmton absorption coefficient p, the linear absoriJtlon 
coefficient tv, and the optical electric field direction in 
GaAs at !• 06 (im and at room temperature. 



J3tcn^f'^aV) 

o (em'h 

Electric 

Hold 

PrcHcnt wort; 

V.02tl-U.040 

0,7-a.O 

tool] 

Geplio* 

D.oas 

* * , 

[0011 

Lcd and Fan“ 

0.23-0,GG 

4,0 


Kteinman at cl. ° 

0.010-U.045 

0.05-4.0 

« • • 

Jayaramon anil Lee'’ 

5,6 

a-4 

• • • 

Olasman ct at." 

C.2 

... 

• ■ • 

nalnUm and Chang’ 

0,02 

... 

• ■ « 

Arsen'ov oJ.' 

H).a 

4.0 


Gasov ct af.“ 

5.0 

... 



“Reference 18. 
“Reference 20, 
“Reference 21. 
‘^Reference 22. 


“Reference 23, 
'Reference 24, 
“Reference 25, 
“Reference 26. 


used was 4 mm thick with the [111] direction nor- 
mal to the optically polished surfaces. For our 
experiment data were taken with the laser electric 
field along the [llOj direction. 

In order to measure any change in the transmis- 
sion coefficient it was necessary to use incident 
irradiances in excess of 1 GW/cm®, Since this 
irradiance is comparable to tlie critical irradiance 
4[. for GaP then the effect of free hole absorption 
must be considered. This was considered in our 
interpretation of the data, and the TP A coefficient 
^ was determined by comparing the experimental 
energy transmission coefficient to the theoretical 
coefficient obtained by integrating Eq. (12) over 
the spatial and temporal profile of the laser pulse. 
The free hole absorption cross section was as- 
sumed to be 2x10*'’' From our data we find 

that cv=0.5 cm"' and d =^0.0002 ±0,0001 cm/MW. 

V. DISCUSSION AND SUMMARY 


Finally we observed a small anisotropy of the 
TPA coefficient in Sample 3. For the laser field 
along the [llO] direction we find that ^ is approxi- 
mately 20^ larger than when the laser field is 
along the [OOl] direction. A similar anisotropy 
was recently observed In Ref. 18. 

A single crystal of high-resistivlty (o > 10“ D cm) 
«-type CdTe was examined for TPA. The crystal 
was 0.49 cm in length and was irradiated along 
the [llO] direction with the electric field polarized 
along the fOOlj axis. The best fit to the experi- 
mental data gives a -0.45 cm”' and j3 =0.025 
±0,005 cm/MW. 

A sample of O.lO-cra-thick />-type ZnTe was ir- 
radiated along the [llO] direction with the laser 
electric field along the [001] direction. The best 
fit to the data gives <v =8 cm”' and ^ =0.008 
±0.004 cm/MW. 

In addition to the zinc-blende crystals described 
above we have also measured the TPA properties 
of CdSe, a wurtzite structure. The sample was 
0.385 cm thick, and the incident laser irradiance 
was directed along the r axis of tins hexagonal 
crystal. The data indicate that cv =0.2 cm*' and 
(3 = 0.030 ±0.005 cm/MW, 

Of the materials considered in this experiment 
GaP is unique for two reasons. Firstly, at the 
laser energy of 1.165 eV per photon, GaP does 
not exhibit a direct-gap two-photon absorption. 
Instead U.e transition needs a phonon to participate 
in the absorption process. Secondly, the only 
crystal Uiat wo had available to us had an orienta- 
tion such that the laser electric field could not be 
polarized along the [001] direction. We found that 
the two-photon absorption coefficient is much 
smaller in tliis material than in any of the other 
materials that we studied. The sample that we 


There are very large discrepancies among the 
various measurements of the two-photon absorption 
coefficient J3 in GaAs as can be seen from Table 
I. Our measurements seem to be in closest 

agreement with those of Ref, 21. It is noteworthy 
that any spatial or temporal fluctuations in the 
laser pulses used for the TPA measurements can 
lead to larger TPA. Thus If the lasers used in 
some of these other experiments did not produce 
smooth reproducible pulses, the inferred value of 
the two-photon absorption coefficient would be 
larger than the true value. 

In Sec. II we discussed the importance of free 
hole absorption and noted that picosecond pulses 
are favored over nanosecond pulses for an accurate 
measurement of the TPA coefficient. The cross 
sections for free hole absorption were estimated 
from the one-photon absorption measurements in 
/(-type materials. The results are as follows: 
1x10”''^ cm“ (GaAs),^^ 4x10”'*' cm^ (CdTe),^“ 
3x10*'’' cm® (ZnTe),®“ and 2x10”''' cm“ (GaP).'° 
There is no known measurement for CdSe. From 
these values the critical irradiance, la, was cal- 
culated and except for GaP, as noted above, only 
data points less than were used to infer the TPA 
coefficient. If we compare our experiments wivli 
others as in Tables I and II we find that our experi- 
mental results for the TP A coefficients are among 
the smallest reported values. Our results for the 
direct-gap semiconductors are typically within a 
factor of 2 of the theoretical resu' ,s obtained using 
the Keldysh model to calculate the TPA coefficient. 

Our experimental results as well as the results 
based on several theoretical calculations are pre- 
sented in Table III, For the direct-gap semicon- 
ductors the TPA coefficient was measured with the 
laser electric field polarized along tiie [OOl] direc- 
tion; thus, second-harmonic generation was pro- 
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TAI3LE II. ENTjcrimcntal values of Iho ftvo-photon 
absorption cnofficlcnt 0 in CdTe, CdSe, ZnTe, and Gal- 
at room tomperaturo. The units arc cm/J.lW. 



CdTe 

CdSo 

ZnTe 

Gap 

Present work 

0.02S 

0.03 

0,008 

2x10''’ 

Bcpko“ 

0.17 




Ralston and Chang” 

0.2 

0.2 



Basov ei al,' 


O.DS 



Lee and Fan'’ 



0.04 


Yee and Chau' 




1.7x10-3 

Bryulmer ct al' 


0.14 




“Reforoncc 18, 

''Reference 24. 

''Reference 20. 

'’Reference 20. 

®Reforonoo 30, 

*F. Bryukner, V. S, DnoprovslUi, and V. U. Khal- 
tatov, Kvant. Electron. ^ 1300 (1074)5 [Sov. J. Quan- 
tum Electron, 4, 740 (1074)]. 


hibited for tliese materials, and the imaginary part 
of the third-order optical nonlinear susceptibility 
Xiui(“ o), tu, - w) at the laser frequency may be 
inferred from the TP A coefficient. 

In Gap several different attenuation mechanisms 
are possible. In addition to indirect TP A and free 
hole absorption it is also possible to have second- 
harmonic generation followed by one-photon ab- 
sorption for the laser electric field polarized along 
the [llO] direction. Finally since the observed at- 
tenuation was very small three-photon absorption 
was also considered. To predict the amount of 
three-photon absorption a simple calculation using 
the Keldysh model shows that three-photon absorp- 
tion is negligible compared to the other attenuation 
mechanisms at the irradiances used in this experi- 
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mont. Moreover, wc can calculate tho effective 
TP A coefficient due entirely to second harmonic 
generation followed by one-photon absorption. We 
find for GaP at 1,064 pm that j3sno = 10*° cm/MW. 
This is approximately two orders of magnitude less 
than, the experimental value that we obtain. We 
conclude, therefore, that the attenuation mechanism 
is a combination of hvo-photon absorption across 
the indirect gap plus linear absorption by the re- 
sulting free holes. If the mlue of the free hole ab- 
sorption cross section is known, then we may infer 
the TP A coefficient 0 from the measured transmis- 
sion change. This value is reported also in Table 
HI. Our experimental value compares very favor- 
ably with the theoretical value given in Ref. 30, 
but is approximately an order of magnitude less 
than the experimental value reported in this same 
reference. 

In Sec. rv we noted that at high irradiances there 
is a departure of the e:q}crimental data from the 
theoretical curve ob>ained from the lower irradi- 
ance data. We noted that one possible explanation 
for this departure would be beam deformation due 
to the self -focusing aspect of tlie real part of the 
third-order nonlinear susceptibility, hi order to 
estimate this effect we have calculated the beam 
intensity radius n(f) at the exit surface of the crys- 
tal for an input radius Wq* The relation between 
these radii is given by the foHowing approxima- 
tion®^ 



The length of the crystal is I, P is the maximum 
laser power, and fto is the laser v;ave vector in 
vacuum. The parameter P„ is the critical power 
for self-focusing and is given by 

P,r=r®/8w®H2 . (16) 

Here c is the vacuum speed of light, ft' is the laser 


TABLE III. Summary of the experimental results obtained in this work, and the results of several tijcoretical models. 


Material 

Index of 
refraction 

0 (cm/MW)“ 

(cm^/erg) 

1” 

j3 (cm/MVV), 
2' 

The-iry 

v'’ 

4' 

GaAs 

3.475 

(2,8±0.5)x10-3 

l.exior” 

1.7x10"* 

3.1x10** 

0.2x10** 


CdTe 

2.818 

(2. 5*0, 5) X 10-3 

l.lxlo'” 

2. lx 10"* 




CdSe 

2.539 

(3,0*O.6)x10-3 

1,'Oxio-” 

1.9x10"* 




ZnTe 

2.777 

(8. 0*4.0) X 10*3 

3.3x10-** 

a.sxitr* 

2.5x10-* 

3.2x10"* 


GaP 

3.106 

(2*l)xl0-^ 





1.9x10"^ 


“If more tlian one sample was measured, 
”Calculated from Eq. (41) of Ref. 7. 
'Reference 20, degenerate valance band, 
'’Reference 20, degenerate valance band, 
'Reference 30. 


the average value Is given. 

no exciton. 
with exciton. 
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angular frequency, and Wg is the intensity depen- 
dent nonlinear index of refraction. For seK-focus- 
ing produced by an electronic nonlinearity, the in- 
tensity-dependent nonlinear index of refraction may 
he estimated by Wang’s rule.^^ This is the contri- 
bution due to tlie bound electrons, and as an ex- 
ample, for GaAs we estimate that «2 1.6x10’^° 
esu and F„aG70 W. Additional contributions to 
«2 may come from the nonparabolicity of the con- 
duction band and the number of conduction-band 
electrons,’^'®* However, for the intensities at 
which these effects become important, there is a 
decrease in the linear index of refraction tliat is 
linear in the number of excited free electrons. The 
net result would be a defocusing of the beam. Us- 
ing the maximum possible values of the effective 
;;2 and Eqs. (15) and (16) we find that the area of 
the laser pulse at the exit surface is less than 10% 
smaller than at the incident surface. This is true 
even for the longest sample at the highest power. 
Moreover, we have estimated the effects of ther- 
mal self -focusing due to the temperature depen- 
dence of the linear index of refraction. We find 
that these effects are also small. Finally, the 
beam size was observed using an ii Jrared-absorb- 
ing-vislble-emitting nhosphor. At no time was 
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there a discernible change in the beam size at 
either high or low powers. We conclude, therefore, 
that pulse deformation due to self-focusing effects 
was not a serious problei in the interpretation of 
any of the experiments reported here. 

In conclusion, we have measured tlie TP A coef- 
ficient for four direct gap and one indirect gap 
semiconductors. Our experimentally measux’ed 
values of the TP A coefficient are significantly less 
than most previous measurements, but in reason- 
able agreement with several tlieoretical models. 

For the direct gap semiconductors we have related 
our measurement of the TPA coefficient to tlie 
imaginary part of the third-order nonlinear suscep- 
tibility of the material. 
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